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•  Turbulence  frequency  cuioir  (.«"*) 
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-  Nusselt  number  on  the  turbine  blade  surface 
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3  •  X-tvpe  hot-wire  bisector  angle 

±  •  Offset  angle  of  hot-wire  bisector  from  the  tunnel  reference  direction 

•  Turbulence  frequency  bandwidth  (s“l) 

±1  •  Time  interval  {  #  ) 

±u,  •  Thermal  boundary  layer  thickness  (m) 

£3  •  Knthalpy  thickness  of  a  boundary  layer  (m) 

A*  »  Conduction  thickness  (m) 

Subscripts 

Unless  otherwise  noted,  subscripts  arc  designated  as  follows: 

i  -  designates  the  i  th  blade  surface  section 

z  •  pertains  to  no  jet-grid  in  the  flow 

0  -  Station  0  in  the  tunnel,  y  =  15.75 

02  -  Station  02  in  the  tunnel,  f  —  25.G7 

1  -  pertains  to  wire  1 

2  -  pertains  to  wire  2 

cc  -  pertains  to  the  free-stream 
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Abstract 

Heal  transfer  for  *  turbine  cascade  I*  examined  for  turbulence  scale  effects.  The 
turbulence  integral  scale  and  microscale  lengths  of  the  free-stream  flow  are  controlled  by 
air*jet  injection  through  a  grid  placed  in  the  free-stream  flow.  Air  is  injected  into  the 
flow  in  three  primary  directions:  co-flow,  cross-flow,  and  counter-flow,  at  several  Injection 
pressures.  Results  arc  obtained  for  heat  transfer  on  the  blade  without  the  jet-grid  in  the 
free-stream  Row.  Surface  pressures  on  the  turbine  blade  arc  examined  for  three  injection 
directions.  Surface  beat  transfer  on  the  blade  and  turbulence  intensity  and  turbulence 
scale  results  are  obtained  for  two  grid  locations,  three  air  injection  orientations,  and  three 
injection  orifice  diameters  over  a  range  of  injection  pressures.  Turbulence  Integral  scale 
and  microscale  lengths  are  measured  in  two  free-stream  Row  coordinates.  The  turbulence 
intensity  decay  is  determined.  Heat  transfer  on  the  turbine  blade  is  compared  with  the 
turbulence  integral  scale  and  microscale  lengths. 

Results  show  the  integral  scale  length  is  primarily  dependent  on  the  diameter  of  the 
jet-grid  tube,  but  may  be  controlled  to  a  lesser  degree  by  changing  the  orientation  of  the 
jet-grid  injection  to  a  co-flow  or  counter-flow  direction;  or  by  changing  the  secondary  air- 
jet  injection  pressure.  The  turbulence  microscale  is  primarily  dependent  on  the  location 
of  the  jet-grid  within  the  flow,  hut,  to  a  lesser  degree,  controlled  by  the  orientation  of  the 
jet-grid  injection,  jet  orifice  diameter,  ami  injection  pressure. 

In  general,  jet-grid  injection  produces  a  lower  turbulence  intensity  of  the  free-stream 
than  a  non-blowing  grid,  larger  turbulence  microscale,  and  lower  surface  beat  transfer  co¬ 
efficient.  The  turbulence  integral  scale  is  nearly  constant  with  turbulence  intensity  and 
its  length  is  on  the  order  of  the  outside  diameter  if  the  jet-grid  tube.  The  Stanton  num¬ 
ber  and  the  Nusselt  number  increase  with  increasing  turbulence  intensity  and  decreasing 
turbulence  microscale  length.  Turbulent  boundary  layer  heat  transfer  predictions  compare 
favorably  on  the  pressure  surface  of  tiie  turbine  blade,  except  at  the  leading  edge  of  the 
turbine  blade  for  moderate  free-stream  turbulence  levels  (<I-G  percent).  Laminar  boundary 
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layer  heat  transfer  predictions  compare  favorably  on  the  turbine  blade  suction  surface  for 
low  free-stream  turbulence  levels  (0-1.3  percent). 
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TURBULENCE  SCALE  EFFECTS  ON  HEAT  TRANSFER 
IN  A  LINEAR  TURBINE  CASCADE 


/.  INTRODUCTION 


1.1  General 

The  accurate  knowledge  of  heat  transfer  distribution  on  gas  turbine  blades  is  of 
practical  interest  in  the  design  of  turbotnachincry.  Accurate  prediction  of  convective  heat 
transfer  within  the  turbine  blade  cascade  is  of  particular  importance  to  the  design  of 
gas  turbine  airfoil  design,  which,  in  turn,  affects  the  overall  engine  cooling  design,  cycle 
efficiency,  and  hardware  durability.  Documentation  of  heat  transfer  data  for  turbine  blades 
is  relatively  scarce  as  compared  to  similar  documentation  available  for  flat  plate  geometries. 
Documentation  of  turbulence  scales  in  relation  to  turbine  blade  heat  transfer  is  nonexistent. 

Heat  transfer  within  a  gas  turbine  blade  cascade  has  been  studied  for  turbulence 
intensity  and  Reynolds  number  effects  (Consigny,  1082),  (Priddy,  1985).  Other  studies 
have  concentrated  on  turbulent  boundary  layer  development  over  the  turbine  blade  surface 
and  related  heal  transfer  to  the  blade  (Blair,  1989a),  (Priddy,  1985).  Blair  (Blair,  19S9b) 
reports  a  20-30  percent  increase  in  effective  heat  transfer  rates  for  high  turbulence  level 
(8-10  percent)  free-stream  flow  over  low  turbulence  frec-stream  flow.  More  recently,  heat 
transfer  on  a  flat  plate  for  various  turbulence  levels  generated  by  “jet-grid”  injection  of  a 
secondary  flow  into  the  primary  flow  has  been  studied  (Young,  1989).  Results  indicate  an 
increase  of  up  to  45  percent  in  Stanton  numbers  may  be  achieved  with  high  turbulence 
free-stream  flow  as  compared  to  low  turbulence  frec-stream  flow.  Additionally,  the  same 
study  found  the  Stanton  number,  for  the  flat  plate,  to  be  inversely  proportional  to  the 
turbulence  microscale. 

This  thesis  continues  in  the  manner  of  the  latter  study  for  a  turbine  cascade.  Specif¬ 
ically,  turbulence  integral  scales,  microscales,  turbulence  intensity  and  their  relation  to 
surface  heat  transfer  is  investigated  for  a  turbine  blade  cascade. 
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1.2  Objectives 

The  objectives  of  this  thesis  are  as  follows: 

(1)  Determine  the  turbulence  level  induced  by  jet-grid  injection  at  two  downstream 

* 

locations.  Determine  the  effects  of  injection  pressure,  injection  orifice  diameter,  and  injec¬ 
tion  direction  of  a  secondary  flow  on  frcc-strcam  turbulence. 

(2)  Measure  turbulence  energy,  integral  scales  and  microscales  for  each  (low  condition 
imposed  on  the  frcc-strcam  flow  by  jet-grid  injection.  Determine  the  relationship  between 
turbulence  intensity  and  turbulence  integral  scale  and  microscale  lengths  for  the  jet-grid 
configurations  given  in  objective  1. 

(3)  Measure  the  surface  pressures  on  the  turbine  blade  for  the  various  jet-grid  injec¬ 
tion  configurations. 

(•1)  Determine  the  heal  transfer  on  the  turbine  blade  in  terms  of  the  Stanton  and 
Nusselt  numbers  for  each  jet-grid  injection  configuration  given  in  objective  1.  Determine 
the  relationship  between  surface  heat  transfer  and  turbulence  scales. 

1.3  Method 

A  brief  overview  of  the  experimental  methods  used  is  now  presented. 

Turbine  blade  surface  pressures,  and  heat  transfer  were  measured  for  various  tur¬ 
bulence  levels.  FYce-sircam  turbulence  levels  were  controlled  by  varying  the  location  of 
the  jet-grid,  the  jet-grid  injection  pressure,  injection  orifice  diameter,  and  the  direction 
of  injection  with  respect  tc  the  mainsteam  flow.  A  total  of  60  different  frcc-strcam  (low 
conditions  were  tested,  along  with  3  baseline  free-stream  flow  configurations.  Each  test 
configuration  is  identified  by  a  three  place  letter,  number  or  combination  thereof.  The 
free-stream  flow  conditions  tested  arc  listed  in  Table  1.1  and  Table  1.2.  Table  1.1  is  a  list 
of  flow  configurations  tested  with  the  turbulence  jet-grid  placed  at  the  first  station  (des¬ 
ignated  as  station  0)  ahead  of  the  turbine  blade  row.  Its  upstremn  location  is  denoted  by 
the  nondiinensional  distance  where  i  is  a  linear  distance  and  b  is  the  outside  diameter 
of  the  turbulence  jet-grid  tube.  Table  1.2  lists  the  flow  configurations  testa  with  the  tur¬ 
bulence  jet-grid  placed  at  the  second  station  (designated  as  station  02)  upstream  from  the 
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turbine  blade  row.  Configurations  tested  without  any  jet-grid  placed  in  the  free-streaui  are 
denoted  with  the  grid  location  (£)  denoted  appropriately  as  N/A.  St  ream  wise  secondary 
air  injection  is  denoted  as  co-flow,  injection  perpendicular  to  the  frcc-strcam  is  denoted  as 
cr-flow,  and  injection  against  the  frcc-strcam  is  denoted  as  ct-flow  in  Tables  1.1  and  1.2. 

The  derivation  of  turbulence  scales  and  the  theory  of  laminar  and  turbulent  bound¬ 
ary  layer  heal  transfer  arc  discussed  in  Chapter  2.  The  experimental  apparatus  used 
to  perform  the  experiment  and  instrumentation  are  discussed  in  Chapter  3.  Chapter 
4  presents  detailed  experimental  procedures  used  to  measure  the  turbine  blade  surface 
pressures,  the  frcc-strcam  turbulence  and  turbulence  integral  and  microscale  lengths,  and 
measure  the  heat  transfer  on  the  turbine  blade  surface.  The  results  of  this  study  are  dis¬ 
cussed  in  Chapter  5.  The  results  of  the  blade  surface  pressure  measurements  are  presented 
as  nondimensional  pressure  coefficients  (C*,,).  Doth  turbulence  integral  scale  length  and 
microscale  length  are  calculated  for  the  primary  tunnel  reference  directions.  Surface  heat 
transfer  on  the  blade  is  presented  in  the  form  of  nondimensional  Nusselt  numbers  (;Y«/J 
and  Stanton  numbers  (S/). 
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Table  1.1.  Test  configurations  for  the  turbine  cascade  wind  tunnel-station  0 


station 

f 


no  grid 


15.75 


15.75 


15.75 


15.75 


injection 

direction 


co*  flow 


co- flow 


co-flow 


co-flow 


co-flow 


co-flow 


cr-flow 


cr-fiow 


cr-flow 


cr-flow 


cr-flow 


cr-flow 


cr-flow 


cr-flow 


cr-flow 


cr-flow 


cr-flow 


orijtcc 
diameter  D 
mm 


N/A 


N/A 


.66 


.66 


0.66 


0.66 


0.66 


1.32 


1.32 


1.32 


2 


1.98 


1.98 


injection 

pressure 

kPa 


N/A 


N/A 


83.5 


156.3 


2-15.8 


300. 


106.2 


i 


157.5 


2-13.1 


296.3 


78.0 


160.5 


233.9 


291.8 


83.-1 


15S.2 


2-17.1 


303.3 


80.0 


153. -I 


238.9 


293.3 


150.2 


228.7 


cn 

15.75 

0.66 

83.9 

C12 

15.75 

eses 

0.66 

163.5 

C13 

15.75 

ct-flow 

0.66 

238.6 

CM 

15.75 

0.66 

295.1 

ct-flow 

1.32 

81.1 

C22 

15.75 

ct-flow 

1.32 

155.0 

C23 

15.75 

1351291 

1.32 

238.8 

C24 

15.75 

takiEEM 

1.32 

292.-1 

C31 

15.75 

ct-flow 

1.98 

80.7 

C32 

15.75 

ct-flow 

1.98 

1-13.2 

C33 

15.75 

ct-flow 

1.9S 

241.1 

1  C3-1  |  15.75 

ct-flow 

1.98 

293.2 

Table  1.2.  Test  configurations  for  the  turbine  cascade  wind  tunnel-station  02 


//.  THEORY 


2.1  Turbulence  Integral  Scale  ami  Microscale  Lengths 


The  turbulence  integral  scale  (A)  an<?  microscale  (A)  arc  of  interest  from  the  stand* 
point  of  turbulence.  The  integral  scale  define*  the  *izc  of  the  large  turbulent  eddies  in 
a  turbulent  flow.  The  microscale  is  a  measure  of  the  smallest  measurable  homogeneous 
group  of  fluid,  where  all  particles  within  the  group  share  the  same  velocity.  Figure  2.1 
shows  the  concept  of  turbulence  integral  scales  and  microscales.  Note  that  the  turbulence 
integral  scales  arc  comprised  of  microscales. 


For  analysis  purposes,  consider  the  velocity  of  a  fluid  in  a  single  direction. 

The  velocity  in  the  X*direction  (u)  may  be  considered  as  consisting  of  a  steady  state 
component  («)  and  a  fluctuating  component  (u')  of  velocity. 


ii  =  u  +  u' 


(2.1) 


In  the  same  manner  that  vibrations  in  a  solid  occur  at  different  frequencies  and 
strength,  turbulence  in  a  fluid,  such  as  air,  occurs  at  different  'Vequencies  and  with  different 
strengths  of  turbulence.  These  frequencies  are  measurable  and  related  to  the  turbulence 
scale  lengths  by  the  spectral  distribution  of  the  tnrbulei.ee.  If  /?„*(«)  defines  the  spectral 
distribution  of  u,:  in  the  frequency  bandwidth  between  n  And  n  +  i/n,  then 


(2.2) 


In  Equation  2.2,  the  quantity  dn  is  the  unit  frequency  interval.  £u»(n),  the  turbulence 
l>owcr  per  frequency  bandwidth  (An  ),  is  obtained  from  the  Fourier  transform  of  u'2  from 
the  time  domain  to  the  frequency  domain  in  the  following  manner.  . 

For  a  discrete  function,  Parseval’s  theorem  is  written  (Brigham,  198G)  as 


A’-l 


i«  —  ■  ,  A'-l  . 

(«<»'>) 


(2.:i) 


k=0 


n=0 


2-1 


2-2 


where 

//(«') 


=  spectral  function  of «'  from 
the  Fourier  transform  of  u' 


or  rewriting  Equation  2.3, 

„■  fwi1 
A’  J 

Cowiaring  Equation  2.4  with  Equation  2.2,  £„*(n)  must  be  related  to  as 


(2,1) 


KAn)  = 


«(»•) 


;V 


1 

An 


(2.o  J 


Vmu:  tKt*  »h^ry  of  Fourier  Transforms,  the  longitudinal  correlation  fun*  tion  (/(*))  of 
the  special  distribution  of  turbulence  is  related  to  the  spectral  power  <lv!r.;ty  (Hinze,  1959). 
The  function  f[x)  is  related  to  £„*(»)  as 

/(*)  sir  £„•(«)  COS  <tn  (2.0) 

Uu  Jo  \  «  / 

By  inverse  transform  relation,  Equation  2.6  becomes 


(2.7) 


The  primary  derivation  ol  the  longitudinal  correlation  function  is  for  a  time  series. 
In  Equations  2.6  and  2.7,  time  ( l )  has  been  replaced  using  Taylor’s  hypothesis.  Namely, 


/  =  ~ 
u 


As  n  -»  0,  then  Equation  2.7  becomes 


lim  - )  =  r  f(x)dxs. 

n-°  \A ui2EA")J  Jq 


(2.S) 


where 
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a  X-coinponeni  turbulence  integral  scale 


Am 


The  microscale  (A«)  is 


Differentiating  Equation  2.6  twice  yields 


(2.9) 


(2.10) 


In  experimental  measurements  of  turbulence,  the  integrals  in  Equations  2.9  anil  2.10 
are  replaced  by  summations  over  a  finite  frequency  domain,  £„*(”)  *s  computed  using 
Equation  2.5,  ami  »  is  a  finite  frequency  bandwidth  defined  as 


1 

"  =  (AT,  .AD 


(2.11) 


where 

.VM  =  total  sample  size 

Jit  =  sample  time  interval 

The  foregoing  analysis  is  based  on  a  similar  analysis  presented  by  Ccbecci  and  Smith 
(Cebcci,  1974).  In  a  similar  manner,  the  turbulence  scales  are  determined  for  any  other 
coordinate  direction  of  interest. 

2.2  Turbulence  Intensity  Decay 

The  decay  of  turbulence  in  a  frec-stream  flow  is  a  measure  of  the  dissipation  of  the 
turbulence  microscale  lengths.  The  dissipation  is  in  the  form  of  turbulence  microscale 
growth,  where  the  microscales  dissipate  their  energy  by  combining  with  other  microscales. 
Studies  by  Young  and  Han  (Young,  1989)  and  Blair  (Blair,  1983a)  indicate  the  microscale 
dissipation  is  related  to  the  size  of  the  turbulence  generating  grid.  Blair  offers  the  following 
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correlation  for  dissipation  of  the  turbulence  level  [Tu  ), 


Tu  =  0.7$ 


(2.12) 


Equation  2.12  is  applicable  only  where  the  free-stream  turbulence  is  uniform.  Diair 
determined  uniform  free*strcam  turbulence  is  not  reached  until  approximately  ten  grid 
tube  diameters  downstream  of  the  turbulence  generating  grid;  therefore  Equation  2.12  is 
only  applicable  for  J  >  10. 

The  studies  of  Young  and  Han,  found  excellent  agreement  with  Equation  2.12  with 
the  jet-grid  place  in  the  flow,  but  with  no  injection.  With  secondary  flow  injection  through 
the  jet  grid,  the  turbulence  levels  were  generally  higher  than  those  predicted  by  Equation 
2.12.  The  dissipation  of  the  integral  scale  was  not  correlated.  In  a  turbine  cascade  wind 
tunnel,  the  turbulence  decay  has  not  been  correlated. 


2.H  Laminar  1 1 oil  Transfer  on  a  Surface  with  n  Pressure  Gradient 

Laminar  heat  transfer  (or  flow  over  a  consta.il-tomperature  body  of  arbitrary  shape 
is  developed. 

The  rate  of  growth  of  any  thermal  boundary  layer  thickness  A: /,  is  a  function  of 
local  parameters  only  (Eckert,  1042). 


ds 


ds 


where 

=  frec-strcam  velocity  on  the  boundary  layer  edge 
«  =  the  distance  along  the  surface  of  the  body 

v  =  local  kinematic  viscosity 

Pr  -  Prandtl  number  of  the  fluid 

Using  the  thermal  conduction  thickness  (A*  =  £),  where  k  is  the  fluid  thermal 
conductivity,  and  h  is  the  convection  heat  transfer  coefficient,  then  the  variation  of  A.t  is 


given  in  nendimensional  form  as 

f(*&  dI” 

v  dz  “  7  \v  ’  ds  ’ 

Assuming  the  function  /  is  the  same  as  for  the  family  o:  wedge  flows,  then  Equation 
2.13  is  readily  evaluated  using  similarity  solutions. 

From  Kays  and  Crawford  (Kays,  1987),  for  flow  over  a  flat  plate  and  Pr  =  0.7, 

A'urtc,-*  =  Ci 


from  which 

h  3i 


where  C‘i  is  a  coefficient  which  depends  on  in,  the  wedge  parameter,  and  Pr.  Employing 
the  wedge  condition. 

*  CV" 

the  following  expressions  are  developed, 


t'cc  _  1  -  »i 

v  As  =  (C,)3 

A43  dUcc  Ml 

1/  d$  ~  (C,)2 


(2.1*1] 


(2.15) 


In  the  expressions  above,  mi  is  the  wedge  flow  parameter  is  defined  as 

mi  S  -X-  (2.10, 

*  * 

where  0W  =  to  the  wedge  half-angle. 

For  Pr  =  0.7,  Equations  2.14  and  2.15  arc  plotted  for  the  exact  wedge  solutions. 
Smith  and  Spaulding  (Smith,  1958)  propose  a  linear  fit  to  the  wedge  flow  solution 
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Figure  ’2.2.  Wedge  How  solution  for  a  laminar  boundary  layer,  (Kay*.  19s«| 


shown  in  Figure  2.2  of  the  form 


- —  =  ll.CS-2.S7~~- 

n  as  a  its 


12.IT » 


Equation  2.17  is  integrated  and  reduced  to 

a  2  —  1  I.GsV/q  l  *'*'ds 
1  =  /•  3.S7 

t  * 


12. Is  I 


Noting  A*  s  j,  Equation  2. IS  is  j»ut  into  the  form  of  a  Stanton  number 
Pr  =  0.7  is  assumed. 


St  =  0.4  IS 


I- •(ptu0-'"' 


|/o(/>tW-*Td.« 


0.5 


(2.19) 


where 

ft  =  local  viscosity 

p  =  local  fluid  density 


Equation  2.19  provides  a  useful  method  for  calculating  the  laminar  heat  transfer  fur  a 


two-dimensional  body.  Application  of  the  laminar  heat  transfer  theory  to  the  turbine  blade 
Is  necessarily  restricted  to  those  areas  along  the  blade  that  exhibit  laminar  or  laminar-like 
flow  conditions. 


2,4  Heat  Transfer  on  an  Arisymetric  Hody  with  Arbitrary  Varying  Ftxcsttxtm  Velocity 
and  Surface  Temperature 

Prediction  of  surface  heat  transfer  for  a  turbine  blade  has  been  suggested  (Ambrok,  19.'?  \ 
to  include  the  arbitrarily  varying  pressure  gradient  on  an  axisymmetric  body  with  a  tur¬ 
bulent  boundary  layer.  The  model  docs  not  account  for  complex  flow  phenomena,  such 
as  a  secondary  vortex  core,  or  laminar  separation  bubbles  that  may  be  present  within  a 
turbine  cascade. 

for  a  constant  surface- temperature,  constant  free- stream  velocity,  the  Stanton  num¬ 
ber  may  be  expressed  in  the  form 

$.*  =  ah;*  v>m 


where 

C*  =  a  coefficient 

a  Reynolds  number  based  on  surface 
distance  $  raised  to  the  -«  power 


The  energy  integral  for  a  flat  plate  is  given  as 


SI  = 


ds 


12-21) 


where 

Aj  =  enthalpy  thickness 


Substituting  Equation  2.20  into  Equation  2.21,  and  integrating  yields 


A?  = 


CU^ns'-n 
i/"(l  -  n) 


(2.22) 


where 
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=  the  local  frees*  ream  velocity 
=  local  kinematic  viscosity  of  the  fluid 


Ms 


U 


Equation  2.22  is  substituted  into  Equation  2.20'to  jive  the  Stanton  number  as  a  function 
of  the  Reynolds  number  with  respect  to  Aj  rather  than  s. 

Si  a  (2.23) 

Equation  2.23  does  not  presume  any  previous  history  of  the  turbulent  boundary  layer  on 
the  surface. 

For  turbulent  boundary  layers  on  a  flat  plate,  from  (Kays,  10S7) 

C  a  0.02S7  />r“0H 
and 

it  «  0.2 

Pr  ss  Prandtl  number  of  the  fluid  How 
Evaluation  of  Equation  2.23  using  the  given  values  of  C  and  n  give 

St  -  0.0 l25/,r“wffcA9"°’a*  (2.2*1) 

In  order  to  evaluate  the  local  surface  heat  flux  in  terms  of  the  local  surface  temper* 
mure  difference,  the  integral  energy  equation  is  used. 

?  =  ^^[A,«p„u„(3r,  -  r„)|  (MW 

where 

q'  =  heat  flux  per  unit  surface  area 

Cp  =  specific  heat  at  constant  pressure 
Poo  =  local  free-stream  density 

2’co  =  local  froo-stream  temperature 

T,  =  local  surface  temperature 

R  =  axisymmetric  radius 
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Noting 


(2.20) 


6<~('n-7«)/>ccCP(^ 
then  Equation  2.25  may  be  algebraically  manipulated  into  the  form 


-  T^d,  =  d[pwy«A2  «('/',  -  rcojjris  (2.27) 


where 

/i  =  local  fluid  viscosity 

Equation  2.27  may  be  integrated  directly  if  we  assume  at  the  stagnation  point  Utc  -  0 
where  $  -  0.  Using  Equation  2.2G  with  Equation  2.27  gives 


St  =  Q.02S7Pr-(M 


R^HT,  -  7'co)0,25//0,2 
Hi  -  T^p^ds)™ 


(2.2S) 


where  C  and  «  are  the  same  numerical  values  previously  proposed. 

For  the  turbine  blade,  R  is  a  constant,  and  />«  is  relatively  constant  for  low  Mach 
numbers.  Equation  2.2$  is  further  simplified  for  the  above  conditions  as 


St  =  0.02$77,r-°--‘ 


(T,  -  7«)0‘*si,0-2 

I/o(7j  - 


(2.29) 


Equation  2.29  may  be  evaluated  for  turbulent  flow  over  an  axisymmctric  body  with  a 
pressure  gradient.  Others  (Kays,  19S7)  have  found  Equation  2.29  to  be  in  good  agreement 
for  a  variety  of  applications,  including  turbine  blades.  However,  it  fails  badly  for  strongly 
accelerated  flows  where  A*  >  1  x  10“G.  K  is  defined  as  the  acceleration  parameter 


K  = 


v  dUoo 

VI  ds 


(2.30) 


For  very  strongly  accelerate-)  flows,  K  >  3.5  x  10"6,  the  entire  turbulent  boundary  layer 
begins  to  rclaminamc  as  the  viscous  sublayer  thickens,  engulfing  the  entire  boundary  layer 
(Kays,  1970);  hence,  this  theory  no  longer  applies  for  this  condition. 
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III.  EXPERIMENTAL  APPARATUS 


EXPERIMENTAL  APPARATUS 


3.1  Turbine  Cascade  Wind  Tunnel 


Testing  was  conducted  in  a  large  scale,  low-aspect  ratio,  linear  turbine  cascade  wind 
tunnel  shown  in  Figure  3.1.  The  facility  is  owned  by  the  Air  Force  Institute  of  Technology 
(AFIT)  Aeronautical  Engineering  Department,  Wright- Patterson  AFB,  Ohio.  Specific 
tunnel  dimensions  and  tunnel  flow  characteristics  arc  listed  in  Table  3.1.  The  wind  tunnel 
is  a  draw-down  tunnel,  powered  by  a  Buffalo  Forge,  Model  No.  BL  *  305,  centrifugal  fan. 
The  fan  was  powered  at  a  constant  velocity  by  a  14.9  kW,  2-10  volt  electric  motor.  The 
tunnel  wind  velocity  was  manually  controlled  with  a  series  of  variable  angle  inlet  vanes 
which  were  set  at  an  angle  opposite  to  the  direction  of  revolution  of  the  fan;  hence,  tended 
to  reduce  flow  swirl  normally  associated  with  centrifugal  fans. 


3.2  'Turbine  Blade  Cascade  and  1'urbinc  Blade  Pivfdc 

A  photograph  of  the  turbine  blade  cascade  lest  section  is  shown  in  Figure  3.2.  The 
turbine  blade  cascade  consists  of  four  geometrically  identical  turbine  blade  sections.  A 
reference  plane  was  arbitrarily  established  at  the  leading  edge  of  the  turbine  blades  and 
a  reference  axis  system  established  with  the  origin  in  the  reference  plane.  Figure  3.3  is 
a  schematic  of  the  turbine  blade  test  section,  the  reference  axes,  and  the  reference  plane. 
The  Z-axis  is  into  the  paper,  with  Z  =  0  corresponding  to  the  point  located  on  the  inside 
of  the  upper  cndwall  of  the  tunnel.  The  X-axis  was  defined  in  the  direction  parallel  to  the 
incoming  flow  incidence  angle.  The  Y-axis  lies  in  the  reference  plane,  at  a  134.7  degree 
angle  measured  in  a  clockwise  direction  from  the  norma]  to  the  reference  plane.  The  wind 
tunnel  flow  direction  was  adjured  with  inlet  side  boards  so  the  flow  incidence  angle  to  the 
turbine  blade  row  was  44.7  degrees,  as  measured  in  a  clockwise  direction  from  the  normal 
of  the  reference  plane.  This  incidence  angle  was  chosen  to  match  earlier  tests  using  an 
identical  turbine  blade  profile  (Moore,  1984c),  (Langston,  1977). 
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Figure  3.1.  The  A  FIT  low-aspect  ratio,  linear  turbine  cascade  wind  tunnel 


Tabic  3.1.  Turbine  cascade  wind  tunnel  dimensions  and  How  parameters 


ll'iMrf  Tunnel  Dimensions  and  Flow  Parameters 

nozzle  dimension 

0.53G  m  x  0.536  m 

nozzle  contraction  ratio 

11.5:1 

nozzle  length 

1.025  in 

inlet  height 

0.114  m 

inlet  width 

0.218  in 

inlet  length 

0.450  m 

test  section  height 

0.114  m 

tunnel  exit  length 

0.650  m 

tunnel  inlet  to  exit  angle 

75.0° 

inlet  flow  velocity  range 

—T")l  in 

exit  flow  velocity  range 

mwiT'nni 

diffuser  length 

1.922  m 

freestream  turbulence  level 

1.26  %  ** 

♦*  Measured  at  £  =  15.75,  $  =  8.00,  and  £  a  0.00 


Two  locations  upstream  of  the  reference  plane  were  chosen  for  placement  of  the  jet- 
grid  device.  The  first  position  (station  0)  was  at  a  nondimensional  distance  of  |  =  15.75 
ahead  of  the  reference  plane,  where  x  is  the  linear  distance  and  b  is  the  jet-grid  tube 
diameter.  The  second  position  (station  02)  of  the  jet-grid  was  at  £  =  25.07  upstream  of 
the  reference  plane.  Figure  3.3  shows  the  position  of  each  station  relative  to  the  reference 
plane. 

The  turbine  blade  true  chord  was  114.3  mm,  with  a  unity  chord  to  span  aspect  ratio. 
Based  on  the  blade  axial  chord  length  (c  )  the  range  in  the  exit  Reynolds  number  [Ilce  ) 
was  3.0-7.0  x  105.  The  turbine  blade  nomenclature  is  given  in  Figure  3.4.  The  flow  turn 
angle  through  the  blade  row  was  105  degrees. 

3.3  Turbulence  Generating  Jet-Grid  Device 

Free-stream  turbulence  was  generated  with  a  jet-grid  device.  The  jet-grid  device 
allows  the  user  to  vary  the  intensity  of  freestream  turbulence  by  varying  the  injection 
pressure  (secondary  mass  flow),  and  by  varying  the  direction  of  the  injection  jets.  Ad¬ 
ditionally,  the  turbulence  intensity  may  be  varyied  through  placement  of  the  jet-grid  at 
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Figure  3.2.  Turbine  blade  cascade  lest  section 
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Table  3.2.  Jet-grid  injection  device  specifications 


Jet-Grid  Specification* 

number  of  tubes 

10 

tube  spacing  (center  to  center) 

25.40  mm 

tube  outer  diameter  (6) 

6.350  mm 

number  of  injection  orifices/tubc 

7 

total  number  of  orifices 

70 

smallest  orifice  diameter 

0.66  mm 

medium  orifice  diameter 

1.32  mm 

largest  orifice  diameter 

1.98  mm 

(lilfcrent  stations  relative  to  the  test  section.  Pscvious  studies  on  flat  plate  surface  heat 
transfer  have  employed  such  a  device  with  remarkable  success  (Young,  10S9).  The  jet-grid 
device  employed  by  Young  and  Han  (Young,  1980}  was  limited  to  secondary  flow  injection 
in  a  downstream  or  upstream  flow  direction. 

This  study  used  a  jet-grid  device  consisting  of  ten  vertical  aluminum  tubes  with  a 
scries  of  evenly  spaced  orifices  drilled  along  the  longitudal  axis  of  each  tube.  The  tubes 
were  placed  directly  into  the  free-stream  flow  through  the  wind  tunnel  endwall.  The  wind 
tunnel  endwall  was  drilled  to  accept  the  tubes.  Each  tube  was  secured  in  place  using  a 
simple  holding  bar  attached  to  the  outer  surface  of  the  endwall.  Each  tube  was  blocked  on 
one  end,  with  the  other  end  attached  to  flexible  tubing.  Figure  3.5  shows  a  photograph  of 
the  jet-grid  device  and  Table  3.2  lists  specifications  of  the  jet-grid  device. 

Three  primary  injection  directions  were  used  during  this  study;  streamwise  with  the 
free-stream  flow,  p  crpendiducar  to  the  free-stream,  and  opposing  the  free-stream,  and  were 
designated  throughout  this  study  as  co-flow,  cross-flow,  and  counter-flow,  respectively.  A 
schematic  of  the  jet-grid  injection  directions  is  given  in  Figure  3.6. 

All  feeder  lines  to  the  jet-grid  were  attached  to  a  pressure  manifold  which  allowed 
even  distribution  of  the  secondary  flow  to  each  tube  prior  to  injection.  The  injection 
pressure  was  measured  with  an  Endevco,  Model  1675-2A,  685,000  kPa  differential  pressure 
transducer,  centrally  located  along  the  length  of  the  manifold.  A  single  (2.54  cm  inner 
diameter)  high  pressure  air  line  was  attached  to  each  end  of  the  manifold.  The  incoming 
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Figure  3.5.  Jet-Grid  injection  device 
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Figure  3.(5.  Schematic  of  orifice  injection  directions 

air  was  triple-filtered  and  dried  using  a  simple  mechanical  air-moisture  separator.  Filtering 
of  the  air  was  necessary  for  protection  of  the  hot-wire  probe  and  to  prevent  contamination 
of  the  hot-wires  and  the  surface  temperature  test  surface.  Two  separated  sources  of  high- 
pressure  air  were  needed  to  provide  adequate  mass  flow  through  the  larger  diameter  orifices 
and  maintain  a  steady  injection  manifold  pressure.  The  primary  air  source  was  low-volume, 
high-pressure  shop  air;  the  secondary  air  source  was  high-volume,  pressurized  air  from  an 
Ingersoll-Rand,  Model  55*1  E9,  three  stage  air  compressor,  with  a  1,200  liter  primary  storage 
tank  and  a  *1,000  liter  secondary  dump  tank. 

3.4  Hot-Wire  Anemometers 

Instantaneous  local  velocity  was  measured  using  a  Thermo  Systems  Model  No.  100 
(IFA-100)  Intelligent  Flow  Analyzer,  constant-temperature  anemometer  controller  and  a 
Thermo  Systems,  Model  No.  1240-T1.5,  serial  no.  7713,  X-type  tungsten  hot-wire  probe, 
also  called  an  X-wirc  probe.  This  X-wire  was  a  cross-flow  type  probe  which  allowed  entry 
into  the  wind  tunnel  from  the  endwall  with  minimum  interference  with  the  free-stream 
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flow.  The  probe  has  two  0.0038  mm  diameter  tungsten  wires,  offset  from  each  other  by 
approximately  90  degrees.  The  tungsten  hot* wire  probe  was  chosen  over  more  durable 
hot-film  probes  because  of  its  excellent  frequency  response  to  turbulence  in  the  flow.  The 
exact  intersection  angle  of  the  hol*wircs  was  measured  with  a  microscope.  Each  hot-wire 
output  signal  was  amplified  and  filtered  using  the  Thermo  Systems,  Mode)  No.  15S,  signal 
conditioner  prior  to  input  into  the  high-speed  data  acquisition  system. 

3.5  Surface  Pressure  7'cst  Surface 

The  surface  pressure  test  surface  was  an  aluminum  turbine  blade  instrumented  with 
23  static  pressure  ports.  The  leading  edge  static  pressure  port  was  located  at  37.8  percent 
of  the  span.  The  trailing  edge  static  pressure  port  was  at  the  GO  percent  span  position. 
Remaining  pressure  ports  were  on  the  imaginary  line  extending  from  the  leading  edge 
pressure  port  to  the  trailing  edge  pressure  port  on  both,  the  pressure  ami  suction  surfaces. 
This  staggered  distribution  of  static  pressure  ports  was  selected  to  alleviate  some  of  the 
effects  of  leading  orifices  disturbing  the  flow  on  the  turbine  blade  surface  prior  to  the 
static  pressures  being  measured  by  subsequent  pressure  taps.  The  pressure  port  taps  were 
connected  to  steel  tubing  which  extends  through  the  plexiglass  cndwall  of  the  wind  tunnel 
test  section.  Each  pressure  tap  was  connected  to  the  Scanivalvc  dS  port  selector  with 
flexible  plastic  tubing.  A  G.895  kPa  differential  pressure  transducer  measured  the  static 
pressures  at  each  pressure  port  on  the  blade.  Figure  3.7  is  a  photogragh  of  the  surface 
pressure  instrumented  turbine  blade.  The  staggered  distribution  of  the  pressure  taps  is 
visible  along  the  pressure  surface  in  this  figure.  The  location  of  each  pressure  tap  in  terms 
of  a  nondimcnsional  chord  length  (  f  )  and  a  nondimonsional  surface  distance  (  \  )  is 
given  in  Table  3.3.  The  distance  along  the  surface  of  the  blade  from  the  leading  edge  of 
the  turbine  blade  is  s ;  c  is  the  turbine  blade  true  chord  length. 

3.6  Heat  Tmnsfcr  Test  Surface  and  Thermocouples 

The  heat  transfer  test  model  was  constructed  of  an  expanded  urethane  foam  core 
with  a  stainless  steel  foil  surface  wrapped  around  the  exterior  of  the  foam  core.  The  test 
turbine  blade  surface  was  heated  by  passing  an  electric  current  through  the  stainless  steel 


3-10 


Table  3.3.  Location  of  static  pressure  taps  on  the  turbine  blade 


pressure  tap  no. 

£ 

c 

9 

€ 

pressure  sur/uce  (}>■<) 
tucthn  (urfacc  (tf) 

r~  -  j 

0.01739 

ps 

2 

0.0C522 

0.08772 

PS 

3 

0.1244 

0.1991 

P* 

*1 

0.1913 

0.3245 

ps 

5 

Msm 

0.4174 

ps 

G 

S 

0.6184 

ps 

7 

0.0000 

0.8026 

ps 

0.7GS7 

0.8974 

ps 

9 

0.S522 

0.5763 

P* 

10 

0.9235 

1.0290 

ps 

11 

1.0000 

1.1815 

ps 

12 

0.9026 

1.1167 

ss 

13 

0.7613 

1.0088 

ss 

14 

0.6670 

0.9035 

ss 

15 

0.5817 

0.8026 

ss 

10 

0.*1S70 

0.6974 

ss 

17 

0.3S52 

0.6404 

ss 

IS 

0.3322 

ss 

19 

MMO! 

0.3991 

ss 

20 

■TOTiW! 

0.2737 

ss 

21 

0.03565 

0.1298 

SS 

22 

0.0000 

0.07460 

ss 

23 

0.008C96 

ss 

True  chord  c 
Pressure  tap  no.  1 
Pressure  tap  no.s  2-10 
Pressure  tap  no.  11 
Pressure  tap  no.s  12  -23 


=  1 1*1.3  min 
=  leading  edge  port 
=  pressure  surface  ports 
=  trailing  edge  port 
=  suction  surface  ports 
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Figure  3.7.  Turbine  blade  instumcnted  for  surface  pressure  measurements 

foil.  A  50  amp,  2.5  volt  DC  current  was  needed  to  drive  the  temperature  of  the  turbine 
blade  test  surface  high  enough  to  provide  a  measurable  temperature  difference  between 
free-stream  flow  temperature  and  the  blade  surface  temperature.  The  urethane  core  was 
cast  from  a  mob!  of  the  pressure  instrumented  turbine  blade  in  order  to  assure  identical 
blade  geometry  with  the  pressure  test  surface.  A  small  amount  of  the  outer  surface  of  the 
urethane  foam  core  was  removed  to  accept  the  stainless  steel  foil  and  maintain  the  proper 
blade  dimensions.  The  dimensions  of  the  stainless  steel  foil  were  2.54  x  10”3  mm  thick. 
2G3.2  mm  in  length,  and  110.0  mm  in  width.  The  width  of  the  stainless  steel  foil  was  3 
mm  less  than  the  span  of  the  foam  core  in  order  to  prevent  an  electrical  short  between  the 
electrically  heated  test  surface  and  the  tunnel  endwalls.  Copper  bus  bars  were  attached  to 
each  end  of  the  stainless  steel  foil,  then  the  foil-bus  bar  assembly  was  glued  to  urethane 
foam  core  turbine  blade  using  unexpanded  urethane.  The  foil  surface  was  continuous  from 
one  bus  bar  on  the  pressure  surface  trailing  edge  to  the  other  bus  bar  on  the  suction 
surface  trailing  edge.  The  stainless  steel  foil  was  not  connected  at  the  trailing  edge.  A 
5  mm  section  of  urethane  foam  separated  the  bus  bars  and  prevented  electrical  contact 
of  the  bus  bars.  Figure  3.8  shows  the  turbine  blade  constructed  for  surface  heat  transfer 
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Figure  3.S.  Heat  transfer  test  turbine  blade 


tests.  The  copper  bus  bars  extended  beyond  the  span  of  the  turbine  blade  in  order  to  pass 
through  the  plexiglass  wind  tunnel  endwall.  The  test  turbine  blade  was  secured  within  the 
test  section  with  two  bolls  passing  through  each  endwall  into  the  side  of  the  turbine  blade. 
Metal  inserts  were  placed  within  the  urethane  foam  core  to  support  the  attachment  holts. 

Temperature  was  measured  using  23  iron-conslantan  (Type  J),  3G-gunge  thermocou- 
pics  attached  to  the  underside  of  the  foil  surface.  Along  the  interior  of  the  foam  blade,  five 
Type  J  thermocouples  were  placed  along  the  camber  line  and  directly  in  line  with  the  sur¬ 
face  thermocouples.  The  interior  thermocouples  were  used  to  account  for  heat  conduction 
into  the  turbine  blade.  Each  thermocouple  was  calibrated  for  temperature  sensitivity  and 
accuracy  using  an  ice  bath.  The  accuracy  of  the  thermocouples  was  ±  1  degree  Celsius. 
Each  thermocouple  was  attached  to  the  stainless  steel  foil  using  a  low-current  arc  welder. 
A  small  amount  of  silver-tin,  Sn  63,  solder  was  applied  at  each  thermocouple  contact  with 
the  steel  foil  for  added  support  of  the  thermocouple.  The  thermocouples  were  located  on 
the  heat  transfer  test  turbine  blade  at  the  same  position  as  the  pressure  port  taps  were 
located  on  the  surface  pressure  instrumented  turbine  blade.  Thus,  Table  3.3  gives  the  loca¬ 
tion  of  the  thermocouples  for  the  heat  transfer  test  turbine  blade  as  well  as  the  location  of 
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the  pressure  port  tnps  on  the  surface  pressure  instrumented  turbine  blade.  Thermocouples 
in  the  interior  of  the  turbine  blade  were  numbered  as  2*1.2$,  from  leading  edge  to  trailing 
edge,  respectively.  During  the  final  construction  phase  of  the  turbine  blade,  three  ther¬ 
mocouples,  numbers  1, 2  and  24  were  broken  and  remained  inactive  throughout  the  tests. 
All  other  thermocouples  remained  active.  All  thermocouple  leads  exit  the  turbine  blade 
around  the  quarter  chord  point  and  pass  through  the  tunnel  endwall.  Each  thermocouple 
lead  was  connected  to  the  high-speed  data  acquisition  system. 

3.7  fostumzntntion 

A  block  diagram  of  the  instrumentation  control  is  shown  in  Figure  3.9. 

A  Hewlett-Packard  High-Speed  Data  Acquisition  System,  Model  No.  HP3S52A,  con¬ 
trolled  the  majority  of  instrumentation  used  during  the  experiment.  The  HP3S52A  is  a 
mainframe  with  a  10-bit  processor  and  2  megabyte  internal  memory  buffer.  The  backplane 
is  configured  for  separate  component  inputs.  The  following  components  were  used  during 
this  study: 

(1)  HP44701A,  Integrating  Voltmeter. 

(2)  HP44702H,  High-Speed  Voltmeter. 

(3)  HP44713,  High-Speed  FET  Multiplexer. 

(4)  HP4472S,  Relay  Actuator. 


The  Integrating  Voltmeter  sampled  pressure  and  temperature  signals  input  to  the 
High-Speed  FET  Multipcxer.  A  second  High-Speed  Fet  Multipexer  was  used  in  conjunction 
with  the  High-Speed  Voltmeter  connected  with  a  dedicated  ribbon  cable  for  high-speed 
sampling  rates  (100,000  signals  per  second  per  channel).  The  high-speed  sampling  rates 
were  necessary  in  order  to  digitize  the  hot-wire  probe  output  signals  accurately.  The  Relay 
Actuator  controlled  the  Scanivalvc  drive  motor  and  the  power  supply  used  to  heat  the 
stainless  steel  foil  surface  of  the  turbine  blade. 

A  Zenith  24$  personal  computer  with  a  National  GPII3,  IEEE-4S8,  internal  interface 
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board  was  used  to  control  the  mainframe  and  transfer  data  to  floppy  disk  for  permanent 
storage.  All  data  acquisition  software  was  written  using  Microsoft  QuickBASIC. 

Measurement  instrumentation  used  my  be  classified  as  pressure  mcasu.ment  devices, 
turbulence  measurement  devives,  and  temperature  measurement  devices,  and  arc  discussed 
as  follows. 

The  wmd  tunnel  frcc-strcam  total  and  static  pressures  were  measured  with  aK.fc  E. 
Company,  0.63  cm  U-lubc  manometer,  using  distilled  water  as  a  pressure  indicator  medium. 
A  Scanivalve  II,  48  port  driver  motor  with  a  6.895  kPa  differential  pressure  transducer  was 
used  to  measure  surface  pressures  on  the  turbine  blade.  The  jet-grid  injection  manifold 
pressure  was  measured  with  an  Eudcvco,  Model  1675-2A,  6S95  kPa  pressure  transducer. 
Each  pressure  transducer  output  signal  was  calibrated  for  pressure  using  a  known  applied 
pressure.  Calibration  equations  were  entered  into  the  high-speed  data  aquisition  system 
where  pressure  transducer  output  signals  (volts)  were  converted  to  pressure  units  (kPa) 
and  transferred  to  permanent  storage  on  a  floppy  disk. 

Turbulence  measurements  were  made  with  the  Thermo  Systems  Intelligent  Flow  Ann- 
ly:ct\  IFA-100.  Hot-wire  probe  output  signals  were  amplified  and  filtered  using  the  Thermo 
Systems,  Model  No.  15S,  Signal  Conditioner.  Each  hot-wire  probe  was  calibrated  for  flow 
speed  and  angle.  Temperature  calibration  of  the  hot-wires  was  not  performed.  Instead,  the 
hot-wires  were  recalibrated  for  angle  and  velocity  calibration  when  the  ambient  tempera¬ 
ture  varied  more  than  5  degrees  Celcius  from  the  temperature  of  the  current  calibration. 
Appendix  C  lists  the  hot-wire  calibration  procedure.  All  hot-wire  signals  were  read  by 
the  High-Speed  Voltmeter  and  stored  into  the  data  acquisition  system  memory  buffer,  and 
later  transferred  to  floppy  disk  for  permanent  storage. 

Temperature  was  measured  using  iron-constantan  (Type  J)  thermocouples.  The 
high-speed  data  acquisition  system  converts  thermocouple  voltages  to  temperatures  using 
an  internal  ice-point.  A  single  Type  J  thermocouple  and  an  Omega  Digicator  was  used  to 
visually  indicate  the  frcc-strcam  temperature. 

A  Hewlett-Packard,  Model  No.  HPG727,  100  ampere,  20  volt  DC  Power  Supply, 
supplied  a  continuous  DC  current  through  the  stainless  steel  foil  surface  of  the  heat  transfer 
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lest  turbine  blade.  The  current  was  measured  with  the  power  supply  ampere  gunge.  A  3* 
place,  Hewlett-Packard  Digital  Multimeter  was  used  to  measure  the  heated  blade  resistance 
and  the  power  supply  voltage. 
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1  -  Blade  surface  pressure  ports 

2  «  Frec-stream  pitot  tube 

3  -  Hot-wire  probe  (2  channels) 

4  -  Electrically  heated  blade  surface 

5  -  Frec-stream  temperature  sensor 

6  -  Turbine  blade  thermocouples 


Figure  3.9.  Block  diagram  of  instrumentation  and  control  structure 
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IV.  EXPERIMENTAL  PROCEDURES  AND  DATA  REDUCTION 


LI  Turbine  Cascade  Test  Section  Parameters 

The  flow  within  the  wind  tunnel  was  adjusted  to  match  equivalent  test  conditions 
studied  in  earlier  tests  performed  on  the  same  turbine  blade  profile  (Moore,  1‘JSlc).  The 
primary  areas  of  concern  arc  the  cascade  inlet,  cascade  exit,  and  the  wind  tunnel  flow 
speed,  c.g.,  the  exit  Reynolds  number. 

4.1. /  Frcc-Strcatu  Flout  Inlet  Incidence  Angle  and  Velocity  Profile  The  free- stream 
flow  incidence  angle  was  measured  along  the  span  of  the  wind  tunnel  using  the  following 
procedure. 

A  dual  passage  pressure  probe  was  inserted  into  the  inlet  of  the  wind  tunnel.  The 
passages  of  the  pressure  probe  were  connected  one  to  each  side  of  a  U-lnbe  water  manome¬ 
ter.  Water  or  any  fluid  of  similar  density  gives  excellent  sensitivity  to  pressure  differences 
exerted  on  either  side  of  the  probe.  The  pitot  probe  was  rotated  along  the  tunnel  Z-axis 
until  the  differential  height  of  the  manometer  was  nulled.  The  corresponding  flow  incidence 
angle  was  read  from  an  attached  protractor  to  within  ±  1  degree.  The  flow  incidence  angle 
was  set  to  the  desired  value  using  adjustable  inlet  side  boards.  During  the  present  study, 
the  flow  incidence  angle  was  adjusted  to  *1*1.7  degrees  measured  from  the  reference  plane 
to  the  incoming  free-stream  velocity  vector.  See  Figure  3.3  for  the  layout  of  the  test  sec¬ 
tion  geometry.  Results  of  three  traverses  at  three  Z-axis  positions  are  presented  in  Figure 
*1.1.  All  traverses  were  at  the  same  position  upstream  of  the  reference  plane  (j-  =  17.5, 
b  =  jet-grid  tube  diameter  (6.35  mm)).  The  results  of  Figure  *1.1  represent  the  optimum 
flow  conditions  after  the  inlet  side  boards  were  adjusted.  From  Figure  *1.1,  it  can  be  seen 
the  flow  incidence  angle  was  not  completely  uniform;  however,  the  largest  deviation  from 
desired  incidence  angle  did  not  occur  directly  ahead  of  the  turbine  blade  cascade  section 
where  surface  pressures  and  surface  heat  transfer  were  measured. 

The  speed  of  the  incoming  flow  was  measured  for  uniformity  along  the  span  of  the 
cascade.  A  1.6  mm  diameter  pitot  probe  was  inserted  into  the  flow,  and  the  flow  speed 
was  measured  using  a  U-tube  manometer.  The  speed  of  the  flow  was  computed  using  the 
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Figure  4.1.  Turbine  cascade  (low  incidence  angle.  Measurement  station:  £  =  15.75 


IlernouUi  incompressible  (low  equation. 


where 


P0  —  /«  + 


=  free-stream  total  pressure 
=  free-stream  static  pressure 
=  free-stream  air  density 
=  free-stream  flow  velocity  (speed) 


The  quantity  \pocUca2  is  the  dynamic  pressure  and  is  denoted  as  q w.  Manipulation  of 
Equation  4.1  gives 

U03=\—(P0-P00)Y  (4.2) 

tPoo  J 

Figure  4.2  presents  the  results  of  a  typical  traverse  along  the  span  of  the  wind  tunnel 
inlet.  The  free-stream  flow  speed  exhibited  the  same  trend  as  the  flow  incidence  angle 
results  presented  above. 
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Figure  4.2.  Frcc-strcam  flow  velocity  along  the  inlet  span 


Table  4.1.  Wind  tunnel  inlet  boundary  layer  measurements 


tunnel  surface 

f 

X 

6 

& 

t 

boundary  layer 
thickness  (mm) 

side  board 

17.5 

20.1 

8.0 

5.0 

side  board 

17.5 

22.4 

8.0 

4.0 

cud  wall 

17.5 

0.0 

0.0 

3.5 

end  wall 

17.5 

0.0 

1S.0 

4.1 

The  boundary  layer  along  the  inlet  side  boards  and  wind  tunnel  end  walls  was  mea* 
sured  with  the  pitot  probe.  Boundary  layer  measurements  for  the  turbine  cascade  wind 
tunnel  inlet  are  presented  in  Table  4.1. 

f.1.2  Frcc-Strcam  Flow  Exit  Angle  and  Velocity  Profile  The  flow  exiting  the  turbine 
cascade  was  measured  for  periodicity  of  the  flow.  The  flow  exit  speed  was  approximately 
85  m/s  (average).  At  this  speed,  the  flow  angle  probe  used  in  the  inlet  incidence  angle 
measurements  was  overly  sensitive  to  any  fluctuation  in  the  angle  of  incoming  flow;  hence, 
flow  exit  angle  was  not  measured.  Rather  exit  speed  was  measured  with  a  pitot  probe  for 
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Figure  4.3.  Exit  velocity  profile  of  turbine  cascade  wind  tunnel 


periodicity.  The  exit  speed  was  found  to  be  periodic,  with  slowed  How  regions  immediately 
behind  the  trailing  edge  of  the  turbine  blade.  From  Figure  4.3,  the  How  speed  varied  from 
C3.0  m/s  along  one  side  of  the  exit  section  to  91.3  m/s  along  the  opposite  side  of  the  exit. 
Results  are  given  for  a  single  traverse  at  a  selected  station  ( j  =  -32.0,  §  a  S.0).  The 
wake  of  turbine  blades  2  and  3  arc  readily  apparent  as  regions  of  slowed  flow.  Although 
the  How  exit  velocity  is  non-uniform  along  the  span  of  the  tunnel  exit,  in  the  region  of 
turbine  blade  3,  where  blade  surface  pressures  and  surface  heat  transfer  was  measured,  the 
exit  velocity  was  nearly  uniform  (except  in  the  wake). 


Neither  the  flow  exit  angle  nor  the  flow  exit  speed  were  adjusted  in  this  study. 
Although  specific  flow  exit  speed  profiles  were  not  presented  by  Moore  and  Ransmayr, 
they  conclude  the  exit  flow  exhibits  a  “strongly  nonuniform  spanwise  distribution  of  losses” 
(Moore,  1984c).  The  word  ”spanwise”  as  it  was  used  by  Moore  applies  to  the  span  of  the 
turbine  blade,  not  to  the  span  of  the  cascade  lest  section.  It  is  likely  flow  exit  conditions 
exhibited  during  this  study  were  similar  to  exit  flow  conditions  of  other  linear  turbine 
cascade  studies. 


f.1.3  inlet  and  Ezit  Reynolds  Number  Variation  The  free-stream  velocity  increas'd 
dramatically  as  the  flow  passed  through  the  turbine  cascade.  Some  portions  of  the  exiting 
flow  were  over  100  percent  faster  than  the  incoming  flow.  The  exit  Reynolds  number  (A’cc) 
was  chosen  as  the  free-stream  velocity  figure  of  merit  for  this  study.  The  exit  Reynolds 
number  is  defined  as 


(4.3) 


where 

o  =  exit  free-stream  velocity 

c  =  axial  chord  length  of  the  turbine  blade 

i/©8  ss  exit  free-stream  kinematic  viscosity 


As  discussed  in  the  previous  section,  the  location  of  the  measurement  of  the  exit  flow 
velocity  was  critical  since  the  exiting  flow  velocity  profile  was  highly  nonuniform.  The  flow 
exit  velocity  was  measured  in  the  geometric  center  of  the  tunnel  exit.  The  station  location 
was  £  =  -32.0,  *  =  2-1.0,  and  §  =  $.0. 

The  exit  Reynolds  number  for  this  wind  tunnel  varied  from  3.0-5.2  x  105.  All  tur¬ 
bulence  and  heat  transfer  measurements  in  this  study  were  made  at  Rce «  5  x  IQ5. 


Turbine  Elude  Surface  Pressure 

The  pressure  coefficient  (Cj, )  on  the  turbine  blade  is  computed  as 


P-Poo 

9» 


(4,1) 


where 

P  =  static  pressure  on  the  turbine  blade 

Poo  =  free-stream  static  pressure 

9co  =  free-stream  dynamic  pressure 

The  turbine  blade  static  pressures  were  measured  with  a  differential  pressure  trans¬ 
ducer.  Static  pressures  along  the  surface  of  the  instrumented  blade,  the  free-stream  total 
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pressure  (P0  ),  ami  free- stream  static  pressure  were  scanned  using  the  *18  port  Scanivalve 
port  selector.  Pressures  were  read  into  the  Hewlett-Packard  data  acquisition  system  as 
voltage  signals.  Each  static  pressure  signal  was  sampled  400  times,  averaged  to  a  single 
value,  then  the  average  output  signal  was  converted  to  a  pressure  using  an  internally  stored 
calibration  equation  which  was  previously  determined  for  this  specific  pressure  transducer. 
The  coefficient  of  pressure  was  computed  using  Equation  4.4.  Subroutine  PRESS.  DAS  of 
MAINDRV.BAS  listed  in  Appendix  E  automated  the  aforementioned  pressure  measure¬ 
ment  process. 

4.8  Jet-Grid  Injection 

The  jet-grid  was  introduced  into  the  free-stream  flow  at  two  upstream  stations.  The 
stations  were  denoted  as  station  0  (  y  =  15.75);  and  as  station  02  (£  =  25.07).  The  upper 
end  wall  of  the  tunnel  was  drilled  to  accept  the  jet-grid  tubes.  Flow  was  injected  as  either 
co-llow,  cross-flow,  or  counter-flow,  at  a  pressure  range  from  70  kPa  -  310  kPa,  depending 
on  the  test  being  performed.  Tables  1.1  and  1.2  lists  configurations  tested. 

Secondary  flow  was  slowly  introduced  to  the  jet-grid  pressure  distribution  manifold 
until  the  desired  pressure  was  obtained.  Filtered  shop  air  (6.55  x  105  kPa)  provided  the 
primary  feed  flow  to  the  pressure  distribution  manifold.  A  pressure  regulator  maintained 
constant  flow  pressure.  Distribution  manifold  pressures  were  measured  with  an  embedded 
0.985  x  10°'  kPa  pressure  transducer  ami  pressures  were  recorded  by  the  Hewlett-Packard 
data  acquisition  system. 

The  larger  diameter  orifice  tubes  required  a  second  source  of  feed  flow.  An  Ingcrsoll- 
Rand  three  stage  compressor  with  5,200  liter  storage  capacity  provided  a  second  source 
of  high-volume  flow.  During  the  highest  demand  for  secondary  flow,  a  steady  injection 
pressure  could  be  maintained  for  approximately  2  minutes,  using  both  feed  sources  of 
secondary  air  flow. 

The  injected  air  flow  temperature  differed  from  the  free-stream  temperature  by  as 
much  as  «  4  degrees  Celsius.  For  computations  involving  the  free-stream  temperature, 
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the  flow  temperature  was  measure*!  downstream  of  the  jet-grid.  Heat  transfer  calculations 
were  based  on  the  "mixed”  flow  temperature. 


4-4  Time  Discrete  Velocity  and  Turbulence  Intensity 

The  time  discrete  velocity  of  the  air  flow  was  measured  using  an  X-wire  probe.  The 
X-wire  selected  was  a  Thermo  Systems  Model  No.  1240-TI.5.  This  hot-wire  anemometer 
used  a  pair  of  0.0050S  nun  diameter  tungsten  wires  which  offered  excellent  response  to 
fluctuations  in  the  flow  and  allow  measurement  of  turbulence  microscale  lengths  on  the 
same  order  of  magnitude  as  the  wire  diameter.  A  brief  discussion  of  the  velocity  equations 
and  the  turbulence  equations  used  in  this  study  follows.  Appendix  C  gives  further  details 
on  the  calibration  of  the  hot-wire  probe  for  vclocity/turbulcnce  measurement. 

The  geometry  of  the  hot-wires  in  normal  configuration  is  shown  in  figure  4.4.  This 
is  the  same  geometry  assumed  for  calibration  of  the  X-wire  probe.  From  Figure  *l.-l 


=  f-*»-/3  (d.C) 

If  we  let  o  =  !-*.,  then  8\  and  82  may  defined  as 

8\  —  o  4*  ft  (*l.f ) 

82  —  Q  —  ft  (4.S) 

The  effective  velocity  each  hot-wire  measures  is  defined  for  each  wire  as 

f/X<//2  =  tfooW*+Mn2*t)  01.9) 

-  U<x>7{ cos2  82  +  *2sinJ02)  Of. 10) 


Note  Equation  4.9  and  Equation  4.10  represent  two  equations  for  two  unknowns,  namely 
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=  local  flow  velocity  vector 
U  2  =  Wire  1  and  Wire  2,  respectively 

=s  X-wire  bisector 

=  angle  subtended  by  l )  and  either  X-wire 
=  angle  subtented  by  Uw  and  li 
=  normal  to  tl'j 
=  normal  to  IFj 

=  angle  subtended  by  and  n) 

=  angle  subtended  by  (/w  ami  nj 
=  angle  subtended  by  D  and  the  reference  X-axis 


Figure  4.4.  Geometry  of  X-typc  hot-wire  probe  -  End  flow,  normal  configuration 


Vti  ami  0.  Solving  Equations  4.9  anti  4.10  for  U yields 

,  _ _ _ —Vi - —  (.j.H) 

cos3  Q\  +  k\  sin3  9\  cos3  #*  +  *2  sin3  #* 

From  Equations  4.9  and  4.0,  #j  and  9 *  are  recognized  as  complementary  angles. 

Hence,  cos3#*  =  sin3#*  and  sin3#*  =  cos3#j.  Substitution  of  these  values  into  Equation 
4.1 1  yields  with  some  manipulation 


U\ </;3(cos3  #*  *f  H**sin3  #*)  =  f/*<//3(sin3  #*  +  k\  cos3  #*)  (4.12) 

Dividing  Equation  4.12  by  cos3#*  and  solving  for  #* 


1  +  =  l',,„5(un  +  i|) 

-  »>,„’)  =  t'j,,,5*.'!  -  y.,„5 

iV*.  -  Vi.i 


8;  k  Inn"1 


Combining  Equation  4.S  with  Equation  4.13  gives 


. . easa 


(4.13) 


(4.14) 


The  correct  root  of  Equation  4.14  must  be  determined  from  the  geometry  of  the  How. 
From  Figure  4.4,  a  clockwise  rotation  from  the  bisector  (J?)  defines  a  positive  incidence 
angle  (0).  This  specific  flow  geometry  is  labeled  case  1.  Case  2  defines  a  counterclockwise 
rotation  from  the  bisector  with  a  negative  incidence  angle  (  -0). 

In  case  1,  V ,c//  >  l/2<// 

In  case  2,  tf2<//  >  Ult{/ 

The  term  within  the  parenthesis  of  Equation  4.14  represents  the  physics  of  the  inci¬ 
dence  angle  of  the  flow  as  it  approaches  each  wire.  The  incidence  angle  of  the  time  discrete 
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velocity  vector  indicates  the  term  within  the  parenthesis  of  Equation  4.14  must  he  positive 
if  >  Ui'fj  and  negative  if  V7,n  >  V\,ir  Equation  4.1*1  is  rewritten  to  include  the 
physics  of  the  flow  and  account  for  the  two  roots  of  Equation  4.14. 


0  =  o-  tan"1  bjit 

1  ,nk 

where 

sgn  =  +1 

if  U\,n  >  U7ttl 

i 

it 

!> 

if  V\,n  <  U7,n 

(4.15) 


Equation  4.15  is  solved  at  each  time  step.  The  X-wire  calibration  equations  provide 
V\,tf  and  during  each  sampling  interval,  and  o  is  known  for  each  X-wire;  therefore,  J 
is  computed  for  each  time  step.  Equations  4..r>  and  4.G  are  solved  for  0\  and  07  at  each  time 
step,  respectively.  Equation  4.7  or  Equation  4.8  is  now  solved  for  UK  at  each  time  step. 
In  essence,  the  time  discrete  local  velocity  is  determined,  where  both  angle  and  magnitude 
of  the  local  flow  is  known  at  each  time  step. 

If  A'  velocity  measurements  are  made  by  the  X-wire  probe  over  a  finite  time  interval, 
then  the  time  discrete  velocity  is  decomposed  into  local  velocity  components,  ii  and  v,  for 
each  time  step ». 

u,  =  U<x,  cos(#  -  A) 

e,  =  sin(/?;  -  A)  (4.10) 

=  the  angle  subtended  by  the  X-wire  bisector  and 
and  the  wind  tunnel  reference  X-axis.  Positive 
A  is  in  the  same  sense  as  positive  0 

The  mean  velocity  components  are  given  as 


where 

A 


r*  Vi 

TJ= 

N 


(•US) 


The  fluctuating  velocity  components,  u'  and  v',  arc  defined  for  each  lime  step, »,  as 

u'{  =:«;  —  ¥  (4.19) 

v'r^Vi-V  (4.20) 

The  root-mcan-squarc  of  the  fluctuating  velocity  components,  arc  given  as 

<«» 


*-(=£)* 


(4.22) 


The  x-  and  ^components  of  turbulence,  Tur  and  Tuv.  respectively,  are  defined  as 


Tu,[%)  =  '£  •  100% 

(4.20) 

'/«„(%)  =  |  •  100% 

(4.24) 

where 

V  -  the  vector  sum  of  17  4-  ? 

The  consolidated  turbulence  (7’u  )  is  computed  as 

n.  (JV+7V)i 

(4.25) 

4.5  Turbulence  Integral  Scale  and  Microscale  Measurements 

The  turbulence  integral  scale  and  microscale  lengths  were  determined  from  the  time- 
discrete  local  velocity  measurements.  The  method  of  analysis  follows  the  turbulence  scale 
discussion  presented  by  Cebecci  and  Smith  (Cebcci,  1974).  All  turbulence  scales  were 
resolved  to  the  tunnel  reference  coordinate  axis. 
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The  X-compoucnt  turbulence  integral  scale  (A„)  an<l  the  Y-componcnl  turbulence 
integral  scale  (Au)  arc  computed  as 


Au  =  lim  | 

n-0  ' 

(4.20) 

II 

ill 

a 

< 

&*<■>) 

(4.27) 

where 

n  =  turbulence  fre<picncy 

Ku'in)  -  X-componcnl  turbulence  power  per  frequency  bandwidth 
£v*(m)  =  Y-component  turbulence  power  per  frequency  bandwidth 
The  limit  is  taken  as  the  turbulence  frequency  («)  approaches  zero. 


The  X-componcnt  of  turbulence  microscale  (A„)  and  the  Y-component  of  turbulence 
microscale  (Au)  are  determined  as 


An 


f  Tr  «,2> 

i 

JT*’  «3£u'(")d». 

(rW) 

1 

A 2-* ) 

Jo . n'Gv.(n)d» 

01.28) 


01-29) 


where 


Nmar  =  highest  measured  frequency  of  the  turbulence 


The  turbulence  power  per  frequency  bandwidth  was  determined  from  the  Fast  Fourier 
Transform  (FFT)  of  the  fluctuating  velocity  components  (u'  and  o').  The  highest  frequency 
of  the  turbulence  used  for  the  microscale  calculations  in  this  study  was  arbitrarily  set  at 
19,000  Hz.  This  cut-oflT  frequency  eliminated  a  spike  in  the  FFT  data  at  20,000  Hz,  which 
was  attributed  to  instrumentation  noise.  In  addition,  during  the  experiment  it  became 
apparent  the  physical  limit  of  turbulence  frequency  in  this  tunnel  was  around  15,000  Hz. 
The  magnitude  of  the  frequency  components  is  nearly  zero  for  frequencies  beyond  13-15,000 
Hz. 
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A  program,  TUBAT.BAS  was  written  to  compute  the  turbulence  integral  scale 
lengths  and  microscale  lengths.  Program  TUBAT.BAS  is  included  in  Appendix  E. 

4-6  Surface  Beat  Tranfer  Measurement 

f.6.1  Turbine  Blade  Surface  Temperature  Measurement  The  heat  transfer  test  tur¬ 
bine  blade  was  heated  using  the  external  power  supply.  A  50  amp  current,  at  a  potential 
of  approximately  2.5  volts  was  conducted  through  the  stainless  steel  foil  surface.  The  tur¬ 
bine  blade  was  allowed  to  come  to  steady  state  conditions.  The  steady  state  temperature 
condition  on  the  blade  was  reached  in  approximately  five  minutes.  During  this  time,  tne 
temperature  of 'he  blade  was  constantly  monitored  by  the  HP3852A  Data  Acquisition  Sys¬ 
tem  in  order  to  prevent  overheating  of  the  turbine  blade  surface.  Overheating  will  cause 
the  underlying  urethane  foam  to  out-gas  and  separate  the  foil  surface  from  the  urethane 
core. 

Once  steady  state  conditions  were  reached,  the  HP3S52A  Data  Acquisition  System 
read  each  thermocouple  voltage  and  converted  the  voltage  to  temperatures.  A  single  scan 
of  the  thermocouples  was  completed  by  the  HP3852A  in  less  than  0.3  seconds.  Prior  to 
making  a  temperature  scan  of  the  thermocouples,  the  external  power  supply  which  heats 
the  turbine  blade  surface  was  momentarily  shut  ofi";  otherwise,  the  current  would  short 
through  the  thermocouple  leads  into  the  FET  multiplexer  and  damage  the  PET  junction. 
After  reading  all  thermocouples,  the  HP3S52A  turned  the  power  supply  back  on.  The 
entire  process  of  turning  the  power  ofT,  reading  the  thermocouples,  and  turning  the  power 
supply  back  on,  was  completed  in  about  0.5  seconds.  A  change  in  the  turbine  blade  surface 
temperature  was  not  detectable  until  1.3  seconds  after  the  power  supply  was  shut  ofT;  hence 
turning  the  power  supply  off  momentarily  had  no  effect  on  the  temperature  measurements. 

The  temperature  readings  for  each  test  configuration  listed  in  Tables  1.1  and  1.2  were 
stored  on  disk  for  post-processing.  Subroutine  TRUN.BAS  of  the  main  data  acquisition 
code  was  used  to  control  the  HP3852A  during  temperature  measurements.  A  source  code 
listing  of  TRUN. BAS  is  given  in  Appendix  E. 
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4.6,2  Stanton  and  Nusselt  Numbers  on  the  Diode  The  Stanton  number  (St  )  ami 
Nusselt  number  (Afu  )  for  the  turbine  blade  surface  heat  transfer  are  computed  from  a 
power  balance  on  tlte  turbine  blade  surface.  The  turbine  blade  surface  is  considered  as  a 
scries  of  adjoining  sections.  The  sum  total  of  all  the  i  surface  sections  make  up  the  total 
surface  area  of  the  heat  transfer  surface  of  the  turbine  blade.  At  each  i  tii  section  of  the 
turbine  blade  surface,  a  specific  power  balance  was  applied  as 


9t/«f(rie«/,  =  QcanJuctioni  "b  ^convection,-  9ro<fi«(ion, 


(•t.30) 


and 


where 

titled  neat, 

tlconduttion, 

^convection, 

(/radiation, 

1 

Rt 

kfoam 

Ai 

Li 

Tint, 


tJdtetrtent, 

t/cauJudwn, 

Veonvrdioni 

(/radiation, 


i2n' 


•i, 


(T,-r 


hjoamA, 

L, 

h,Ai(T,t-T^) 


■  Ml, 


.) 


(•1.31) 


=  electrical  power  input  to  each  blade  section 
=  conductive  heat  loss  at  each  section 
=  convective  heat  loss  at  each  section 
=  radiative  heat  loss  at  each  section 
=  electrical  current  (amperes) 

=  electrical  resistance  of  the  total  foil  surface  area  (ohms) 

=  thermal  conductivity  for  urethane  foam  (0.020  ) 

=  area  of  the  itli  section  of  the  turbine  blade  foil  surface 
=  linear  distance  from  the  surface  of  the  blade  to  the  interior 
point  in  the  turbine  blade  where  the  temperature  is  measured 
=  section  surface  temperature 

=  internal  temperature  of  the  turbine  blade  corresponding 
to  each  surface  area  section 
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=  convection  heat  transfer  coonicient  for  each  section 
ss  Stcfan-Holtznunn  constant  (5.67  x  10”8  W/m  •  K) 
=  cmissivity  of  the  stainless  steel  foil  (0.17) 


h, 

a 

c 


The  thertna!  conductivity  value  is  for  rigid  urethane  foam  at  300  K;  the  cmissivity  value 
is  for  highly  polished  stainless  steel  foil  at  300  K  (Incropcra,  1985). 

Substituting  the  set  of  Equations  4.31  into  Equation  4.30,  canceling  the  section  area 
(/t<)  from  all  terms,  ami  solving  for  /»,  gives 

+  Vir.  -  rM.)  -?v) 

hi  = - - - - -  (4.32) 

“  1  <s> 

Using  /»,,  the  Nussclt  number  obtained  for  each  blade  section. 


where 

•s,  =  distance  along  the  surface  of  the  blade  from  the  leading  edge 

of  the  blade  to  the  specific  thermocouple  location  within  area  i 
Aw,  =  thermal  conductivity  of  air  evaluated  at  the  mean  fluid 
temperature  (  Tmtan  ) 


T  - 

1  meant 


(4.34) 


The  Stanton  number  (St)  for  each  blade  section  is  computed  as 


St{  = 


Nu; 

PrRe„ 


(4.35) 


where 

Pr  =  Prandtl  number  of  air 

Rc3i  =  Reynolds  number  based  on  surface  distance  s 


4-15 


Table  *1.2.  Experimental  measurement  uncertainties 


ETjxrimcnial  Reading  | 

variable 

range 

variation 

uncertainty 1  %  ) 

length  measurements 

— 

— 

±0.50 

mm 

0.50 

manometer  ( P0  -  /'«) 

13.71 

cm  //jO 

±  0.13 

cm 

0.95 

injection  pressure 

20G.85 

MPa 

±  1.379 

MPa 

0.07 

blade  surface  temperature  (T) 

mm 

mm 

2. GO 

T«, 

27.50 

•c 

EM 

3.GI 

1 

50.00 

amp 

±0.2 

amp 

0.40 

resistance  (  R  ) 

0.13 

n 

0.001 

n 

0.70 

emissivity  (<) 

0.17 

— 

1.80 

Pr 

— 

0.10 

thermal  conductivity  (  A*„fr  ) 

0.02G3 

W/mK 

— 

0.10 

thermal  conductivity  (  k/Mm  ) 

0£2G 

W/mK 

— 

_ 

0.10 

Stanton  au<l  Nussclt  numbers  were  computed  from  the  blade  temperatures  using 
subroutine  TIIUN.RAS of  the  main  data  acquisition  program  listed  in  Appendix  E. 

./.7  Er/Ki-hncnlnl  and  Calculated  Uncertainties 

Measurement  uncertainties  are  arc  listed  in  Table  *1.2.  An  example  of  an  uncertainty 
calculation  for  a,  measurement  is  shown  for  the  turbine  blade  surface  temperature  (  T  ). 

[n  oM » 

unecrtuintxj{T)  -  "  ’JJ  -  =  2.G0 % 

*JO»0 

Instrumentation  accuracy  is  given  in  Table  4.3. 
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Table  4.3.  Instrument  accuracy  specifications 


Instrument  Accuracy  \ 

instrument 

range 

uncertainty  (  7o  ) 

HIM4701  Voltmeter 

30  mV 

0.010 

Hl’44701  Voltmeter 

3  V 

O.OOS 

HP44702H  Voltmeter 

2.5G  V 

0.005 

1IP44713  Thermocouple  (  T  ) 

30  °C 

0.333 

1IP3S52A  internal  Clock 

10  /is 

0.75 

TSI  I FA- 100 

BU 

0.02 

Scani valve  transducer 

2500  k!*a 

0.8G 
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Using  the  uncertainties  listed  above,  the  uncertainties  of  calculated  variables  were 
computed.  Table  4.4  lists  the  uncertainties  of  the  calculated  variables.  An  example  of  a 
calculated  uncertainty  is  given  for  the  exit  Reynolds  number  (  Hec  ). 


&Kcf  =  [(1.0)a  +  (I.31)2  4-  (0.5)a]  ^  =  1.72% 


Table  4.4.  Uncertainty  of  calculated  variables 


r.  RESULTS  AND  DISCUSSION 


$.!  Effect  of  the  Jet-GrUl  on  Turbine  Wade  Surface  Pressures 

The  effect  of  (he  Jet-grid  on  the  surface  pressure  distributions  was  examined  for  co¬ 
flow,  cross-flow,  and  counterflow  injection  of  the’  secondary  flow  over  a  range  of  injection 
pressures. 

A  baseline  case  was  established  with  the  surface  pressures  measured  on  the  turbine 
blade  without  any  turbulence  grid  installed  in  the  free-stream  flow.  The  surface  pressures 
are  presented  as  nondimcnsional  coefficients  of  pressure  (Cr)  along  the  nondimensinnalifccd 
chord  (x/c)  of  the  turbine  blade. 

Figure  5.1  shows  the  pressure  distribution  on  the  turbine  blade  surface  over  a  range  of 
the  wind  tunnel  exit  Reynolds  numbers.  The  surface  pressure  distribution  determined  by 
Langston  (Langston,  1977)  on  an  identical  turbine  blade  profile  is  included  for  comparison. 
Figure  5.1  demonstrates  the  independence  of  the  pressure  coefficient  from  the  exit  Reynolds 
number.  The  pressure  side  of  the  turbine  blade  exhibits  a  stagnation  point  (C',,  =  1)  at  the 
0.7  percent  chord.  Along  the  pressure  surface  of  the  turbine  blade,  the  pressure  coefficients 
measured  compared  favorably  with  Langston's  data  except  near  the  trailing  edge  of  the 
turbine  blade.  The  suction  surface  pressure  coefficients  exhibited  significant  deviations 
from  Langston's  results,  most  noticeably  in  the  pressure  distribution  measured  near  the 
leading  edge  of  the  suction  surface.  The  free-stream  incidence  angle  was  the  same  for  both 
studies;  therefore,  this  difference  in  pressure  distribution  between  the  present  study  and 
Langston’s  results  is  attributed  to  the  secondary  flow  effects  within  the  cascade,  and  the 
actual  location  of  the  pressure  port  taps  on  the  turbine  blade. 

Studies  of  others  (Moore,  19S*lc,  Langston,  1977)  indicate  strong  secondary  flow 
effects  are  present  within  linear  turbine  cascades.  A  strong  passageway  vortex  impinges  on 
the  suction  surface  and  tends  to  decrease  the  pressures  measured  on  the  suction  surface. 
The  secondary  flow  effects  are  more  pronounced  for  low  aspect  ratio  turbine  cascade  wind 
tunnels.  The  present  study  used  a  turbine  cascade  with  a  unity  aspect  ratio,  thus  secondary 
flow  effects  arc  expected.  Also,  because  the  turbine  blade  used  in  the  present  study  was 
less  than  half  the  size  of  Langston’s  model,  the  secondary  flow  within  a  linear  turbine 
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Figure  5.1. 


Pressure  distributions  on  the  turbine  blade  for  several  exit  Reynolds  numbers. 
No  turbulence  grid  installed 


cascade  is  expected  to  have  a  more  pronounced  effect  on  the  smaller  scale  turbine  cascade 
and  the  turbine  blade  surface  pressures  would  be  affected. 

For  all  remaining  turbine  blade  surface  pressure  measurements  and  pressure  coefli- 
clout  calculations,  the  free-stream  total  pressure  and  dv itamic  pressure  were  measured  with 
the  pitot  tube  in  a  different  position  in  the  wind  tunnel  than  the  position  used  to  measure 
the  free-stream  total  and  static  pressures  for  the  results  shown  in  Figure  5.1.  Relocation 
of  the  pitot  probe  was  necessary  in  ordr*r  to  accommodate  the  jet-grid  device.  The  new 
position  of  the  pitot  probe  was  midway  between  turbine  blades  Number  2  and  Number  8. 
15.0  mm  upstream  of  the  reference  plane,  and  midway  between  the  two  tunnel  cudwalls. 
Since  the  free-stream  accelerated  as  it  entered  the  cascade  passageway,  the  measured  free- 
stream  dynamic  pressure  (//»)  was  greater  at  this  pitot  probe  location  than  the  free-stream 
dynamic  pressure  measured  at  the  previous  location  where  tiie  results  siiown  in  Figure  5.1 
were  obtained.  Due  to  the  different  dynamic  pressures  associated  with  the  two  different 
pitot  probe  locations  in  the  wind  tunnel,  the  no-grid  baseline  pressure  distribution  shown 
in  Figures  5.2  •  5.8  is  not  the  same  as  the  no-grid  baseline  pressure  distribution  shown  in 


Figure  5.1.  However,  the  results  shown  in  Figures  5.1  -  5.S  have  the  same  exit  Reynolds 
number. 

The  pressure  distribution  on  the  turbine  blade  with  the  turbulence  grid  installed  at 
two  upstream  locations,  with  no  blowing,  is  shown  in  Figure  5.2.  The  pressure  distribution 
on  the  turbine  blade  surface,  without  the  turbulence  grid  installed,  is  included  for  com¬ 
parison.  On  the  pressure  surface,  the  pressure  coefficient  remained  constant.  The  suction 
surface  pressure  coefficients  decreased  with  the  jet-grid  installed,  especially  near  the  trail¬ 
ing  edge  of  the  blade.  Subsequent  heat  transfer  tests  and  comparison  to  laminar  boundary 
layer  heat  transfer  theory  and  turbulent  boundary  layer  heal  transfer  theory  indicated 
the  higher  turbulence  intensity  generated  by  the  turbulence  grid  in  the  free-stream  flow 
induced  a  turbulent  boundary  layer  on  the  suction  side  which  was  better  able  to  negotiate 
the  curvature  of  the  suction  surface. 

The  free-stream  incidence  angles  were  measured  with  the  jet-grid  installed  in  the 
flow.  For  the  no  injection  cases  shown  in  Figure  5.2  there  was  a  possibility  of  some  small 
incidence  angle  change;  however,  it  was  too  small  of  a  change  to  accurately  measure. 
Incidence  angle  changes  less  than  0.75  degrees  were  within  the  uncertainty  of  the  incidence 
angle  measurements. 

Figures  5.3, 5.*t.  and  0.5  show  the  effect  on  surface  pressure  distributions  on  the  tur¬ 
bine  blade  with  jet.grid  injection  in  different  directions,  over  a  range  of  injection  pressures 
(«  75-300  kPa).  The  same  location  upstream  of  the  cascade  was  used  for  the  results  of 
tlmso  three  figures  (£•  =  15.75). 

Doth  co-flow  and  cross-flow  exhibited  very  little  effect  on  the  pressure  distributions 
on  the  turbine  blade  compared  to  the  case  with  no  injection.  However,  the  counter¬ 
flow  injection  case,  shown  in  Figure  5.5,  had  a  significant  effect  overall  on  the  turbine 
blade  pressure  distribution.  With  counter-flow  injection,  the  pressure  coefficient  increased 
slightly  over  most  of  the  turbine  blade  pressure  surface  as  compared  to  the  no  injection 
case  (solid  line).  The  pressure  coefficient  seems  to  approach  stagnation  levels  around  | 
=  0.2-0.3.  A  slight  decrease  in  Cp  occurred  near  the  trailing  edge  of  the  pressure  surface 
which  indicates  the  flow  was  accelerated  in  this  region.  Conversely,  counter-flow  injection 
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Figure  5.2. 


Pressure  distribution  on  the  turbine  blade  with  the  jet-grid  installed  at  tvvu 
upstream  locations.  No  secondary  How  injected.  llce  -  5.0  x  10s 
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Figure  5.3.  Pressure  distribution  on  the  turbine  blade  with  co-flow  injection  at  the  near¬ 
upstream  location.  ]ice  =  5.0  x  10s 
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Figure  5.4.  Pressure  distribution  on  the  turbine  blade  with  cross* flow  injection  at  the 
near-upstream  location.  ltce  -  5.0  x  10* 
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Figure  5.5.  Pressure  distribution  on  the  turbine  blade  with  counter-flow  injection  at  the 
near-upstream  location.  Rce  =  5.0  x  10s 


decreased  the  pressure  coeflicent  on  the  blade  suction  surface  dramatically.  The  injection 
direction  was  the  dominant  factor  in  changing  the  turbine  blade  pressure  distribution.  The* 
injection  pressure  had  less  effect  on  the  pressure  distribution  on  the  turbine  blade  than 
the  injection  direction.  The  flow  incidence  angle  was  measured  for  this  counter-flow  case 
and  found  to  have  changed  from  -1*1.7  degrees  to  43.0  degrees  which  might  account  for  the 
decrease  in  pressure  coefficients  on  the  suction  surface.  The  Row  incidence  angle  for  the 
other  two  injection  directions  remained  at  44.7  degrees.  The  effect  of  counter-flow  injection 
on  the  secondary  vortex  core  within  the  cascade  passageway  is  not  known,  but  some  effect 
is  suspected.  This  area  deserves  further  study. 

With  the  .jet-grid  placed  at  the  far-upstream  location  (f  =  25.67),  the  pressure 
distributions  on  the  turbine  blad ,  were  similar  to  the  pressure  distributions  on  the  turbine 
blade  with  the  jet-grid  placed  at  station  0.  The  effect  on  the  blade  pressure  distribution 
for  co-flow,  cross-flow,  and  counterflow  injection  at  this  upstream  location  is  displayed  in 
Figures  5.G,  5.7,  and  5.S,  respectively.  For  all  three  injection  directions  tested  at  station  02. 
the  effect  of  jet-grid  injection  on  the  pressure  coefficient  was  confined  to  the  suction  surface 
only.  The  pressure  distributions  shown  in  Figures  5.6,  5.7,  and  5.S  are  nearly  identical. 
At  this  jet-grid  position,  counter-flow  injection  shown  in  Figure  5.S  did  not  have  a  same 
effect  as  counter-flow  injection  at  station  0  shown  in  Figure  5.5.  This  result  became  very 
important  during  subsequent  turbine  blade  heat  transfer  measurements. 

5.2  Turbulence  intensity  and  Decay 

The  turbulence  intensity  of  the  free-stream  flow  was  measured  for  all  test  config¬ 
urations.  Equations  4.23,  4.24,  and  4.25  were  used  to  compute  the  X-component,  V- 
component,  and  the  consolidated  turbulence  intensity,  respectively.  The  turbulence  inte¬ 
gral  length  scales  (Au,  A„)  and  turbulence  microscales  (Au,  Au)  were  computed  from  the 
turbulence  energy  spectrum  of  the  free-stream  flow.  Figure  5.9  shows  the  turbulence  en¬ 
ergy  spectrum  (Eu<)  as  a  function  of  turbulence  frequency  for  the  free-stream  turbulence 
measured  along  the  X-axis  without  the  jet-grid  installed  (Test  XXX).  Figure  5.10  shows 
the  Y-axis  component  of  the  spectrum  of  turbulence  energy  (£„»)  for  the  same  test.  The 
vonKarman  spectrum  (Batchelor,  1953)  of  turbulence  predicts  a  -5/3  slope  of  the  turbu- 
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Pressure  distribution  on  the  turbine  blade  with  co-flow  injection  at  the*  far* 
upstream  location,  //c,-  =  5.0  .\  105 


0.0  0.2  0.4  0.6  0.5  t.O 


3.0 


2.0 


1.0 


0.0 


-1.0 


-2.0 


-3.0 


-4.0 


-5.0 


Figure  5.7.  Pressure  distribution  on  the  turbine  blade  with  cross-flow  injection  at  the 
far-upstream  location.  Rce  =  5.0  x  10s 
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Figure  5.8.  Pressure  distribution  on  the  turbine  blade  with  counter-flow  injection  at  the 
far-upstream  location.  Rce  =  5.0  x  105 


lwicc  energy  spectrum  if  the  free-stream  turbulence  is  homogeneous.  The  -5/3  slope  is 
shown  as  a  solid  line  in  Figure  5.9.  Both  the  X-axis  and  the  Y-axis  component  of  the 
turbulence  energy  spectrum,  Figures  5.9  and  5.i0,  respectively,  were  too  scattered  to  draw 
any  conclusions  as  to  the  fit  of  the  theoretical  prediction.  This  is  an  indication  that  the 
turbulence  was  not  uniform  when  the  jet  grid  was  not  placed  into  the  free-stream. 

The  turbulence  energy  spectrum  for  the  jet-grid  at  station  0  (f  =  15.75)  ,  without 
secondary  injection  (Test  000),  is  shown  in  Figures  5.1!  and  5.12.  The  turbulence  energy 
spectrum  for  the  jet-grid  placed  at  station  02  (|  =  25.G7),  without  secondary  injection 
(Test  002),  is  shown  in  Figures  5.13  and  5.14.  With  the  jet-grid  in  place  at  cither  station, 
the  slope  of  the  turbulence  energy  spectrum  matched  the  -5/3  slope  very  well,  which 
indicated  homogeneous  free-stream  turbulence  was  achieved. 
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Figure  5.10.  v’-turbulencc  energy  spectrum— Test  XXX.  Theory  (Batchelor,  1953) 
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Figure  5.12.  v'-turbulenco  energy  spectrum— Test  000.  Theory  (Batchelor,  1953) 
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Figure  5.13.  i/'-turhulence  energy  spectrum— Test  002.  Theory  (Batchelor,  1953) 


Figure  5.14.  v'-turbulence  energy  spectrum— Test  002.  Theory  (Batchelor,  1953) 
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Plots  of  turbulence  energy  spectrum  for  all  other  test  configurations  arc  given  in 
Appendix  B.  For  all  other  test  configurations,  the  turbulence  energy  spectrum  was  well 
correlated  to  the  theoretical  prediction.  This  is  an  indication  of  consistency  of  the  turbu¬ 
lence  measurement  process  and  the  uniformity  of  the  turbulence  in  the  free-stream. 

The  results  of  the  turbulence  measurements  arc  tabulated  in  Tables  5.1  and  5.2.  In 
Table  5.1,  the  test  designators  A,  B,  and  C  indicate  co-flow,  cross-flow,  and  counter-flow 
injection  at  station  0,  respectively.  In  Table  5.2,  test  designators  D,  E.  and  F  indicate  co¬ 
flow,  cross-flow,  and  counter-flow  injection  at  station  02,  respectively.  Note  the  extremely 
low  turbulence  level  measured  during  test  AM  as  compared  to  other  tests  at  the  same 
station  and  the  ''ory  large  value  of  A»  as  compared  to  the  value  of  A„.  The  results  of 
Test  AM  are  not  included  in  the  following  discussion  of  turbulence  intensity  since  sub¬ 
sequent  testing  indicated  the  hot-wire  probe  was  not  measuring  free-stream  conditions, 
but  was  being  impinged  upon  directly  by  a  jet  of  injected  flow.  Also,  for  Tests  F21,  F22, 
and  F2.'t;  the  Y-component  integral  scales  and  Y-component.  microscales  were  virtually 
indistinguishable. 

The  turbulence  intensity  was  plotted  against  the  jet-grid  injection  pressure  for  co. 
flow,  cross-flow,  and  counter-flow  injection  at  both  stations.  Figure  5.15  shows  the  variation 
of  turbulence  intensity  for  co-flow  injection  at  both  jet-grid  locations  in  the  free-stream 
flow.  The  higher  levels  of  turbulence  intensity  were  measured  for  the  jet  grid  at  station  0 
and  the  lower  turbulence  levels  were  measured  at  station  02.  Two  conclusions  arc  made 
from  Figure  5.15.  First,  the  turbulence  intensity  was  relatively  constant  over  the  range 
of  injection  pressures  tested  for  co-flow  injection  for  all  three  jet-grid  orifice  diameters. 
Second,  the  greatest  variation  in  turbulence  intensity  was  due  to  the  location  of  the  jet- 
grid  in  the  free-stream  flow  with  respect  to  the  leading  edge  of  the  turbine  cascade. 

The  turbulence  intensity  variation  for  cross-flow  and  counter-flow  test  configurations 
are  shown  in  Figure  5. 1C  and  Figure  5.17,  respectively.  From  Figures  5.15,  5.16,  and  5.17, 
the  most  noticable  trend  is  that  secondary  injection  did  not  have  a  large  effect  on  the 
free-stream  turbulence,  with  cross-flow  injection  (Figure  5.1C)  having  the  least  effect  on 
the  free-stream  turbulence  of  all  the  injection  directions  testou.  The  results  of  cross-flow 
injection  are  extremely  important  in  that  the  same  turbulence  l-.vel  is  maintained  over  the 
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Table  5.1.  Turbulence  intensity  and  turbulence  scale  results  for  turbine  cascade  wind 
tunnel  test  configurations— station  0 
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Table  5.2.  Turbulence  intensity  and  turbulence  scale  results  for  turbine  cascade  wind 
tunnel  test  configurations— station  02 
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Figure  5.15. 


Effect  of  jet-grid  injection  pressure  on  free-stream  turbulence  intensity  for 
co- flow  injection 


range  of  injection  pressures,  but  from  Tables  5.1  and  5.2,  for  tests  designated  with  either 
a  II  or  an  E  (cross-llow  injection),  the  turbulence  microscale  lengths  are  not  the  same  and 
to  a  lesser  degree,  the  turbulence  integral  scale  lengths  arc  not  the  same.  Hence  cross-flow 
injection  allows  one  to  observe  the  effect  of  turbulence  scale  lengths  on  turbine  blade  heal 
transfer  measurements  at  a  constant  turbulence  level.  Further  discussion  on  this  effect  is 
delayed  until  heat  transfer  results  arc  discussed. 

The  variation  of  turbulence  intensity  with  injection  pressure  is  shown  in  Figures  5.15, 
5.1G,  and  5.17  for  co-flow,  cross-flow,  and  counter-flow  injection,  respectively.  The  lower 
turbulence  levels  arc  the  result  of  jet-grid  injection  at  station  02,  while  the  higher  levels 
of  turbulence  arc  the  result  of  injection  at  station  0.  From  Figures  5.15  and  5.1G,  the 
turbulence  intensity  was  relatively  constant  for  all  injection  pressures,  and  for  all  injection 
orifice  diameters,  using  co-flow  and  cross-flow  injection.  However,  Figure  5.17  shows  the 
frec-stream  turbulence  intensity  increased  1-2  percent  with  increasing  injection  pressure 
for  both  1.32  mm  and  1.98  mm  injection  orifice  diameters,  but  only  with  the  jet-grid  at 
station  0.  A  1-2  percent  change  in  turbulence  intensity  represents  a  small  change  in  the 
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Figure  5.17.  Effect  of  jet-grid  injection  pressure  on  free-stream  turbulence  intensity  for 
counter-flow  injection 
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overall  turbulence  level,  and  for  the  most  part,  the  injection  pressure  seemed  to  have  a  small 
effect  on  turbulence  intensity.  This  result  is  in  contrast  to  the  results  of  Young  and  Han 
(Young,  1989)  where  they  determined  jet-grid  injection  increased  the  turbulence  intensity 
of  the  free-stream  flow  by  an  amount  of  up  to  45  percent  over  the  non-injection  case.  The 
difference  in  results  may  be  due  to  several  factors.  Voting  and  Han  used  supersonic  jot-grid 
injection,  while  this  study  was  limited  to  sonic  injection,  except  during  test  A 13,  where 
sonic  injection  is  thought  to  have  occurred.  Young  and  Han  (Young,  19S9)  found  that  the 
turbulence  level  remainied  relatively  constant  until  they  used  supersonic  injection.  Only 
with  supersonic  injection  were  they  able  to  obtain  the  large  increase  in  turbulence  intensity. 

An  analysis  of  the  free-stream  turbulence  decay  is  presented.  The  turbulence  iitten»il> 
decay  behind  the  jet-grid  at  station  9  (j  =  10.75)  compares  favorably  with  the  turbulence 
decay  correlation  used  by  Hlair  (Hlair,  1983a).  The  free-stream  turbulence  level  varied 
from  8.9  percent  using  the  Q.GG  mm  diameter  injection  orifice  and  co-flow  injection  (Test 
AIM),  up  to  11.8  percent  with  the  jet-grid  at  station  0  with  no  jet-grid  injection  (Test  000). 
The  correlation  used  by  Hlair  (Hlair,  1983a)  is  given  by  Equation  2.12.  Using  Equation 
2.12,  the  free-stream  turbulence  level  at  a  non-dimensional  distance  of  £  =5  25.07  behind 
the  jet-grid  is  computed  as  10.9  percent  which  closely  matched  the  results  of  the  present 
study. 

At  station  02  (j  =  25.07),  the  free-stream  turbulence  level  varied  from  *1.8  percent 
using  the  1.32  mm  injection  orifice  and  cross-flow  injection  (Test  E21),  up  to  0.0  percent 
using  the  0.GC  mm  injection  orifice  and  co-flow  injection  (Test  Dll).  Equation  2.12  predicts 
a  7.7  percent  turbulence  level.  This  represents  a  10  percent  difference  between  computed 
and  measured  turbulence  levels.  One  reason  for  the  difference  between  the  computed 
and  measured  turbulence  levels  may  be  different  test  conditions  existed  between  Hlair’s 
experiment  and  the  present  experiment.  Blair’s  correlation  was  developed  from  tests  on 
a  flat  plate  where  the  wind  tunnel  free-stream  velocity  was  constant  along  the  length  of 
the  tunnel.  In  contrast,  in  the  linear  turbine  cascade  wind  tunnel  of  this  study,  the  free- 
stream  accelerated  as  the  flow  approached  the  leading  edge  of  the  turbine  cascade  where 
the  turbulence  intensity  was  measured  with  a  hot-wire  probe. 

Turbulence  is  proportional  to  the  fiuctating  component  of  velocity  (t/,  1/).  As  the 
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free- stream  velocity  increases,  the  fluctuating  component  of  that  velocity  remains  relatively 
constant  or  even  decreases  under  accelerating  flow  conditions,  while  the  mean  velocity  in* 
creases  and  the  measured  turbulence  intensity  decreases.  The  difference  between  acceler* 
aiing  frcc-strcam  and  steady  frcc-stream  conditions  may  account  for  the  difference  in  the 
predicted  turbulence  level  from  Equation  2.J2  and  the  actual  turbulence  level  measured 
with  the  jet-grid  at  station  02. 

5.2.1  1  urbutcncc  Integral  Scale  The  X-componcnt  integral  scale  (A„)  length  was 
on  the  same  order  as  the  outside  diameter  of  the  jet-grid  injection  tube  (6  =  0.035  mm) 
(See  Figures  5.1$,  5.10,  and  5.20.  As  seen  in  Figures  5.21, 5.22,  and  5.23,  the  Y-component 
integral  scale  (A,.)  was  an  order  of  magnitude  smaller  than  the  X-componeni  integral  scale, 
except  for  the  case  where  it  is  thought  that  a  jet  of  injected  flow  impinged  directly  on  the 
hot-wire  anemometer  and  distorted  the  test  results  (Test  AM). 

The  variation  of  integral  scale  length  was  examined  for  dependence  on  jet-grid  orifice 
diameter,  injection  pressure,  and  injection  direction.  Figures  5.1$,  5.19,  and  5.20  show  ef¬ 
fect  on  the  X-componcnt  integral  scale  for  the  thret?  orifice  diameters  tested  over  the  same 
range  of  injection  pressures  for  co-flow,  cross-flow,  and  counter-injection,  respectively,  at 
station  0  (|  =  15.75).  The  results  of  Figures  5.1$,  5.19,  and  5.20  were  repeatable  for  sim¬ 
ilar  frce-stream  conditions.  Uncertainly  in  defining  the  zero  frequency  turbulence  energy 
(/♦„*(>»)  as  ii  —  0)  was  minimimized  by  averaging  all  turbulence  energies  between  0  and 
1000  Hz  frequencies;  this  method  provided  consistent  results  for  the  computation  of  the 
turbulence  integral  scale  lengths  listed  in  Tables  5.1  and  5.2.  The  turbulence  scale  mea¬ 
surements  (both  integral  scale  end  microscale)  had  a  calculated  uncertainty  of  2  percent. 

Comparing  the  results  given  in  Figures  5.18, 5.19,  and  5.20,  the  X-componcnt  integral 
scale  remained  nearly  constant  over  the  range  of  injection  pressures.  The  only  exception 
was  a  30  percent  increase  in  integral  scale  length  which  occurred  for  co-flow  injection  at  an 
injection  pressure  of  approximately  300  *  ’’a,  Figure  5.18.  From  Figure  5.20,  counter-flow 
injection  increased  slightly  the  X-componci.t  integral  scale  over  the  case  with  the  jet-grid  in 
place,  but  without  injection.  The  larger  diameter  injection  orifices,  generally  produced  the 
larger  X-component  integral  scales.  Young  and  Han  (Young,  1989)  propose  counter-flow 
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injection  increases  the  elective  diameter  of  the  jet-grid  tube  that  the  frcc-stream  sees  as  it 

flows  through  the  grid.  Ttcse  results  support  this  proposal  since  the  larger  X-component 

integral  scale  lengths  were  obtained  for  the  counter-flow  configuration. 

* 

Figures  5.21,  >.22,  and  5.23  show  the  efTcct  on  the  Y-component  integral  scale  for  the 
three  orifice  diameters  tested  over  the  same  range  of  injection  pressures  for  co-flow,  cross¬ 
flow,  and  counter-injection,  respectively,  at  station  0.  The  Y-component  integral  scales 
were  nearly  independent  of  injection  orifice  diameter.  The  Y-component  integral  scale  also 
seemed  insensitive  to  injection  pressure  except  for  one  case.  Figure  5.21  shows  a  dramatic 
increase  (  S00  percent)  of  integral  scale  size  for  the  0.65  mm  diameter  injection  orifice  at 
about  300  kPa.  As  discussed  previously,  this  datum  represents  the  results  of  Test  AM 
where  supersonic  injection  velocities  may  have  occurred.  Since  injection  pressures  were 
measured  and  not  flow  rates,  it  is  not  possible  to  accurately  determine  the  exit  velocity. 
However,  the  high  jet-grid  injection  pressures  make  it  highly  probable  that  the  flow  was 
supersonic  during  Test  AM. 

Similar  to  the  results  for  the  X-componcnt  and  Y-component  integral  scales  with  the 
jet-grid  at  station  0,  results  of  the  X-component  and  Y-component  integral  scales  were 
obtained  for  the  jet-grid  located  at  station  02  (£  =  25.67).  Plots  of  the  X-component 
integral  scale  versus  injection  pressure  for  all  three  diameters  are  given  in  Figures  5.2*1, 
5.25,  and  5.2G  for  co-flow,  cross-flow,  anti  counter-flow,  respectively.  The  Y-component  of 
integral  scales  measured  over  the  injection  pressure  range,  with  the  jet-grid  at  station  02, 
are  shown  in  Figures  5.27, 5.2S,  and  5.29  for  co-flow,  cross-flow,  and  counter-flow  injection 
directions,  respectively. 

Most  noticeable  from  Figures  5.2-1, 5.25,  and  5.2C  is  the  result  that  the  X-component 
integral  scale  length  was  on  the  order  as  the  jet-grid  tube  diameter  (0.635  mm).  This  is 
the  same  result  obtained  with  the  jet-grid  at  station  0.  Hence,  the  turbulence  integral 
scale  did  not  decay  downstream  of  the  jet-grid.  In  contrast  to  the  results  obtained  for  the 
integral  scale  length  with  the  jet-grid  at  station  0,  Figure  5.20,  counter-flow  injection  at 
station  02,  shown  in  Figure  5.26,  did  not  increase  the  X-component  integral  scale  length. 
The  efTcct  ofinjection  orifice  diameter  on  X-component  integral  scale  length  was  the  same 
as  for  the  jet-grid  placed  at  cither  station.  On  the  other  hand,  the  Y-component  integral 
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scale  length  became  more  dependent  on  Injection  orifice  diameter  than  with  the  jet-grid 
at  station  0.  A  larger  variation  of  the  Y-component  integral  scale  corresponding  to  each 
orifice  diameter  is  visible  in  Figures  5.27,  5.28,  and  5.29  as  compared  to  the  variation  of 
the  Y-componcnl  integral  scale  obtained  with  the  jet-grid  at  station  02  shown  in  Figures 
5.21,  5.22  and  5.23. 
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Figure  5.  IS. 


Effect  of  injection  pressure  ami  orifice  diameter  on  X-component  integral 
scale  length  for  co-flow  injection  at  £  =  15.75 


The  effects  of  jet-grid  injection  orifice  diameter  on  the  X-component  and  Y-component 
integral  scale  were  cross-plotted  against  their  respective  free-stream  turbulence  level  com¬ 
ponent.  Figures  5.30-5.32  present  the  X-component  turbulence  integral  scale  lengths  as 
a  function  of  the  respective  X-component  of  turbulence  intensity  (Tur)  for  each  of  the 
three  jet-grid  injection  diameters.  Each  plot  shows  the  integral  scale  for  all  injection  pres¬ 
sures,  for  the  jet-grid  at  both  locations  in  the  free-stream  flow.  The  injection  directions 
are  denoted  by  either  a  solid  symbol  for  co-flow  injection,  a  hatched  symbol  for  cross-flow 
injection,  or  a  open  symbol  for  counter-flow  injection.  The  X-component  integral  scale 
remained  relatively  independent  of  the  free-stream  turbulence  intensity.  Figures  5.30-5.32 
show  the  same  integral  length  scale  was  obtained  throughout  the  range  of  turbulence  lev- 
els  tested  (5-13.9  percent),  except  in  the  case  of  co-flow  injection,  at  the  10-12  percent 
turbulence  level,  shown  in  Figure  5.32.  In  this  case,  the  solid  star  symbol  represents  the 
datum  obtained  from  Test  AM,  where  it  is  thought  supersonic  injection  through  the  jet- 
grid  occured.  All  other  results  are  for  subsonic  injection  through  the  jet-grid.  Comparison 
between  each  of  these  three  figures  shows  the  injection  orifice  diameter  had  no  distinct 
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Figure  5.19.  Effect  of  injection  pressure  and  orifice  diameter  on  X-component  integral 
scale  length  for  cross-flow  injection  at  j-  =  15.75 
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Figure  5.20.  Effect  of  injection  pressure  and  orifice  diameter  on  X-component  integral 
scale  length  for  counter-flow  injection  at  |  =  15.75 
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Figure  5.21. 


Fifed,  of  injection  pressure  and  orifice  diameter  on  Y-component  integral 
scale  length  for  co-llow  injection  at  £  =  15.75 


Figure  5.22.  Effect  or  injection  pressure  and  orifice  diameter  on  Y-component  integral 
scale  length  for  cross-flow  injection  at  |  =  15.75 
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Figure  S.23. 


Effect  of  injection  pressure  ami  orifice  diameter  on  V-component  integral 
scale  length  for  counter-flow  injection  at  £  =  15.75 
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Figure  5.24.  Effect  of  injection  pressure  and  orifice  diameter  on  X-coinponent  integral 
scale  length  for  co-flow  injection  at  f  =  25.07 
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Figure  5.25.  lifted  of  injection  pressure  ami  orifice  diameter  on  X-componcnt  integral 
scale  length  for  cross-flow  injection  at  £  =  25.67 
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Figure  5.2G.  EfTect  of  injection  pressure  and  orifice  diameter  on  X-component  integral 
scale  length  for  counter-flow  injection  at  f  =  25.67 
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Figure  5.27. 


Effect  of  injection  pressure  and  orifice  diameter  on  Y-component  integral 
scale  length  for  co-flow  injection  at  £  =  25.07 
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Figure  5.28.  Effect  of  injection  pressure  and  orifice  diameter  on  Y-component  integral 
scale  length  for  cross-flow  injection  at  |  =  25.G7 


5-2G 


s 

E 


< 


SO.O  100.0  150.0  200.0  250.0  300.0 


1.0 

0.9 

0.0 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0.0 


pniT!’TT'Tw,mTi'F|'iTi'ivrt^Tiri‘!'rTn,n,iT''nTn,>i'n 


CounUr-flo*  Injection 

Alt  r  T«*U 

i/b  «  25.67 


OOOOOD  «  0.66  mm 

aaaa*  n  *  |.3z  mm 

MWD  ■  1.08  mm 


50.0  100.0  I  SO.O  200.0  250.0  300.0 


Injection  Pressure  (RPa) 


1.0 


0.9 


o.n 


0.7 


0.6 


0.5 


0.4 


0.3 


0.2 


0.1 


0.0 


Figure  5.29.  lifFcct  of  injection  pressure  ami  orifice  diameter  on  Y-component  integral 
scale  length  for  counterblow  injection  at  f  =  25.07 


elfect  on  the  X-component  integral  scale  length.  This  reinforces  the  observation  that  the 
X-component  turbulence  integral  scale  length  Is  primarily  a  function  of  the  jet-grid  tube 
diameter  and  does  not  dissipate  downstream  of  the  jet-grid. 

Lastly,  from  Figures  5.30-5.32,  the  different  injection  orientations  produced  integral 
scales  that  arc  grouped  together  for  particular  turbulence  levels.  As  an  example,  Figure 
0.31  shows  the  X-component  integral  scales  for  cross-flow  injection  (hatched  symbols),  are 
all  plotted  in  the  same  region  at  both  the  5-G  percent  turbulence  level  and  the  10-12 
turbulence  level.  This  phenomena  may  prove  to  be  useful  for  selecting  a  precise  turbulence 
level  once  the  overall  range  turbulence  intensity  for  all  types  of  jet-grid  injection  has  been 
obtained. 

Figures  5.33-5.35  show  the  Y-component  integral  scale  as  a  function  of  the  Y* 
component  turbulence  intensity  for  each  of  the  three  jet-grid  injection  orifice  diameters. 
These  figures  show  the  Y-component  integral  scale  data  was  more  randomly  distributed 
than  the  X-component  integral  scales  shown  in  Figures  5.30  -  5.32.  For  example,  this 
randomness  of  data  is  readily  apparent  when  one  compares  Figure  5.30  with  Figure  5.33, 
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Figure  5.30. 


Meet  of  X*component  turbulence  intensity  on  the  X-eouiponent  integral 
scale  length  for  D  =  O.GO  nun 
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Figure  5.31. 


EfTect  of  X-coinponent  turbulence  intensity  on  the  X-component  integral 
scale  length  for  D  =  1.32  min 
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Figure  5.32.  Effect  of  X-component  turbulence  intensity  on  the  X-component  integral 
scale  length  for  I.9S  nun 


which  represent  the  X  and  Y  components  of  turhnlence  integral  scale  for  the  same  tests. 
Also,  as  the  injection  orifice  diameter  increased,  the  overall  range  of  values  possible  for  the 
Y-component  turbulence  integral  scale  length  increased.  The  grouping  of  results  obtained 
for  the  Y-component  integral  scale  lengths  from  co-llow,  cross-flow,  and  counter-flow  is  not 
as  readily  apparent  from  these  plots  as  it  was  in  Figures  5.30  -  5.32. 

5.2.2  Turbulence  Microscale  The  turbulence  microscale  size  was  affected  by  jet-grid 
location,  injection  pressure,  and  injection  orifice  diameter.  The  location  of  the  jet-grid  had 
the  most  effect  on  the  microscale  size  and  the  injection  orifice  diameter  had  the  least  effect 
on  the  microscale  size.  Comparing  the  results  tabulated  in  Table  5.1  with  those  results 
of  Table  5.2  for  the  turbulence  microscales,  the  X-component  microscale  was  on  the  order 
of  5-10  fiin  with  jet-grid  injection  at  station  0  (|  =  15.75)  and  d.2-9.9  times  larger  with 
jet-grid  injection  at  station  02  (f  =  25.67)  and  the  Y-coinponcnt  microscale  was  on  the 
order  of  30-60  /un  with  jet-grid  injection  at  |  =  15.75  and  9.*1-10.3  times  larger  with  the 
jet-grid  injection  at  station  02. 

From  Tables  5.1  and  5.2,  the  X-component  turbulence  microscale  is  much  smaller 
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Figure  5.3<1.  Effect  of  Y-component  turbulence  intensity  on  the  Y-componont  integral 
scale  length  for  D  =  1.32  nun 
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Figure  5.35.  Kireci  of  Y-compouent  turbulence  intensity  on  the  Y-component  integral 
scale  length  for  1.9S  mm 


than  the  Y-component  turbulence  microscale,  whereas  the  the  results  of  the  previous 
section  show  the  X-component  integral  scale  length  was  greater  than  the  Y-componont 
integral  scale  length.  This  result  is  expected  since  the  smallest  fluctuations  in  the  free¬ 
st  ream,  as  it  passes  through  the  jet-grid,  arc  in  the  X  direction.  Shearing  of  the  wake 
behind  each  jet-grid  tube  creates  a  high  frequency  disturbance  in  the  free-stream,  in  the 
X-dircdion.  The  high  frequency  disturbances  due  to  shearing  between  the  wake  and  the 
free-stream  produce  small  turbulence  microscales.  The  Y-component  microscale  is  due  to 
Y-dircction  fluctuations  of  the  flow  in  the  wake  of  the  jet-grid  tubes.  The  Y-direction 
fluctuations  are  primarily  dependent  on  the  oscillating  components  of  the  wake,  which  oc¬ 
cur  at  a  much  lower  frequencies  than  the  frequencies  associated  with  shearing  of  the  How 
in  the  X-direction.  As  a  result  of  the  low  frequency  disturbances  to  the  free-stream,  the 
Y-component  turbulence  microscales  arc  relatively  large  compared  to  the  X-component 
turbulence  microscales. 

The  X-component  turbulence  microscale  size  (Au)  and  Y-component  turbulence  mi¬ 
croscale  size  (A„)  are  tabulated  for  each  test  in  Tables  5.1  and  5.2.  The  X-component 
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microscale  size  for  co-flow,  cross-flow,  am!  counter-flow  jet-grid  injection  at  both  jet-grid 
locations,  as  a  function  of  jet-grid  injection  pressures,  is  presented  in  Figures  5.36-5.41. 
The  Y-componcnt  microscale  size  for  co-flow,  cross-flow,  and  counter-flow  jct-gri<l  injec¬ 
tion  at  both  jet-grid  locations,  as  a  function  of  jet-grid  injection  pressures,  is  presented  in 
Figures  5.42-5.47. 

Figure  5.3G  shows  a  pronounced  variation  of  the  X-component  microscale  for  co¬ 
flow  injection  with  the  smallest  diameter  orifice  (0.66  mm)  at  station  0,  especially  at 
the  higher  injection  pressures.  Again,  similar  to  the  related  effect  on  integral  scale  for 
this  test,  this  increase  in  turbulence  microscale  to  over  4  times  its  nominal  value  for  lower 
injection  pressures  was  probably  due  to  supersonic  injection.  Young  and  Han  (Young.  1089) 
experienced  similar  results  in  their  tests  where  supersonic  injection  pressures  were  used  to 
increase  the  microscale  length.  As  seen  in  Figure  5.37  for  cross-flow  injection  and  Figure 
5.38  for  counter-flow  injection,  the  microscale  was  relatively  constant  over  a  75-300  kl'a 
injection  pressure  range. 

At  station  02,  co-flow  injection  shown  in  Figure  5.39  produced  up  to  a  SI  percent 
variation  in  X-component  microscale  size  for  the  1.98  mm  injection  orifice.  Cross-flow 
injection  at  station  02,  Figure  5.40,  and  the  X-component  microscales  for  counter-flow 
injection  at  station  02  are  shown  in  Figure  5.41.  Both  cross-flow  and  counter-flow  injection 
produced  a  variety  of  X-component  microscale  sizes,  depending  on  injection  pressure  and 
orifice  diameter. 

The  Y-componcnt  microscales  for  co-flow  injection  at  station  0  are  shown  in  Figure 
5.42.  The  large  increase  in  Y-component  microscale  size  for  the  0.66  mm  diameter  injection 
orifice  is  attributed  to  supersonic  injection  conditions  that  may  have  been  present.  Y- 
coinponent  microscales  as  a  result  of  cross-flow  injection  at  station  0,  Figure  5.43,  varied 
in  size  depending  on  the  injection  pressure  and  orifice  diameter.  Counter-flow  injection 
at  station  0  increased  the  Y-component  microscale  as  the  injection  pressure  increased,  as 
shown  in  Figure  5.44. 

The  variation  of  the  Y-component  microscale  as  a  function  of  injection  pressure, 
with  the  jet-grid  at  station  02  is  shown  for  co-flow  injection  in  Figure  5.45.  With  the 
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Figure  5.3G. 


Effect  of  Injection  pressure  and  orifice  diameter  ou  X-compononl  microscale 
length  for  co-llow  injection  at  £  =  15.75 
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Figure  5.37.  Effect  of  injection  pressure  ami  orifice  diameter  on  X-component  microscale 
length  for  cross-flow  injection  at  f  =  15.75 
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Figure  5.3  S. 


Filed  of  injection  pressure  ami  orifice  diameter  on  X-component  microscale 
length  for  counter-flow  injection  at  £  =  15.75 
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Figure  5.39.  Effect  of  injection  pressure  and  orifice  diameter  on  X-componont  microscale 
length  for  co-flow  injection  at  j  =  25.67 
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Figure  5.40.  Effect.  of  injection  pressure  and  orifice  diameter  on  X-componeni  microscale 
length  for  cross- flow  injection  at  £  =  25.07 
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Figure  5.41.  Effect  of  injection  pressure  and  orifice  diameter  on  X-component  microscale 
length  for  counter-flow  injection  at  |  =  25.07 


5-35 


SO.O  100.0  150.0  200.0  250.0  300.0 

Injection  Pressure  (KPa) 


100.0 

90.0 

00.0 

70.0 

60.0 

50.0 

40.0 

30.0 

20.0 

10.0 

0.0 


Figure  5.42.  Effect  of  injection  pressure  and  orifice  diameter  on  Y-component  microscale 
length  for  co-flow  injection  at  £  =  15.75 


jet-grid  at  station  02,  no  supersonic  injection  occurred  and  there  was  no  accompanying 
large  increase  in  the  Y-component  turbulence  microscale  as  seen  in  Figures  5.3G  and  5.42. 
Figure  5.4G  indicates  cross-flow  injection  produced  a  relatively  small  change  in  the  Y- 
componenl  microscale  size  for  the  two  smallest  injection  orifices  with  a  small  increase 
in  the  Y-component  microscale  size  for  the  largest  injection  orifice  at  a  150  kPa  injection 
pressure.  For  counter-flow  injection  at  station  02.  Figure  5.47,  the  Y-component  microscale 
length  varied  significantly  depending  on  the  injection  orifice  diameter.  In  this  case,  the  1.32 
mm  injection  orifice  diameter  produced  microscales  twice  the  size  of  the  microscale  sizes 
produced  by  other  injection  orifice  diameters  for  the  same  injection  pressure  and  injection 
direction.  The  reason  for  the  large  increase  of  the  Y-component  microscale  size  for  the 
1.32  mm  orifice  using  counter-flow  injection  is  not  known. 

Cross-plots  of  turbulence  microscale  as  a  function  of  their  respective  component  of 
turbulence,  for  each  injection  orifice  diameter,  are  presented  in  Figures  5.48-5.53.  The 
jet-grid  injection  pressure  and  orientation  are  designated  on  each  figure  accordingly.  As 
seen  in  each  of  these  six  figures,  there  was  a  large  variation  in  the  turbulence  microscale 
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Figure  5  .>13.  FI  Feet  of  injection  pressure  and  orifice  diameter  on  Y-component  micro>ca!«* 
length  for  cross-flow  injection  at  £  =  15.75 
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Figure  5.44.  Effect  of  injection  pressure  and  orifice  diameter  on  Y-component  microscale 
length  for  counter-flow  injection  at  f  =  15.75 
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Figure  5.45.  Fifed  of  injection  pressure  and  orifice  diameter  on  Y*compoitem  microscale 
length  for  codlow  injection  at  ?  =  25.G7 
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Figure  5.*IG.  Effect  of  injection  pressure  and  orifice  diameter  on  Y-componcnl  microscale 
length  for  cross-flow  injection  at  f  =  25.G7 
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Figure  5.‘17. 


Flfec!  of  injection  pressure  sunt  orifice  diameter  on  Y-com  Donenl  micro, >calo 
length  for  counter-flow  injection  at  f  =  25.67 


at  the  low  turbulence  levels  ami  a  small  variation  of  the  turbulence  scales  at  the  high 
turbulence  levels.  (This  result  indicates  that  in  the  study  of  turbulence  microscale  clfcci- 
on  heat  transfer,  the  jet-grid  should  be  placed  far  upstream  of  the  turbine  cascade  in 
order  to  get  the  largest  variation  of  the  microscale  length.)  Also  seen  in  these  figure-, 
both  the  X-componeut  and  Y-compoueni  microscale  lengths  are  inversely  proportional  to 
the  turbulence  intensity.  This  is  in  contrast  to  the  results  shown  in  the  cross-plots  of 
the  turbulence  integral  scale  length  with  their  respective  component  ofiurbulcncc  intensity 
(Figures  5.30  -  5.35).  These  particular  figures  also  serve  to  highlight  the  differences  in 
injection  orientation.  For  example,  from  Figure  5.-IS,  co-flow  injection  produced  nearly 
the  same  size  X-componcnt  microscales  centered  around  the  10  percent  turbulence  level, 
cross-flow  injection  produced  similar  microscale  sizes  in  the  region  of  1 1  percent  turbulence 
level,  and  counterflow  produced  similar  size  microscales  at  12  percent  turbulence  level.  The 
same  trend  is  observed  for  the  X-component  microscale  plots  of  Figure  5.<19  and  Figure 
5.50.  Injection  pressure  seems  to  have  a  more  random  effect  on  the  turbulence  microscale 
size  than  the  turbulence  integral  scale  length  in  thtt  the  microscales  are  not  the  same  size 
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Figure  5.dS. 


Variation  of  tin?  X-componout  turbulence  microscale  with  the  X-compoiiom 
of  turbulence:  O  =  Q.tiG  tnm 


for  a  injection  direction  at  a  particular  level  of  free*. stream  turbulence.  This  was  not  the 
case  for  the  turbulence  integral  scale  lengths  (See  Figures  r».riO  -  5.32). 

At  the  lower  turbulence  level  percent),  counterllow  injection  produced  the  largest 
size  microscales;  Figures  5 .*19  and  5.52  demonstrate  the  effect  of  conmerflow  injection  for 
an  orifice  diameter  of  1.32  mm.  One  possible  explanation  for  the  effect  of  coumcr-llow 
injection  on  the  turbulence  microscale  size  is  that  although  counter-flow  injection  may  have 
a  large  shearing  effect  on  the  free-stream,  any  increased  turbulence  level  that  one  might 
gel  from  an  increased  shearing  due  to  counter-flow  injection  decays  as  the  free-stream 
passes  through  the  jet-grid.  Further  study  in  the  area  of  turbulence  scale  production  and 
dissipation  is  necessary  before  firm  conclusions  may  be  formulated. 


5.3  Turbine  Blade  Beat  Transfer 

5.3.1  Baseline  'Turbine  Blade  Beat  Transfer  Tests  Heat  transfer  on  the  turbine 
blade  was  measured  for  all  test  configurations  given  in  Tables  1.1  and  1.2,  except  test 
numbers  AM,  A15,  A24,  A34,  BM,  B2d,  CM,  C2*l,  and  C34.  The  results  of  these  tests 
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Figure  5.51.  Variation  of  the  Y*component  turbulence  microscale  with  the  Y-component 
of  turbulence:  D  =  0.06  mm 
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wore  confmetl  to  turbulence  measurements  only.  The  results  are  presented  in  the  form  of 
Stanton  number  (St)  and  Mussel l  number  (;Vn)  on  the  turbine  blade.  The  Stanton  number 
and  Nusscll  number  on  the  turbine  blade  without  the  jet-grid  in  the  free* stream  flow  are 
designated  by  subscript  x.  With  the  jet-grid  at  station  0  (£  =  15.75),  the  Stanton  number 
and  Nusselt  number  are  designated  with  a  subsript  0  Lastly,  with  the  jet-grid  at  station 
02  (f  s  25.07).  the  Stanton  number  and  Nussclt  number  are  designated  by  subscript  02. 

The  Stanton  and  Nussclt  numbers  on  the  turbine  blade  with  no  jet-grid  in  the  free- 
stream  How  arc  presented  in  Figures  5.5*1  and  5.55,  respectively.  The  Stanton  and  Nusselt 
numbers  with  the  jet-grid  at  station  0  are  presented  in  Figures  5.5G  and  5.57,  respectively. 
The  Stanton  and  Nussclt  numbers  with  the  jet-grid  at  station  02  arc  presented  i:i  Figures 
5.5S  and  5.59.  These  three  heat  transfer  tests  form  a  baseline  for  comparison  with  all.  other 
heat  transfer  tests. 

The  results  of  all  other  tests  are  presented  as  the  ratio  of  the  Stanton  numbers  on 
the  blade  to  the  Stanton  number  on  the  blade  for  a  baseline  test  configuration.  Nusselt 
numbers  for  all  other  tests  are  presented  as  a  ratio  of  Nusselt  numbers  in  the  same  manner 
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as  {lie  Suwon  number  ratios  were  formed.  This  method  of  presentation  gives  a  graphical 
comparison  for  each  heat  transfer  test  to  the  baseline  heat  transfer  tests.  True  values  for 
Stanton  and  Nussclt  numbers  may  be  calculated  by  multiplying  the  Stanton  number  ratio 
or  Nusselt  number  with  the  respective  baseline  Stanton  or  Nusselt  number. 
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Figure  5.51.  Stanton  number  on  the  turbine  blatle  surface  with  no  jet-grid  installed.  Tim 
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Figure  5.55.  Nusselt  number  on  the  turbine  blade  surface  with  no  jet-grid  installed.  Test 
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Figure  5.57.  Nusselt  number  on  the  turbine  blade  surface  the  jet-grid  at  station  02  (f  = 
15.75).  Test  000 
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Stanton  number  on  the  turbine  blade  surface  with  the  jet-grid  at  station  02 
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Figure  5.59.  Nussolt  number  on  the  turbine  blade  surface  the  jet-grid  at  station  02  (j  = 
25.67).  Test  002 
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All  plots  of  the  Stanton  number  and  Nusselt  number  arc  presented  as  a  function  of 
the  nondimcnsional  distance  (j)  along  the  surface  of  the  turbine  blade.  The  results  from 
the  pressure  surface  arc  given  on  the  left  side  of  each  plot  and  results  front  the  suction 
surface  arc  given  on  the  right  side  of  each  plot.  The  negetive  |  values  have  no  physical 
significance,  but  are  a  result  of  the  plotting  package  used  to  create  the  plots. 

5.3.2  Comjmrison  to  Theoretical  Predictions  for  Heat  Transfer  Figure  5.5-1  shows 
the  Stanton  number  for  the  turbine  blade  surface  without  the  jet-grid  installed.  The  free- 
stream  turbulence  level  was  1.26  percent  and  the  turbulence  microscale  in  the  X-direction 
was  computed  as  CSC  //in.  The  turbulence  integral  scale  length  was  9.-I  mm.  Computation 
of  the  Y-direction  turbulence  microscale  and  integral  scale  failed  to  produce  a  valid  result 
for  the  turbulence  microscale.  This  was  a  direct  result  of  the  turbulent  energy  spectrum 
for  the  Y-direction  ,  which  was  highly  scattcrred.  Recall  that  a  spectrum  that  does  not 
approach  a  -5/3  slope  indicates  a  condition  of  non-uniform  turbulence. 

The  laminar  boundary  layer  heat  transfer  prediction  of  the  Stanton  number 
(Kays,  19S7)  is  included  on  Figure  5.5-1  as  well  as  the  turbulent  turbulent  boundary  layer 
prediction  for  Stanton  number  on  the  turbine  blade  (Ambrok,  1957). 

The  Stanton  number  measured  on  the  pressure  surface  matches  well  with  the  turbu¬ 
lent  boundary  layer  prediction  of  Stanton  number  except  near  the  leading  edge  of  the  tur¬ 
bine  blade  at  about  j  =  -0.2.  In  this  region,  it  has  been  suggested  (Consigny,  1982)  that 
on  the  pressure  suface,  immediately  downstream  of  the  leading  edge,  the  flow  is  strongly 
decelerated  over  a  short  distance  which  results  in  a  laminar  separation  bubble  and  a  lower 
heat  transfer.  The  low  level  heating  rate  was  followed  immediately  by  a  high  heating  rate 
due  to  the  reattaching  turbulent  boundary  layer.  The  surface  pressure  distribution  shown 
in  Figure  5.2  indicates  the  flow  was  highly  stagnated,  hence,  a  laminar  boundary  layer  ex¬ 
ists  in  this  region  .  Figure  5.55  shows  the  Nusselt  number  was  approximately  400  for  this 
region,  hence  heat  transfer  was  low  compared  to  the  theoretical  Nusselt  number  computed 
from  the  turbulent  boundary  layer  theoretical  prediction  of  the  Stanton  number.  The  tur¬ 
bulent  boundary  layer  prediction  assumes  a  fully  turbulent,  attached  boundary  layer  and 
is  unable  to  predict  the  Stanton  number  in  the  region  of  the  pressure  surface  leading  edge. 


In  the  region  of  the  trailing  e<!gc  on  the  pressure  side,  Figure  5.51,  the  measured 
Stanton  number  is  low.  Other  studies  on  turbine  blade  heat  transfer  (Blair,  19$9a). 
(Consigny,  1982)  indicate  the  Stanton  number  measured  in  this  region  should  be  as  high  or 
higher  than  the  turbulent  boundary  layer  prediction  of  the  Stanton  number.  One  possible 
reason  for  the  low  value  of  Stanton  number  in  this  region  is  because  of  inaccuracies  in 
accounting  for  heat  conduction  losses  near  the  trailing  edge.  Th*  proximity  of  the  copper 
bus  bars  and  the  thinness  of  the  blade  section  contribute  to  the  heat  conduction  losses  in 
this  region.  Near  the  trailing  edge,  no  interior  thermocouple  was  placed  directly  under  the 
surface  thermocouples.  The  interior  blade  temperature  in  this  region  was  estimated  using 
a  thermocouple  1  cm  ahead  (towards  the  turbine  blade  leading  edge)  of  this  position  which 
might  have  overestimated  the  actual  conduction  losses.  An  estimate  of  the  heat  conduc¬ 
tion  loss  was  computed  using  the  opposing  suction  surface  thermocouple.  Results  of  this 
calculation  indicate  the  Stanton  number  in  this  area  may  actually  be  as  high  as  0.0015. 
This  being  the  ease,  an  error  bar  is  indicated  as  shown  to  account  for  the  uncertainty 
of  the  conductive  heat  losses  in  this  region.  This  same  argument  applies  to  the  trailing 
edge  region  of  the  suction  surface  where  *  =  1.0.  Elsewhere  on  the  blade,  the  blade  was 
thick  enough  to  allow  interior  thermocouples  to  be  placed,  and  there  were  no  heat  sink 
masses  present,  like  the  bus  bars.  The  results  of  the  remaining  tests  are  presented  as  a 
ratio  of  Stanton  numbers  and  as  a  ratio  of  Nussclt  numbers.  The  Stanton  number  ratio 
and  the  Nussclt  number  ratio  are  valid  over  the  entire  surface  except  in  the  region  of  the 
trailing  edge  of  the  turbine  blade,  for  both  surfaces,  where  a  there  is  larger  uncertainty 
in  the  Stanton  or  Nussclt  number  ratio  due  to  the  uncertainty  associated  with  the  heat 
conduction  losses.  The  uncertainty  of  the  Stanton  number  in  the  region  of  the  trailing 
edge  is  shown  as  an  error  bar  in  Figure  5.54. 

On  the  suction  surface  presented  in  Figure  5.54,  the  measured  Stanton  number  was 
best  correlated  by  the  laminar  boundary  layer  heat  transfer  prediction  for  the  majority  of 
the  surface.  Near  the  trailing  edge  of  the  suction  surface  ( *  «  0.7),  the  laminar  boundary 
layer  transitioned  to  a  turbulent  boundary  layer  and  the  Stanton  number  was  best  pre¬ 
dicted  by  the  Ambrok  (Ambrok,  1957)  turbulent  boundary  layer  prediction  of  the  Stanton 
number.  The  decrease  in  Stanton  number  in  the  region  of  *  =  1.0,  on  the  trailing  edge  of 
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ilio  suction  surface,  is  related  to  the  uncertainty  of  computing  the  heat  conduction  in  this 
region. 

With  the  jet-grid  introduced  at  station  0,  the  frcc-strcam  turbulence  increased  to 
1 1.S5  percent  and  the  turbulence  microscales  decreased  dramatically  for  the  low-turbulence 
case.  This  was  the  second  baseline  case  for  heat  transfer  on  the  turbine  blade.  The  results 
of  this  test  (Test  000)  arc  illustrated  in  Figure  5.5G  and  Figure  5.57.  On  the  pressure 
surface,  the  Stanton  number  increased  above  the  zero  free-stream  turbulence,  turbulent 
boundary  layer  prediction,  except  at  the  leading  edge  where  a  laminar  bubble  separated 
the  How  and  lowered  the  heat  transfer  rate.  On  the  suction  surface,  the  boundary  layer 
on  the  blade  was  initially  laminar,  but  started  transitioning  to  a  turbulent  boundary  layer 
immediately  (j  »  0.1).  With  transition  to  a  turbulent  boundary  layer,  the  Stanton  number 
increased  accordingly. 

Figure  5.57  gives  the  Nussell  number  on  the  turbine  blade.  On  the  suction  surface, 
the  Nussell  number  increases  linearly  towards  the  trailing  edge  in  contrast  to  Figure  5.55. 
where  at  j  =  0.5  there  was  a  leveling  of  the  Nussell  number  followed  by  a  rapid  increase 
of  the  Nussclt  number  at  *  =  0.7. 

The  Stanton  number  for  the  turbine  blade  with  the  jet-grid  at  station  02  (jj  = 
25.07)  is  shown  in  Figure  5.5S.  The  free-stream  turbulence  level  was  G.*17  percent  and 
the  turbulence  microscale  lengths  increased  five  times  over  Test  000  microscale  sizes.  The 
turbulent  boundary  layer  heat  transfer  prediction  underestimated  the  Stanton  number  over 
much  of  the  pressure  surface. 

On  the  suction  surface  shown  in  Figure  5.5S,  the  Stanton  number  increased  from  the 
laminar  boundary  layer  levels  of  Stanton  number  to  Stanton  number  levels  predicted  by 
turbulent  boundary  layer  heat  transfer.  The  transition  from  laminar  to  turbulent  boundary 
layer  seems  to  have  begun  at  |  0.3.  In  this  case,  transition  from  laminar  to  turbulent 

boundary  layer  was  delayed  compared  to  the  case  with  the  jet-grid  at  station  0  where 
transition  started  at  *  =  0.1.  The  delayed  transition  was  due  to  the  lower  free-stream 
turbulence  level  associted  with  the  jet-grid  placed  at  station  02.  Comparing  Figures  5.5*1, 
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5.56,  and  5.58,  the  laminar-to-turbulcnt  boundary  layer  transition  point  advanced  toward 
the  leading  edge  as  the  frec-strcam  turbulence  level  increased. 

The  Nussclt  numbers  for  this  configuration  are  given  in  Figure  5.59.  Similar  values 
of  Nussclt  number  on  the  turbine  blade  surface  were  obtained  with  the  jet-grid  at  station 
0.  This  result  was  expected  since  the  frec-strcam  was  moderately  turbulent  for  both  Test 
000  and  Test  002,  and  no  other  test  conditions  were  changed. 

The  ratio  of  Stanton  numbers  and  Nussclt  numbers  for  the  jet-grid  in  at  station  0 
compared  to  the  no  jet-grid  in  the  flow  is  shown  in  Figures  5.60  and  5.61,  respectively. 
From  Figure  5.60,  the  Stanton  number  increased  on  the  pressure  surface  as  much  as  70 
percent  at  the  leading  edge  over  the  low-turbulcnce  test  (no  jet-grid).  The  largest  effect 
was  on  the  suction  surface  where  the  Stanton  number  increase  by  as  much  as  120  percent 
due  to  the  ratio  of  the  turbulent  boundary  layer  Stanton  number  to  the  laminar  boundary 
layer  Stanton  number.  The  ratio  of  the  Nussclt  numbers  shown  in  Figure  5.61  demonstrate 
the  same  increase  in  heat  transfer  rate  on  the  turbine  blade  surface.  This  large  increase  in 
Stanton  and  Nussclt  numbers  was  due  largely  to  the  increased  frec-strcam  turbulence  level. 
Simonich  and  Bradshaw  concluded  in  their  study  of  turbulence  effects  on  heat  transfer  on  a 
flat  plate,  that  the  Stanton  number  increases  with  free-stream  turbulence  at  a  rate  of  about 
five  percent  for  each  percent  increase  in  turbulence  intensity  (Simonich,  1978).  Based  on 
this  estimate,  a  free-stream  turbulence  level  of  11.85  percent  should  produce  a  56  percent 
increase  in  Stanton  number.  The  predicted  increase  in  Stanton  number  was  oxcoded  at 
several  locations  on  both  surfaces  of  the  turbine  blade. 

From  Figure  5.61,  on  the  pressure  surface,  at  |  =  -0.3,  the  rapid  decrease  of  the 
Nussclt  number  ratio  indiactcs  the  heat  transfer  rate  in  this  region  decreased  from  the  test 
without  the  jet-grid  in  the  free-stream  as  compared  to  the  test  with  the  jet-grid  at  station 
0.  One  possible  expalanation  for  this  rapid  decrease  in  Nussclt  number  is  the  laminar  to 
turbulent  boundary  layer  transition  point  moved  its  position  on  the  turbine  blade  surface 
between  the  two  tests.  The  direction  of  movement  of  this  transition  point  seems  to  be 
towards  the  turbine  blade  leading  edge  since  at  *  =  -0.2  the  Nusselt  number  ratio  in 
Figure  5.61  is  1.48  which  was  significantly  greater  than  the  Nussclt  number  ratios  for  the 
remaining  pressure  surface. 
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Figure  5.G1.  Effect  of  the  jet-grid  at  station  0  (f  =  15.75)  on  the  Nussolt  number 
Nusselt  number  ratio.  Test  000 
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Figure  5.62.  Effect  of  the  jet-grid  at  station  02  (£  =  25.67)  on  the  Stanton  number 
Stanton  number  ratio.  Test  002 

Comparison  of  the  Stanton  and  Nusselt  numbers  for  this  test  to  the  Stanton  and 
Nusselt  numbers  for  the  test  with  no  jet-grid  is  presented  in  Figures  5.62  and  5.03.  Fur 
this  case,  from  Figure  5.62,  the  Stanton  numbers  on  the  pressure  surface  increased  about  20 
percent  over  the  no  jet-gri d  Stanton  number  levels  except  at  |  =  -0.2,  where  the  laminar 
to  turbulent  boundary  layer  transition  point  shifted  towards  the  leading  edge  from  its 
position  during  Test  000,  hence  increasing  the  heat  transfer  rate  at  this  poision  on  the 
turbine  blade  surface.  Figu  e  5.63  confirms  the  shift  in  laminar  to  turbulent  boundary 
layer  transition  towards  the  leading  edge  at  ~  =  -0.3  occurred  with  Test  002.  The  ratio 
of  the  Stanton  numbers,  shown  in  Figure  5.62,  and  the  ratio  of  Nusselt  numbers,  shown  in 
Figure  5.63,  on  the  suction  surface  remained  about  the  same  value  in  this  case  (Test  002) 
as  with  the  jet-grid  located  at  station  0  (Test  000). 

5.3.2. 1  Heat  Tmnfcron  tlic  Turbine  Blade  with  Jet-Grid  Injection  The  effect 
of  jet-grid  injection  on  heat  transfer  is  now  considered.  In  all  figures  pertaining  to  this 
section,  multiple  test  results  are  shown  on  a  single  figure.  On  each  figure,  the  injection 
orifice  diameter  and  injection  direction  are  identical  for  all  tests;  however,  the  injection 
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Figure  Effect  of  the  jet-grid  at  station  0  (£  =  25.G7)  on  the  Nusscll  number 
Nusscli  number  ratio.  Test  002 

pressure  was  varied.  The  Stanton  and  Nusscll  numbers  for  the  lowest  jet-grid  injection 
pressure  is  designated  by  square  symbols.  The  results  from  next  highest  injection  pressure 
are  designated  with  triangle  symbols,  and  the  results  from  the  highest  injection  pressure 
are  designated  with  diamond  symbols.  Actual  injection  pressures  for  each  test  are  given 
in  Tables  1.1  and  1.2. 

Figures  shown  in  this  section  gi  the  Stanton  number  ratios  and  the  Nusscll  number 
ratios  measured  on  the  turbine  biade  surface  with  the  jet-grid  at  station  0  (|  =  15.7,’i). 
Only  representative  plots  of  heal  transfer  measurements  in  the  form  of  Stanton  and  Nusscll 
number  ratios  arc  given  in  this  section.  All  other  plots  of  the  Stanton  and  Nussclt  number 
ratios  are  shown  in  Appendix  D. 

First,  co-flow  with  jet-grid  injection  shown  in  Figure  5.G1  had  the  least  effect  on 
the  Stanton  number  on  the  turbine  blade  pressure  surface.  As  may  be  seen  in  Figures 
5.G5,  5.G7,  and  5.G9,  counter-flow  injection  increased  the  Stanton  number  on  the  pressure 
surface  up  to  five  times  the  nominal  value  (with  no  jet-grid  injection).  This  would  suggest 
an  enormous  heat  transfer  rate,  but  the  corresponding  figures  for  Nussclt  number  ratios  for 
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Figure  5.G-I.  Stanton  number  on  the  turbine  blade  for  co-flow  injection.  D  =  0.0G  mm. 
Test  All,  A12  and  A13.  Jet-grid  at  station  0 


counter-flow  injection  ,  Figures  5.GG,  5.6$,  and  5.70,  show  the  Nusselt  number  ratios  are 
unity  or  less.  There  seems  to  be  a  dichotomy  of  results;  however  this  is  not  the  case,  rather 
the  turbine  blade  surface  velocity  distribution  changed  as  jet-grid  injection  was  introduced 
which  changed  the  Stanton  number.  The  Stanton  number  is  inversely  proportional  to 
the  local  Reynolds  number,  lie,,  on  the  blade  surface.  Counter-flow  injection  changed  the 
velocity  distribution  on  the  turbine  blade  pressure  surface  by  stagnating  the  How.  Pressure 
distributions  on  the  turbine  blade  surface  measured  for  counter-flow  injection  support  this 
conclusion  (sec  Figures  5.5  and  5.8).  Since  the  Reynolds  number  is  proportional  to  the 
local  velocity,  any  change  in  velocity  would  result  in  a  change  in  Reynold  number  which 
gives  an  equally  large  change  in  Stanton  number.  The  heat  transfer  coefficient  did  not 
change  and  the  Nusselt  number  for  this  configuration  confirms  this  fact.  Rather  the  local 
velocity  on  the  turbine  blade  surface  changed,  which  affected  the  Stanton  numbers. 

The  dependence  of  the  Stanton  number  on  the  turbulence  microscale  is  more  difficult 
to  determine  because  of  the  Stanton  number’s  inherent  dependence  on  the  local  velocity 
distribution,  however  from  Figures  5.71  and  5.72,  where  the  free-stream  turbulence  was 
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Figure  5.05.  Stanton  number  on  the  turbine  blade  for  counter-flow  injection.  I )  =  O.UtJ 
mm.  Test  Cl  I .  CM2  and  CM3.  Jet-grid  at  station  0 
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Figure  5.G6.  Nussclt  number  on  the  turbine  blade  for  counter-flow  injection.  J)  =  U.tiii 
mm.  Test  Cl  1 ,  CM2  and  CM3.  Jet-grid  at  station  0 
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Figure  5.70.  Nussclt  number  on  the  turbine  blade  for  counter-flow  injection.  D  =  1.9$ 
nun.  Test  C31,  C32  and  C33.  Jct-grid  at  station  0 
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Figure  5.71.  Stanton  number  on  the  turbine  blade  for  cross* flow  Injection.  I)  =  O.CG  mm. 
Test  HU.  M2  ami  M3.  Jet-grid  at  station  0 


matched  for  all  three  tests  shown  in  this  figure,  and  where  the  turbine  blade  surface  pres*  im- 
distribution  from  Figure  5.1  was  constant  for  all  three  tests,  the  increase  in  Stanton  number 
ratio  is  attributable  to  the  respective  decrease  in  X-component  microscale.  Note,  due  to 
room  restrictions  on  the  plot,  only  the  X-component  microscales  have  been  listed,  however 
the  corresponding  Y-componcut  microscales  for  the  X-component  microscales  listed  in 
Figures  5.71  and  5.72  similarly  decreased  in  size  for  these  three  tests. 

The  effect  of  turbulence  microscale  on  the  turbine  blade  Nusselt  number  is  seen  when 
the  free-stream  turbulence  was  matched  for  different  tests.  Cross-flow  injection  pre  ..iced 
measureable  variations  in  the  turbulence  microscale  for  the  same  turbulence  level.  Figures 
5.73,  5.74,  and  5.75  indicate  the  Nusselt  number  increased  on  the  pressure  surface  as  the 
turbulence  microscale  length  decreased  for  the  same  turbulence  level.  The  suction  surface 
Nusselt  number  ratios  were  unaffected  by  the  variation  in  the  microscale  size,  which  was 
expected  since  the  microscale,  which  is  a  measure  of  turbulence  intensity  and  only  affects 
the  turbulent  boundary  layer,  but  the  suction  surface  boundary  layer  was  laminar  or  in 
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Figure  5.72.  Stanton  number  on  the  turbine  blade  for  cross* flow  injection.  I)  =  1.32  mm. 
Test  1)21 , 1122  and  1)23.  Jet-grid  at  station  0 


transition  from  a  laminar  to  a  turbulent  boundary  layer  and  was;  therefore,  unaffected  b\ 
the  variation  in  the  turbulence  microscale  size. 

Figure  5.7(1  shows  the  Mussel i  number  for  both  the  pressure  surface  and  the  surtjun 
surface  was  different  for  all  three  tests  in  this  figure,  but  the  free- si  ream  turbulence  w«* 
only  matched  for  the  two  lower  injection  pressures  shown  in  this  figure.  For  the  two  tests 
of  Figure  5.7(1  where  the  turbulence  level  was  matched,  only  a  small  increase  in  the  Nussell 
number  ratio  is  noticeable  near  the  leading  edge  of  the  pressure  surface.  On  the  suction 
surface,  there  is  no  discernablc  difference  between  the  two  test  results.  The  noticeable 
increase  of  the  Nussclt  number  ratio  for  the  test  with  the  highest  injection  pressure  may  or 
may  not  be  attributible  to  the  decrease  in  microscale  length  since  the  free* stream  turbulence 
level  was  roughly  1  percent  higher  than  the  free-stream  turbulence  level  of  the  other  two 
tests. 

Representative  plots  of  the  Stanton  number  and  Nussclt  number  ratios  for  the  turbine 
blade  with  the  jot-grid  at  station  02  (j  =  25.07)  are  now  presented. 

In  general,  the  Stanton  numbers  ratios  measured  on  the  turbine  blade  were  ui  the 
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Figure  5.76.  Nusselt  number  on  the  turbine  blade  for  cc-flow  injection.  D 
Tesi  A31,  A32  and  A33.  Jet-grid  at  station  0 
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Figure  5."7.  Nusselt  number  on  the  turbine  blade  for  co-flow  injection.  I)  -  1.9s  mm. 
Test  1)31, 1)32  and  1)33.  )et*grid  at  station  02 


order  of  unity  or  slightly  less  than  unity  over  the  pressure  and  suction  surfaces  of  tin* 
turbine  blade.  The  Nusselt  number  ratios  were  also  on  the  order  of  unity,  except  in  the 
case  of  co-flow  injection  through  the  1.9S  mm  diameter  orifice  shown  in  Figure  5.77.  In 
this  plot,  the  Nusselt  number  varied  widely,  and  especially  along  the  presure  surface.  Tin- 
effect  cf  the  turbulence  microscale  on  the  Nusselt  number  is  not.  clear  in  this  plot  since  tin 
turbulence  intensity  was  not  held  constant  for  all  tests  shown  in  this  figure. 

Again,  the  ratios  refer  to  Test  002,  not  Test  XXX.  Although  the  ratio  may  be  unity 
or  less,  the  actual  Stanton  or  Nusselt  numbers  arc  generally  much  high.jr  than  the  Stanton 
and  Nusselt  numbers  on  the  turbine  biade  without  the  jet-grid  installed  in  the  frec-stream 
flow. 

The  most  noticeable  trend  for  this  set  of  figures  is  that  no  large  increase  in  the 
Stanton  number  occurred  for  counter-flow  injection,  Figures  5.78,  5.79,  and  5.80,  as  it  did 
•it  station  0  (Compare  to  Figures  5.05,  5.67,  and  5.69).  The  local  surface  velocities  on 
the  pressure  surface  of  the  turbine  blade  remained  relatively  unafTcctcd  which  is  visible 
from  Figure  5.8.  But  counter-flow  injection  did  cause  the  largest  variation  of  the  Stanton 
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number  as  compared  to  the  other  injection  directions  tested  with  the  jet-grid  at  station 
02.  For  example,  compare  Figure  5.79  with  Figure  5.S1. 

For  all  heat  transfer  tests  made  at  station  02,  the  free-stream  turbulence  and  surface 
velocity  distribution  were  matched  only  in  a  two  cases.  Figure  5.81  and  Figure  5.S2  show 
a  small  increase  in  the  Stanton  number  ratio  and  the  Nusselt  number  ratio,  respect ivoly, 
of  each  test  as  the  turbulence  microscale  size  decreased.  It  is  not  known  why  the  increase 
in  the  Stanton  number  ratios  and  the  Nusselt  number  ratios  was  so  small  for  a  27  percent 
decrease  in  the  tubulence  X-component  microscale  size. 

Figure  5.80  represents  the  second  case  where  both  the  free-stream  turbulence  and 
surface  velocities  on  the  pressure  surface  of  the  turbine  blade  were  matched  for  the  two 
tests  shown.  In  this  case,  the  Stanton  number  ratio  increased  as  the  turbulence  microscale 
decreased,  the  Nusselt  number  ratio  shown  for  this  case  in  Figure  5.83  remained  relatively 
constant  for  both  tests. 
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Figure  5.79.  Stanton  number  on  the  turbine  blade  for  conn  ter*  flow  injection.  I)  =  1.22 
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Figure  5.80.  Stanton  number  on  the  turbine  blade  for  counter-flow  injection.  D  =  1.9$ 
mm.  Test  F31,  F32  and  F33.  Jet-grid  at  station  02 
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Figure  5.82.  Nussclt  number  on  the  turbine  blade  for  co-flow  injection.  D  =  1.32  nun. 
Test  D21, 1)22  and  D23.  Jet-grid  at  station  02 
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17.  CONCLUSIONS  AND  RECOMMENDATIONS 


ti.l  Conclusions 

The  dice!  of  turbulence  scale  lengths  on  heat  transfer  in  a  linear  turbine  cascade  was 
studied.  Of  equal  importance,  methods  were  developed  to  control  the  turbulence  integral 
scale  anti  microscale  lengths  at  a  given  free- stream  turbulence  level.  The  overall  results  of 
this  experiment  arc  presented  as  turbulence  scale  generation  and  turbulence  scale  dfects 
on  the  Stanton  and  Nusselt  numbers  for  the  turbine  blade. 

The  following  conclusions  are  given  for  the  turbulence  scale  results: 

1.  The  turbulence  integral  scale  length  (A)  was  primarily  dependent  on  the  physical 
dimensions  of  the  turbulence  generating  grid  tube.  The  sireamwisc  component  of 
integral  scale  were  on  the  order  on  0.6  mm,  which  corresponds  closely  with  tin' 
outside  diameter  of  the  jet-grid  tube  inserted  into  the  free-stream.  In  this  studv. 
secondary  How  injection  through  a  jet-grid  placed  in  the  free-stream  increased  the 
integral  scale  length  by  as  much  as  *15  percent.  Co-How  injection  increased  and 
decreased  the  integral  scale  lengths  depending  on  injection  pressure  and  injection 
orifice  diameter.  Cross-How  injection  has  negligible  effect  on  the  size  of  the  integral 
scales.  Counter-flow  injection  increased  the  effective  jet-grid  tube  diameter. 

2.  Turbulence  integral  scales  did  not  decay  downstream  of  the  turbulence  generating 
grid. 

3.  The  turbulence  microscale  length  (A)  was  primarily  determined  by  the  location  of 
the  jet-grid  in  the  free-stream  with  respect  to  the  position  of  the  turbulence  measure¬ 
ment  device.  A  large  variation  of  the  microscales  occurred  for  the  jet-grid  location 
when  the  injection  direction  was  changed.  Co-flow  injection  through  the  jet-grid 
tended  to  decrease  the  microscale  lengths.  Cross-flow  injection  produced  a  variety  of 
microscale  lengths.  Counter-How  injection  had  a  relatively  minor  afTcct  on  the  mi¬ 
croscale  lengths.  The  largest  variation  of  the  turbulence  microscales  due  to  injection 
occurred  with  the  jet-grid  placed  at  the  far  upstream  location  =  25.67). 
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•1.  Turbulence  microscales  grew  in  size  as  the  free- stream  flow  moved  downstream  of  the 
turbulence  generating  device.  Turbulence  microscales  measured  at  a  noudimoiiMonnl 
distance  of  §  =  25.67  were  5-S  times  the  size  of  turbulence  microscales  measured  at 
f  =  15.75. 

5.  Turbulence  intensity  varied  inversely  with  turbulence  microscale  size.  The  higher 
free-stream  turblcncc  levels  contained  the  smallest  microscale  sizes.  Turbulence  de¬ 
cay  matched  turbulence  decay  correlations  provided  by  Blair  (Ulair,  19S3a)  for  the 
jet-grid  at  |  =  15.75.  Blair’s  correlation  underestimated  the  turbulence  decay  by  16 
percent  with  the  jet-grid  at  |  =  25.67.  In  this  case,  the  effect  on  the  turbulence  in¬ 
tensity  of  the  free-stream  accelerating  as  it  approached  the  turbine  cascade  is  thought 
to  have  caused  the  measured  free-stream  turbulence  level  to  decrease  beyond  Blair's 
correlation. 

Conclusions  on  the  heat  transfer  results  are  listed  as  follows: 

1.  The  Stanton  number  measured  on  the  turbine  blade  pressure  surface  matches  the  tur¬ 
bulent  boundary  layer  heat  transfer  prediction  presented  by  Ambrok  (Ambrok,  1957) 
for  low-turbulence  ( 1 .26  percent )  free-stream  conditions.  Near  the  leading  edge  of  the 
pressure  surface,  a  laminar  bubble  is  thought  to  have  caused  a  signiglicantly  lower 
Stanton  number  than  what  is  predicted  by  Ambrok’s  theory. 

2.  The  Stanton  number  measured  on  the  suction  surface  is  best  correlated  by  laminar 
boundary  heat  transfer  predictions  (Kays,  19S7)  for  low-turbulence  free-stream  flow. 
An  abrubt  increase  in  Stanton  number  on  the  suction  surface  occurod  around  the 
non-dimensional  surface  distance  of  |  =  0.6,  at  which  point,  the  Stanton  number 
matched  the  turbulent  boundary  layer  prediction.  It  is  thought  laminar  to  turbulent 
boundary  layer  transition  occured  at  this  location  on  the  suction  surface. 

3.  Introduction  of  the  jet-grid  device  produced  significantly  higher  levels  of  free-stream 
turbulence,  which  caused  a  higher  heat  transfer  rate  and  a  higher  Stanton  number 
than  predicted  by  Ambrok’s  theory  on  the  pressure  surface.  Higher  free-stream 
turbulence  advanced  the  laininar-to-turbulent  boundary  layer  transition  point  on 
the  suction  surface.  For  a  11.8  percent  free-stream  turbulence  level,  the  boundary 
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layer  on  the  suction  surface  was  almost  entirely  in  a  transition  state  or  fully  turbulent 
ami  the  Stanton  number  approached  turbulent  boundary  layer  predictions. 

*1.  With  the  free-stream  turbulence  level  held  constant  and  for  a  constant  turbine  blade 
surface  velocity,  the  Stanton  number  increased  as  the  turbulence  microscale  length 
decreased  in  repeated  heat  transfer  tests.  However ,  the  jet-grid  introduced  a  vari¬ 
ation  of  the  local  velocity  along  the  turbine  blade  surface,  which  had  a  large  effect 
on  the  measured  Stanton  number  and  made  it  difficult  to  correlate  the  turbulence 
microscale  lengths  with  the  Stanton  number  for  many  of  the  heat  transfer  tests  per¬ 
formed  in  this  study. 

5.  The  effect  of  the  turbulence  microscale  length  on  the  Nusselt  number  was  best  ob¬ 
served  with  cross-flow  injection  with  the  jet-grid  at  station  0,  where  it  was  observed 
that  the  free-stream  turbulence  level  remained  constant,  the  turbine  blade  surface 
velocities  remained  constant,  the  turbulence  microscales  decreased  with  increasing 
injection  pressure,  and  the  Nusselt  number  increased. 

G.2  Ihconunciulntions 

In  light  of  the  results  of  this  report,  the  following  recommendations  are  made  fur 
advanced  study  of  turbulence  scale  effects  on  heat  transfer  in  a  linear  turbine  cascade: 

Instrumentation  improvements: 

1.  The  turbine  blade  heat  transfer  test  model  and  pressure  surface  model  should  be 
integrated  as  one  test  model.  This  would  allow  the  experimenter  to  simultaneously 
measure  pressure  distributions  and  heat  transfer  on  the  turbine  blade  surface. 

2.  An  additional  Hewlett-Packard  -M713,  2-1-Channel  FET  Multiplexer  should  be  in¬ 
corporated  with  the  present  data  acquisition  system.  This  improvement  would  allow 
more  thermocouples  to  be  used  with  the  test  model  and  allow  a  more  detailed  mea¬ 
surement  of  the  turbine  blade  surface  temperatures. 

3.  Improved  the  jet-grid  device  by  attaching  high-pressure  air  lines,  rather  than  the  low 
pressure  air  lines  presently  used.  This  would  allow  higher  injection  pressures  to  be 
used  and  supersonic  injection  effects  may  be  studied. 


6-3 


•t.  The  test  section  eudwall  should  he  modified  to  allow  easier  access  to  the  test  section 
for  introduction  and  removal  of  the  test  models. 

Areas  of  further  research; 

1.  Investigate  the  free- stream  turbulence  decay  and  microscale  decay  through  the  cas¬ 
cade  passageway. 

2.  Measure  the  effect  of  free- stream  turbulence  level  and  turbulence  microscale  length 
on  the  turbine  blade  boundary  layer. 

It.  Model  turbine  blade  surface  degradation  with  a  rough  surface  model  and  study  the 
effects  on  surface  heat  transfer  and  the  location  of  transition  to  from  laminar  to 
turbulent  boundary  layers  on  both  the  pressure  surface  and  the  suction  surface. 
Rildets  on  the  blade  surface  may  also  be  incorporated  into  the  study. 

>1.  Lastly,  perform  a  flow  visualization  study  in  order  to  determine  the  three-dimensional 
effects  of  the  low-aspect  ratio  test  section.  Use  the  flow  visualization  to  study  the 
effects  of  jet-grid  injection  on  the  local  velocities  on  the  surface  of  the  turbine  blade. 
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Appendix  A.  UlCAT  TR  A  ESTER  DATA  FOR  IUSEUSE  BLADE 

comaviiATioEs 


Heal  transfer  data  for  the  turbine  cascade  with  no  jet-grid  installed:  Test  XXX 


File:  xxx 

Secondary  flow  pressure  = 
Eleotrical  energy  input  = 
y./o 


1.73910E-02 
6 . 52170E-02 
1 . 24349E-01 
1 . 91304E-01 
3 . 15652E-01 
4 . 36522E-Q1 
6 . 00000E-01 
7.68696E-01 
8. 52174E-01 
9 . 23478E-01 
9 . 91304E-01 
9 . 02609E-01 
7.64348E-01 
6 . 66957E-01 
5 . 81739E-01 
4 . 869S7E-01 
3.85217E-01 
3 . 32174E-01 
2 . 19130E-01 
1 . 07826E-01 
3 . 56S20E-02 
0 . OOOOOE+OO 
8 . 69565E-03 


Re  no. 

0. OOOOOE+OO 
1 . 33630E+03 
2.77358E+04 
3.74176E+04 
5.73666E+04 
1.03032E+05 
1.96877E+05 
3.21277E+0S 
4.66597E+05 
5.87590E+05 
7 . 00919E+05 
7 . 17577E+05 
5.77157E+05 
5.59S42E+0S 
5.07137E+05 
4.42468E+05 
4 . 18384E+05 
3 . 51226E+0S 
2 . 68639E+05 
1 . 81815E+05 
7 . 98340E+04 
4 . 24161E+04 
2 . 23173E+04 


0  KPa 
325  Watts 
Husselt 
0. OOOOOE+OO 
0.  OOOOOE+OO 
4.03821E+02 
8.13716E+02 
8.95459E+02 
1 . 05932E+03 
1.42988E+03 
1.29002E+03 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
1 . 19522E+03 
2 . 31694E+03 
2 . 04288E+03 
9 . 04001E+02 
8.40041E+02 
8 . 73489E+02 
7 . 11 123E+02 
5.44172E+02 
2 . 99230E+02 
2. 17413E+02 
1 . 22888E+02 


Stanton 
0. OOOOOE+OO 
0. OOOOOE+OO 
2 . 05935E-02 
3. 07594E-02 
2 . 20784E-02 
1.45424E-02 
1 . 02727E-02 
5. 67935E-03 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
2.92909E-03 
5. 85682E-03 
5 . 69768E-03 
2.88980E-03 
2 . 83992E-03 
3 . 51764E-03 
3.74417E-03 
4 . 23338E-03 
5 . 30149E-Q3 
7 . 24998E-03 
7 . 78844E-03 


A-l 


Heat » 


for  the  turbine  cascade  jet -grid  installed  .  no  blowing .  £  =  15.“ 


Test  000 


File:  000 

Secondary  flow  pressure  = 
Electrical  energy  input  = 
x/c  Re  no. 


1 . 73910E-02 
6 . 52170E-02 
1 . 24343E-01 
1 . 91304E-01 
3 . 15652E-01 
4 . 36S22E-01 
6 . 00000E-01 
7 . 68696E-01 
8. 52174E-01 
9.23478E-01 
9 . 913Q4E-01 
9.02609E-01 
7 . 64348E-01 
6. 66957E-01 
5 .  81739E-01 
4 . 86957E-01 
3.85217E-01 
3. 32174E-01 
2 . 19130E-01 
1 . 07826E-01 
3. 56520E-02 
0 . 00000E+00 
8 . 69565E-03 


0 . OOOOOE+OO 
5 . 28691E+03 
2 . 28461E+Q4 
2 . 3495QE+04 
3.70269E+04 
8.12465E+04 
1 . 71917E+05 
2 . 94917E+05 
4 . 28749E+05 
S.48295E+05 
6. 59944E+05 
6 . 70384E+05 
5 . 70136E+05 
5 . 24245E+0S 
4 . 76330E+05 
4 . 27040E+05 
4 . 06611E+05 
3 . 31049E+05 
2 . 51700E+05 
1 . 67939E+05 
7 . 24332E+04 
3 . 72762E+04 
1.93108E+04 


0  KPa 
325  Watts 
Nusselt 
0. 00000E+00 
0 . Q0000E+00 
5 . 83206E+02 
8.00154E+02 
9 . 89404E+02 
1 . 25927E+03 
1.70049E+03 
1 . 57382E+03 
0. 00000E+00 
0. OOOOOE+OO 
0. 00000E+00 
0. OOOOOE+OO 
2 . 36502E+03 
2.41284E+03 
2 . 20921E+03 
1 . 96938E+03 
1 . 70847E+03 
1 . 30467E+03 
1 . 04757E+03 
7 . 90401E+02 
3 . 99754E+02 
2 . 55815E+02 
1 .4 1271E+02 


Stanton 
0. OOOOOE+OO 
0. OOOOOE+OO 
3 . 61070E-02 
4 . 81702E-02 
3 . 779S2E-02 
2.19226E-02 
1 . 39907E-02 
7 . 54809E-03 
0. OOOOOE+OO 
0. OOOOOE+OO 
0.  OOOOOE+OO 
0.  OOOOOE+OO 
5 . 86728E-03 
6 . 50990E-03 
6 . 56009E-03 
6 . 52292E-03 
5 . 93609E-03 
5 . 57428E-03 
5 . 88682E-03 
6 . 65696E-03 
7 . 80613E-03 
9 . 70678E-03 
1 . 03474E-Q2 


Heat  transfer  data  for  the  turbine  cascade  jet-grid  installed, 
Test  002 


File:  002 

Secondary  flow  pressure  = 
Electrical  energy  input  = 
x/c  Re  no. 


1 . 73910E-02 
6.52170E-02 
1 . 24349E-01 
1. 91304E-01 
3 . 1S652E-01 
4 . 38522E-01 
6 . 00000E-01 
7 . 68696E-01 
8 . 52174E-01 
9. 23478E-01 
9 . 91304E-01 
9. 02609E-01 
7 . 64348E-01 
6. 66957E-01 
5.81739E-01 
4 . 86957E-01 
3 . 85217E-01 
3.32174E-01 
2.19130E-01 
1 . 07826E-01 
3 . 58520E-02 
0 . 00000E+00 
8 . 69565E-03 


O.OOOOOE+OO 
5 . 06586E+03 
2.43590E+04 
3 . 01694E+04 
4 . 80176E+04 
9 . 12736E+04 
1 . 80266E+05 
3 . 01917E+05 
4 . 36054E+05 
5 . 51019E+05 
6 . 59929E+05 
6 . 75457E+05 
5 . 78900E+05 
5 . 24732E+05 
4 . 79838E+05 
4 . 31781E+05 
4 . 09825E+05 
3 . 34915E+05 
2.53604E+05 
1 . 69121E+05 
7 . 35434E+04 
3 . 85336E+04 
2 . 01419E+04 


0  KPa 
325  Watts 
Husselt 
0.00000E+00 
0.00000E+00 
5 . 30052E+02 
7 . 53089E+02 
9.46641E+02 
1. 20913E+03 
1.61313E+03 
1 . 52730E+03 
0.00000E+00 
O.OOOOOE+OO 
O.OOOOOE+OO 
O.OOOOOE+OO 
2 . 62639E+03 
2 . 33533E+03 
2 . 25002E+03 
1 . 98047E+03 
1 . 60245E+03 
1 . 27357E+03 
9 . 64510E+02 
6 . 79014E+02 
3.43152E+02 
2 . 54121E+02 
1.39989E+02 


no  blowing,  f  =  25.0 


Stanton 
O.OOOOOE+OO 
O.OOOOOE+OO 
3.01589E-02 
3 . 53070E-02 
2.78847E-02 
1 . 87373E-02 
1 . 26571E-02 
7.15513E-03 
O.OOOOOE+OO 
O.OOOOOE+OO 
O.OOOOOE+OO 
O.OOOOOE+OO 
6.41705E-03 
6.29494E-03 
6 . 63242E-03 
6.48762E-03 
5 . 53053E-03 
5 . 37859E-03 
5 . 37937E-03 
5 . 67888E-03 
6 . 59969E-03 
9 . 32784E-03 
9 . 83048E-03 


Appendix  B.  TURBULENCE  ENERGY 


This  appendix  contains  plots  of  the  turbulence  energy  spectrum  used  to  calculate 
the  turbulence  integral  scale  lengths  and  turbulence  microscale  lengths  for  each  test.  The 
theoretical  slope  of  the  spectrum  (Batchelor,  1953)  is  shown  for  comparison. 
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Figure  13.30.  Turbulence  energy  spectrum  for  v*:  Test  1321 
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Figure  13.37.  Turbulence  energy  spectrum  for  i/‘:  Test  1)22 
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Figure  13.38.  Turbulence  energy  spectrum  for  v‘:  Test  1322 
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Figure  1M9.  Turbulence  energy  spectrum  for  «  *:  Test  IV23 
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Figure  Turbulence  energy  spectrum  for  «  *:  Test  1121 
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Figure  B.<I2.  Turbulence  energy  spectrum  for  v~:  Test  B2*l 
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Figure  H.<I7.  Turbulence  energy  spectrum  for  u~:  Test  1)33 
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Figure  U.‘1S.  '1'urbulcnce  energy  spectrum  for  v  :  'lest  B33 
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Figure  Turbulence  energy  spectrum  for  i/  *:  Test  1)3-1 
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Figure  B.50.  Turbulence  energy  spectrum  for  i/2;  Test  1)3-1 
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Figure  11.51.  Turbulence  energy  spectrum  for  u  *:  Test  C’ll 
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Figure  B.52.  Turbulence  energy  spectrum  for  v'2:  Test  Cl  I 
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Figure  B.55.  Turbulence  energy  spectrum  for  u  :  Test  CM3 
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Figure  B.57.  Turbulence  energy  spectrum  Tor  u‘:  Test  C‘M 
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Figure  B.oS.  Turbulence  energy  spectrum  for  v  :  Test  CM 
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Figure  B.G1.  Turbulence  energy  spectrum  for  t /*: 
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Figure  B.G2.  Turbulence  energy  spectrum  for  v~:  Test  ('22 
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Figure  R.G-1.  Turbulence  energy  spectrum  for  e 


1 2 


Test  C23 


•  (n.)  (m"/  s") 


io- 

io- 

io- 

10  '* 

io  ** 

io  " 


10*  io’ 

Frequency  (Hz) 


to  • 


Figure  U.G5.  Turbulence  energy  spectrum  for «  *:  Test  C2*l 
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Figure  Il.GG.  Turbulence  energy  spectrum  for  v  *:  Test  C‘2'1 
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Figure  H.C7.  Turbulence  energy  spectrum  for  it  :  Test  031 
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Figure  13.GS.  Turl)ulcnce  energy  spectrum  for  r  *:  Test  C31 
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Figure  H.71.  Turbulence  energy  spectrum  for  «'*:  Test  C'33 
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Figure  B.72.  Turbulence  energy  spectrum  for  v'7:  Test  C33 
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Figure  B.7G.  Turbulence  energy  spectrum  for  i>  :  Test  1)11 
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Figure  H.77.  Turbulence  tuor^v  spectrum  for  »'*:  Test  1)12 
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Figure  B.7S.  Turbulence  energy  spectrum  for  v  *:  Test  D12 
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Figure  B.X3.  Turbulence  energy  spectrum  for  u  :  Test  1)22 
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Figure  H.S5.  Turbulence  energy  spectrum  for  i/*:  Test  1)23 
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Figure  B.8G.  Turbulence  energy  spectrum  for  v  *:  Test  D23 
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Figure  IJ.S7.  Turbulence  energy  spectrum  for  i/‘:  Test  l):tl 
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Figure  H.S5  Turbulence  energy  spectrum  for  t/2:  'lest  1)32 
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Figure  B.90.  Turbulence  energy  spectrum  for  t/2;  Test  D32 
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Figure  B.91 .  Turbulence  energy  spectrum  for  u  :  Test 


Figure  B.92.  Turbulence  energy  spectrum  for  v'2:  Test  D33 
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Figure  B.96.  Turbulence  energy  spectrum  for  u  *:  Test  El2 
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Figure  B.98.  Turbulence  energy  spectrum  for 
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Figure  H.9M.  Turbulence  energy  spectrum  for  u  *:  Test  Iv2l 


Figure  B.100.  Turbulence  energy  spectrum  for  F*:  Test  F'21 
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Figure*  H.101.  Turbulence*  energy  spectrum  for  t/‘:  Test  K‘2'2 
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Figure  H.10'2.  Turbulence  energy  spectrum  for  v'~:  Test  E'22 
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Figure  11.103.  Turbulence  energy  spectrum  for  »/3:  Test  E23 
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Figure  B.IO-J.  Turbulence  energy  spectrum  for  t/2:  Test  E23 
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Figure  I). 108.  Turbulence  energy  spectrum  for  t?  :  Test  E32 
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Figure  B.109.  Turbulence  energy  spectrum  for  u  *:  Test  K33 
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Figure  13.110.  Turbulence  energy  spectrum  for  v  :  Test  E33 
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Figure  B.113.  Turbulence  energy  spectrum  for  u  :  Test  FJ2 


Figure  B.1M.  Turbulence  energy  spectrum  for  t? 
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Figure  11.115.  Turbulence  energy  spectrum  fori/'*:  Test  F13 
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Figure  B.11C.  Turbulence  energy  spectrum  for  t>'2:  Test  F13 
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Figure  H.117.  Turbulence  energy  spectrum  for  u  :  Test  F21 
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Figure  B.118.  Turbulence  energy  spectrum  for  t Test  F21 
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Figure  11.119.  Turbulence  energy  spectrum  for  u  *:  Test  F22 
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Figure  U.120.  Turbulence  energy  spectrum  for  v  :  Test  F22 
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Figure  B.123.  Turbulence  energy  spectrum  for  u  :  Test  F3J 
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Figure  0.12*1.  Turbulence  energy  spectrum  for  v  :  Test  F31 
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Figure  11.125.  Turbulence  energy  spectrum  for  u  :  Test  1*32 
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Figure  I1.12G.  Turbulence  energy  spectrum  for  v  :  Test  F32 
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Figure  11.127.  Turbulence  energy  spectrum  for  «'*:  Ten 
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Figure  H.12S.  Turbulence  energy  spectrum  far  r’2;  Test  F33 
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Appendix  C.  HOT-  H ‘IRE  CA  UliliA T10.\‘ 


The  X-type  hot-wire,  also  called  the  X-wirc,  is  calibrated  for  (low  speed  and  angle  in 
the  turbine  cascade  wind  tunnel. 

The  geometry  of  a  cross-flow  X-wirc  is  given  in  figure  C.l.  This  geometry  is  defined 
arbitrarily  as  a  normal  configuration  for  this  type  of  X-wirc,  that  is,  wire  2  is  towards 
the  right  of  the  oncoming  flow  vector  as  viewed  in  the  direction  of  (/«>.  An  inverse 
configuration  of  the  wire  geometry  is  defined  as  wire  2  to  the  left  of  tf©5. 

C.l  Velocity  Calibration  Measurement* 

The  following  procedure  details  the  velocity  calibration  method  for  the  hot-wiri* 
probe. 

Initial  Probe  and  Probe  Holder  Set-up: 

The  hot-wire  probe  is  inserted  into  the  probe  holder  and  the  angle  of  the  bisector 
with  respect  to  the  probe  holder  is  recorded  or  set  to  0  degrees.  A  simple  protractor  with  a 
pointer  attached  to  the  probe  holder  is  used  for  this  measurement.  The  pilch  angle  (.f)  is 
measured  to  within  ±  0.5  degrees.  The  angle  of  the  probe  longitudinal  axis  with  respect  to 
the  probe  holder  longitudinal  axis  is  the  yaw  angle,  denoted  as  i'\  The  effect  of  varying  v 
on  the  X-wirc  signal  is  determined  to  be  negligible  for  this  type  of  geometry;  for  V'  angles 
between  0-1  degrees. 

The  probe  and  probe  holder  arc  inserted  through  the  wind  tunnel  end  wall  into  an 
area  where  the  free-stream  flow  angle  has  been  determined.  The  bisector  is  set  to  0  degrees 
with  respect  to  the  oncoming  flow  angle  (  (i  =  0  degrees).  Refer  to  Figure  C.l. 

The  Thermo  Systems  IFA-100  is  adjusted  for  the  proper  cable,  probe  holder,  and 
probe  resistance  for  each  wire  of  the  X-type  hot-wire.  Similarly,  the  operating  resis¬ 
tance  for  each  wire  is  set.  Refer  to  the  Thermo  Systems  IFA-100  operator’s  manual 
(TS1  Manual,  1983)  for  detailed  instructions  on  the  proper  setup  of  this  system. 

Speed  Measurement 


C-l 


=  local  flow  velocity  vector 
=  Wire  1  and  Wire  2,  respectively 
=  X-wire  bisector 

=  angle  subtended  by  D  and  either  X-wire 
=  angle  subtented  by  UK  and  B 
-  normal  to  U'j 
=  normal  to 

=  angle  subtended  by  U &  and  rtj 
=  angle  subtended  by  {/w  and  n2 
=  angle  subtended  by  B  and  the  reference  X-axis 

Figure  C.l.  End-flow  X-type  hot-wire  geometry 


The  wind  tunnel  is  started  and  the  free-stream  flow  is  allowed  to  come  to  steady  state 
conditions  at  the  lowest  speed  setting.  The  output  voltage  from  each  hot* wire  is  measured 
400  times  using  the  Hewlett-Packard  Integrating  Voltmeter,  Model  No.  44701  and  aver¬ 
aged  to  a  single  voltage  reading.  The  corresponding  differential  height  in  the  pitot  tube 
manometer  is  measured  and  recorded.  The  wind  tunnel  speed  is  increased  incrementally, 
and  hot-wire  output  signals  measured,  until  a  minimum  of  eleven  readings  per  hot-wire  are 
completed.  This  process  is  automated  using  the  Hewlett-Packard  1IP3852A  High-Speed 
Data  Acquisition  mainframe,  and  the  X-wirc  calibration  program  A" Will  EC  A L /M5 listed 
at  the  cud  of  this  appendix. 

C.2  Angle  Calibration  Measurement# 

Similar  to  the  speed  calibration  procedure,  the  hot-wire  is  calibrated  in  the  wind 
tunnel  at  the  same  location  used  for  X-wire  speed  calibration.  The  wind  tunnel  speed  is 
set  to  a  single  value.  A  tunnel  speed  is  chosen  which  is  representative  of  the  expected  test 
conditions  to  be  used  for  subsequent  pressure  and  heal  transfer  measurements.  The  probe 
is  rotated  from  -10  degrees  to  +10  degrees  in  1  degree  increments.  Hot-wire  output  voltages 
for  each  wire  are  measured  using  the  Hewlett-Packard  Integrating  Voltmeter,  Model  No. 
44701  A.  Each  hot-wire  output  signal  is  measured  400  limes  and  averaged  to  a  single  value 
of  output  voltage.  The  offset  angle  is  recorded  at  each  increment.  A  total  of  21  angle 
increments  and  their  corresponding  hot-wire  output  voltages  are  measured  and  recorded. 

A  QuickBASIC  program  A’ WIHECA LBAS  automates  the  angle  calibration  data 
collection  process  for  X-wire.  A  listing  of  X  WIR  EC  A  L.  BAS  is  included  at  the  end  of  this 
appendix  for  reference.  The  Hewlett-Packard  3852A  High-Speed  Data  Acquisition  System 
measures  and  stores  all  hot-wire  signal  values  in  its  internal  memory  butler.  Permanent 
data  storage  is  on  floppy  disk. 

C.3  Hot-Wire  Probe  Calibration  Equations 

The  governing  equations  for  the  X-wirc  are  given  as  follows: 


C-3 


=  ^(cos’tf,  +  kxsU\78i) 
Vi,!*  =  t/coJ(coss0j  +  i-jsin’tfa) 


(C\l) 

(C.2) 


where 

-  component 

a  component 

=  coefficients 

of  Ucc  perpendicular  to  ll'j 
of  Uoa  perpendicular  to  H'j 
for  wire  1  and  wire  2,  respectively 


The  signal  output  from  each  hot-wire,  in  terms  of  voltage,  is 


=  cj  +«jfrjf// 

=  Ci  +  M?3?//*  + 

(C.-l) 

where 

£1 

ft 

(>1,  b\  and  ci 
«2,  l>2  and  c.2 


-  predicted  wire  1  output  signal  (volts) 
=  predicted  wire  2  output  signal  (volts) 
=  wire  1  velocity  calibration  coefficients 
=  wire  2  velocity  calibration  coefficients 


The  free-stream  flow 
pitot  tube  manometer  as 


velocity  [Uoo)  is  computed  from  the  differential  height  of  the 


(/„=  (— Aa-249V 

\Poo  / 


(C.5) 


where 

Poo  =  free-stream  static  air  density 

A h  =  differential  height  of  manometer  in  inches  of  water 

249  =  conversion  factor  for  inches  of  water  to  kPa  conversion 


Using  Equation  C.5  in  conjunction  with  Equation  C.l  and  Equation  C.2  for  wires  1  and  2, 
respectively,  U\(U  and  U2eiJ  is  computed  for  each  wind  tunnel  speed  measured  with  the 
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liot*wirc  probe.  Values  for  k\  and  fcj  are  arbitrarily  specified  at  this  point  of  the  calibration 
procedure. 

A  quadratic  fit  is  computed  for  Equations  C.l  and  C.2  for  the  independent  variables 
U\  '/ft  and  using  matrix  multiplication  (MATHCAD  Manual,  19S7). 

Matrices  an^  arc  formed 


=  1 


m<J>  =  "■.„* 


B.a>  = 


(C.<i) 


U2<»  =  1 


ili<2>  = 


l!i<3>  =  Ur 


«// 


(C.7) 


where 

ili<,,3,3>  =  columns  1, 2.  and  3  of  Mjl  ,  respectively 

ll';<1,3,3>  =  columns  1,2,  and  3  of  U^,  respectively 


Then  the  coeliicients  «j  ,  b\  ,  and  ci  of  Equation  C.3  are  found  from  the  matrix 
product 


=  (inT-nir,-(«iT 


(C.S) 


where 

WiT  =  transpose  of  IVl 

E\  =  measured  hot-wire  output  voltages  for  wire  1 
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Coefficients  «2 ,  bi ,  and  c?  of  Equation  C.d  are  determined  as 


(0.9) 


where 

\\'iJ  =  transpose  of  Wj 

E-i  -  measured  hot-wire  output  voltages  for  wire  2 

Hence,  Equation  C.8  and  Equation  C.O  give  a  quadratic  fit  of  the  speed  calibration 
data  for  wire  1  and  wire  2,  respectively.  Substitution  of  the  coefficients  determined  from 
Equation  C.S  into  Equation  C.3  give  a  predicted  value  E\ 2.  Similarly,  coefficients  found 
in  Equation  C.O  are  substituted  into  Equation  C.»l  to  give 

The  sum  of  the  squared  error  between  predicted  and  measured  output  signal  for  all 
time  steps,  A',  is  computed  as 


« 

(C.1U) 

A* 

(C.ll) 

where 

A'  s=  total  number  of  speed  calibration  data  measurements 

for  wires  1  and  2 

Ei  =  measured  wire  1  output  voltage 

Ei  =  measured  wire  2  output  voltage 

Coefficients  k\  and  ibj  are  determined  from  the  angle  calibration  data  as  follows. 

Ucc  is  computed  from  Equation  C.5  using  the  tunnel  speed  setting  from  the  angle 
calibration  measurement.  From  Equation  C.l  and  Equation  C.2,  U and  are 
computed.  Equation  C.l  and  Equation  C.2  are  presented  below  for  reference. 
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^^(COS^+W^ 


Substitution  of  {/j,y/  and  into  Equation  C.3  and  Equation  C.*t,  respectively, 
yield  predictions  for  £|2  and  £aJ. 

The  sum  of  the  squared  error  between  the  predicted  output  voltage  and  the  output 
voltage  measured  at  each  offset  angle  increment  for  all  N  values  is  given  by 


=x>v-*r’ 

4 

(CM2) 

.V 

l 

(CM3) 

where 

N  =  total  number  of  angle  calibration  data  measurements  for  wires  1  and 

Ey"  =  measured  wire  1  output  voltage 

£5"  =  measured  wire  2  output  voltage 

Coefficients  k\  and  Jtj  are  chosen  so  as  to  minimize  the  sum  of  the  squared  error  as  calcu¬ 
lated  in  Equations  C.12  and  C.13. 

Actual  implementation  of  the  above  X*wirc  calibration  is  coded  using  MATH CAD, 
version  2.1.  MAT11CAD  allows  the  user  to  sec  plots  of  quadratic  fits  for  both  wire 
equations  while  simultaneously  choosing  coefficients  and  ki  which  minimize  the  error. 

A  sample  printout  of  a  hot-wire  calibration  using  MATJICAD  is  included  in  this  appendix. 
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Sample  calibration  of  a  hot-wire  probe  using  MATH  CAD.  version  *J.l 
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Sample  calibration  of  a  liot*wirc  probe  using  MATllCAl),  version  *2.1 
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Sample  calibration  ©f  a  hot-wire  probe  umiir  MATH  CM),  version  ’2.1 
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C.5  X-TyjK  Hot-Wire  Calibration  Source  Code  Lifting 


*  Hot  wire  (X  type  calibration  program) 

*  High  speed  voltmeter  is  usad  for  signal  measurement. 

*  Written  by  Lallo  Calassi  7/1/89. 

* 

'  This  calibration  program  raads  tha  voltaga  output  by  aach  vira 
'  during  calibration  tasting.  Each  channal  is  raad  500  times  than 

*  averaged  for  a  mean  voltaga. 

I 

*  Tha  program  loops  and  triggara  on  prassing  any  hay. 

CLS 

PRINT . . . " 

PRINT  "X-wira  calibration  program  for  tha  HP  38S2A" 

PRINT  »  HP44701  Intagrating  Voltmatar  in  slot  0600" 

PRINT  "  Disconnact  tha  dadicatad  high-spaad  FET  cabla  from  tna  HP44702B" 
PRINT  SLEEP  4 

5  CLS  :  INPUT  "Is  this  a  (V)alocity  or  (A)ngla  type  calibration?  ",  var$ 
INPUT  "Enter  wire  serial  no.  (data  stored  on  A:\SRI\v.dat)?  ",  dir! 

INPUT  "Filename  to  store  data  on  disk  A?  ",  filename! 

OPEN  "A:\"  +  dir!  +  "\"  +  filename!  +  ".dat"  FOR  OUTPUT  AS  «1 
IF  (var!  «  "V")  THEN 

INPUT  "Current  air  density?  ",  rho! 

END  IF 

PRINT  "a**********************************************" 

PRINT  "" 

VIEW  PRINT  8  TO  22 
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CALL  ibfind("hp3852",  dvn‘/.) 

CALL  ibvrt(dvmX,  "rst") 

10  PRINT  "Press  any  key  to  begin  scan  11 
DO:  LOOP  WHILE  INKEY*  «  "" 

CALL  iburt(dvmX,  "disp  off") 

CALL  iburt(dvmy,,  "real  A<399) ,B(399) ,L,H,AV,S") 

*  This  call  configures  the  HP44701  Integrating  Voltmeter  as  the  measuring 
'  instrument. 

CALL  iburt(dvmX,  "use  600;nplc  1") 

CALL  iburt(dvm'/,,  "conf  dcvjnrdgs  400;azoro  once") 

’  Procede  uith  channel  1  and  channel  2  voltage  scan. 

CALL  ibvrt(dvmy,,  "meas  dev, 321, use  600, into  A;conf  dcv;nrdgs  400") 

CALL  iburt(<lvm‘/.,  "meas  dev, 322, use  600, into  B") 

’Perform  statistcal  averaging  on  the  measurements  stored  in  arrays 

*  A  and  B. 

CALL  ibwrt(dvm*/,,  "stat  L,H,AV,S,A") 

CALL  ibwrt(dvm*/,,  "vread  AV") 

CALL  iburt(dvm‘/,,  "disp  on;  disp  AV") 
rd$  =  SPACE$(16) 

CALL  ibrd(dvm'/,,  rd$) 
chanl!  «  VAL(rd$) 

'  Output  to  the  screen  and  to  the  data  storage  file 
COLOR  0,  5 
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PRINT  "channel  1  voltage  i*:",  chan* ! 

CALL  ibwrt(dvny,,  "stat  L,H,AV,S,B") 

CALL  ibvrt(dvmX,  "vread  AV") 

CALL  ibrd(dvnX,  rdt) 
chan2!  ■  VAL(rd$) 

PRINT  "channel  2  voltage  is:",  chan2! 

IF  (var$  «  "V")  THEN 

INPUT  "Input  the  currant  velocity  setting  ",  vel! 

ELSE 

INPUT  "Input  the  current  angle  setting  ",  vel! 

END  IF 

IF  (var$  «  "V")  THEN 
vel!  »  SQR(vel!  ♦  249!  *  2!  /  rho) 

PRINT  "current  air  speed  is:  ",  vel!,  "  m/sec" 

ELSE 

PRINT  "current  angle  setting  is:",  vel!,  "degrees  from  bisector" 
END  IF 

*  Output  to  data  storage  file  on  disk  A: 

PRINT  il,  vel!,  chanl!,  chan2! 

INPUT  "quit?  (y  or  n)  ",  q$ 

IF  q$  <>  "y"  THEN 
GOTO  10 
ELSE 
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CLOSE  •  ! 


PRINT  "This  stags  of  calibration  completed  -  data  stored  in",  filonamoS 
PRINT  "" 

END  IF 

INPUT  "Do  an  angle  calibration  for  the  current  hot  wire?  (y)  or  (n)",  rst$ 
IF  (rst$  -  "y")  THEN 
CLS  2 
GOTO  5 
END  IF 
END 
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Appendix  D.  ADDITIONAL  II  DAT  Til  ANSI- DR  TESTS— STA  NTON  NUMBER 
RATIOS  AND  NUSSDLT  NUMBER  RATIOS 
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Nusselt  number  on  the  turbine  blade  for  co-flow  injection.  D  =  O.GG  mm. 
Test  All,  A 12  and  A 13.  Jet-grid  at  station  0 
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Figure  D.2.  Stanton  number  on  the  turbine  blade  for  co-flow  injection.  D  =  1.32  mm. 
Test  A21,  A22  and  A23.  Jet -grid  at  station  0 
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Nusselt  number  on  the  turbine  blade  for  co-flow  injection.  I)  =  1.32  mm. 
Test  A21,  A22  and  A23.  Jet-grid  at  station  0 
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Figure  D.4.  Stanton  number  on  the  turbine  blade  for  co-flow  injection.  D  =  1.9S  mm. 
Test  A31,  A32  and  A33.  Jet -grid  at  station  0 
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Figure  D.5.  Stanton  number  on  the  turbine  blade  for  cross-flow  injection.  D  -  1.9$  mm. 
Test  1)31, 1)32  and  1)33.  Jet-grid  at  station  0 
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Figure  D.C.  Stanton  number  on  the  turbine  blade  for  co-flow  injection.  D  =  0.66  mm. 
Test  Dll,  D12  and  D13.  Jet-grid  at  station  02 
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Figure  D.8.  Stanton  number  on  the  turbine  blade  for  cross-flow  injection.  D  =  O.GG  nun. 
Test  Ell,  E12  and  E13.  Jet-grid  at  station  02 
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Nussclt  number  on  the  turbine  blade  for  cross-flow  injection.  0  =  O.CG  mm. 
Test  Ell,  E12  and  EI3.  Jet-grid  at  station  02 
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Figure  D.10.  Nusselt  number  on  the  turbine  blade  for  counter-flow  injection.  D  =  0.66 
mm.  Test  Fll,  F12  and  F13.  Jet-grid  at  station  02 
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Stanton  number  on  the  turbine  blade  for  cross-flow  injection.  1)  =  1.32 
mm.  Test  1221, 1222  and  1223.  Jet-grid  at  station  02 
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Figure  D.12.  Nusselt  number  on  the  turbine  blade  for  cross- flow  injection.  D  =  1.32  mm. 
Test  E21,  E22  and  E23.  Jet-grid  at  station  02 
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Figure  D.M.  Stanton  number  on  the  turbine  blade  for  co-flow  injection.  D  -  1.9S  mm. 
Test  D31,  D32  and  D33.  Jet-grid  at  station  02 
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Figure  D.1G.  Nusselt  number  on  the  turbine  blade  for  cross -flow  injection.  1)  -  1.98  mm. 
Test  E31,  E32  and  E33.  Jet-grid  at  station  02 
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Appendix  E. 


DATA  ACQUISITION  SOURCE  COOK  LISTING 


HP3S52A  Data,  Acquisition  System  Source  Code 
Written  by  lit  Lcllo  Galassi,  1959 
Source  code  language:  Microsoft  QuickBASIC ,  version  *1.5 


'  HP3852A  Data  Acquisition  Program  for  High  Turbulence  Hoat  Transfer 
'  (for  use  on  the  Linear  Turbine  Cascade  Wind  Tunnel,  BLDG  19) 

I 

’  Written  by  Lt  Lello  Galassi,  6/23/89 


1  QuickBASIC  Declarations  Rqv.  C.4 


1  Common  GPIB  status  variables 


COMMON  SHARED  ibstaX,  iberr'/,,  ibent X 


*  GPIB  Subroutine  Declarations 

DECLARE  SUB  IBBNA  (BD*/.,  BDNAMES) 
DECLARE  SUB  IBCAC  (BD*/.,  V*/.) 
DECLARE  SUB  IBCLR  (BD’/.) 

DECLARE  SUB  IBCMD  (BD*/. ,  CMD$) 
DECLARE  SUB  IBCMDA  (BD*/.,  CMD$) 
DECLARE  SUB  IBDMA  (BD*/.,  V*/.) 
DECLARE  SUB  IBEOS  (BD/.,  V'/.) 
DECLARE  SUB  IBEOT  (BD/. ,  VX) 
DECLARE  SUB  IBFIND  (BDNAMEJ,  BD*/.) 
DECLARE  SUB  IBGTS  (BDX,  V*/.) 
DECLARE  SUB  IBIST  (BOX ,  VX.) 
DECLARE  SUB  IBLOC  (BDX) 

DECLARE  SUB  IBONL  (BDX,  VX) 
DECLARE  SUB  IBPAD  (BDX,  VX) 
DECLARE  SUB  IBPCT  (BDX) 

DECLARE  SUB  IBPPC  (BDX,  VX) 
DECLARE  SUB  ibrd  (BD'/,,  rd$) 
DECLARE  SUB  IBRDA  (BD'/,,  rd$) 
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DECLARE  SUB  IBRDF  (BDX,  FLNAHES) 

DECLARE  SUB  IBRDI  (BDX ,  IARRXO,  CNTX) 
DECLARE  SUB  IBRDIA  (BD*/.,  IARRXO ,  CNT*/,) 
DECLARE  SUB  IBRPP  (BDX,  PPRX) 

DECLARE  SUB  IBRSC  (BOX,  VX) 

DECLARE  SUB  IBRSP  (BOX,  SPRX) 

DECLARE  SUB  IBRSV  (BDX,  VX) 

DECLARE  SUB  IBSAD  (BDX,  VX) 

DECLARE  SUB  IBSIC  (BDX) 

DECLARE  SUB  IBSRE  (BDX,  VX) 

DECLARE  SUB  IBSTOP  (BDX) 

DECLARE  SUB  IBTMO  (BDX,  VX) 

DECLARE  SUB  IBTRAP  (HASKX ,  MODEX) 
DECLARE  SUB  IBTRG  (BDX) 

DECLARE  SUB  IBWAIT  (BDX,  MASKX) 

DECLARE  SUB  IBWRT  (BDX,  VRTS) 

DECLARE  SUB  IBWRTA  (BDX,  VRT$) 

DECLARE  SUB  IBWRTF  (BDX,  FLNAMES) 
DECLARE  SUB  IBWRTI  (BOX.  IARRXO,  CNTX) 
DECLARE  SUB  IBWRTIA  (BDX,  IARRXO,  CNTX) 


'  GPIB  Function  Declarations 

DECLARE  FUNCTION  ILBNAX  (BDX,  BDNA«E$) 

DECLARE  FUNCTION  ILCACX  (BDX,  VX) 

DECLARE  FUNCTION  ILCLRX  (BDX) 

DECLARE  FUNCTION  ILCMDX  (BDX.  CHD$,  CNTX) 
DECLARE  FUNCTION  ILCMDAX  (BDX,  CHD*.  CNTX) 
DECLARE  FUNCTION  ILDMAX  (BDX,  VX) 

DECLARE  FUNCTION  ILEOSX  (BDX,  VX) 

DECLARE  FUNCTION  ILEOTX  (BDX,  VX) 

DECLARE  FUNCTION  ILFINDX  (BDNAHEJ) 

DECLARE  FUNCTION  ILGTSX  (BDX,  VX) 

DECLARE  FUNCTION  ILISTX  (BDX,  VX) 

DECLARE  FUNCTION  ILLOCX  (BDX) 

DECLARE  FUNCTION  ILONLX  (BDX,  VX) 

DECLARE  FUNCTION  ILPADX  (BDX,  VX) 

DECLARE  FUNCTION  ILPCTX  (BDX) 

DECLARE  FUNCTION  ILPPCX,  (BDX,  VX) 

DECLARE  FUNCTION  ILRDX  (BDX,  rd$,  CNTX) 
DECLARE  FUNCTION  ILRDAX  (BDX,  rd$,  CNTX) 
DECLARE  FUNCTION  ILRDFX  (BDX,  FLNAHE$) 
DECLARE  FUNCTION  ILRDIX  (BDX,  IARRXO,  CNTX) 
DECLARE  FUNCTION  ILRDIAX  (BDX,  IARRXO,  CNTX.) 
DECLARE  FUNCTION  ILRPPX  (BDX,  PPRX) 
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DECLARE  FUNCTION  ILRSCX  (BOX ,  V/.) 

DECLARE  FUNCTION  ILRSP7,  (BD7, ,  SPRX) 

DECLARE  FUNCTION  ILRSVX  (BDX,  VX) 

DECLARE  FUNCTION  ILSADX  (BDX,  VX) 

DECLARE  FUNCTION  ILSICX  (BDX) 

DECLARE  FUNCTION  ILSREX  (BDX,  VX) 

DECLARE  FUNCTION  ILSTOPX  (BDX) 

DECLARE  FUNCTION  ILTNOX  (BDX,  VX) 

DECLARE  FUNCTION  ILTRAPX  (MASKX,  HODEX) 

DECLARE  FUNCTION  ILTRGX  (BDX) 

DECLARE  FUNCTION  ILWAITX  (BDX,  MASKX) 

DECLARE  FUNCTION  ILWRTX  (BDX,  WRT$,  CNTX) 

DECLARE  FUNCTION  ILWRTAX.  (BDX,  VRT$.  CNTX) 

DECLARE  FUNCTION  ILWRTFX  (BDX,  FLNAHEJ) 

DECLARE  FUNCTION  ILWRTIX  (BDX.  IARRXO,  CNTX) 

DECLARE  FUNCTION  ILWRTIAX  (BDX,  IARRXO,  CNTX) 

DECLARE  SUB  TMON  (dvmX) 

DECLARE  SUB  GPIBERil  () 

DECLARE  SUB  PRUN  (dvmX.  rho!,  dtt,  FileName*,  iflagX,  Temp!) 
DECLARE  SUB  FRUN  (dvmX,  FileNamet) 

DECLARE  SUB  TRUN  (dvmX,  Temp!,  FileName$,  rho!) 

DECLARE  SUB  CONVERT  (VALS! () ,  NNX) 

DECLARE  SUB  FFT  () 

DECLARE  SUB  PLOTT  (xdatl!(),  ydatl!(),  xmin!,  xmax!,  ymin!,  ymax 
xlabj,  ylab$,  pltlab$,  ntorX) 

DECLARE  SUB  SECPRES  (dvmX,  sfa!) 


DIM  chord(24)  AS  SINGLE 


*  Initialize  the  screen 

CLS  :  COLOR  3,0:  LOCATE  7,  12 
PRINT  "HP3852  Controller  Program":  PRINT 
LOCATE  8,  12 
LOCATE  10,  12 

PRINT  "This  program  controls  the  HP3852  Data  Acquisition  System" 
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LOCATE  11,  12 


PRINT  "for  high-turbulence  heat  transfer  measurements." 


PRINT  "" 


SLEEP  5 

CLS  :  LOCATE  7,  7 

INPUT  "Bypass  turn-on  sequence?  ",  ans$ 

IF  ans$  -  "y"  THEN 
GOTO  SS 
END  IF 

PRINT  :  PRINT  "Initialization  of  Program:":  PRINT 

10  INPUT  "HP3852  on?  (y  or  n)",  ans* 

IF  ans$  <>  "y»  THEN 
BEEP 

PRINT  "TURN  ON  HP3852" 

GOTO  10 
END  IF 

20  INPUT  "Blade  and  Scanivalve  power  supplios  on?  (y  or  n)",  ans$ 

IF  anst  <>  "y"  THEN 

BEEP 

PRINT  "TURN  ON  BOTH  POWER  SUPPLIES" 

GOTO  20 
END  IF 


30  INPUT  "Scanivalve  on?  (y  or  n)",  ans$ 

IF  ans$  <>  "y"  THEN 

BEEP 

PRINT  "TURN  ON  SCANIVALVE" 

GOTO  30 
END  IF 


40  INPUT  "IFA  100  on  and  channels  set?  (y  or  n)",  ans$ 

IF  ans$  <>  "y"  THEN 

BEEP 

PRINT  "SET-UP  THE  IFA  100" 

GOTO  40 
END  IF 
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> 


'Read  the  blade  chord  locations  for  pressure  ports  and  thermocouples 
SS  OPEN  "C:\LELLO\chord.dat"  FOR  INPUT  AS  «2  ‘blade  chord  locations  aro 
FOR  J  *  i  TO  23  ‘stored  into  the  data  file 

INPUT  *2,  chord(J)  ‘chord.dat  on  the  C:  drive. 

'PRINT  "chord(",  J,  CHORD(J) 

NEXT  J 
CLOSE  12 


’  Initial  Test  parameters: 

CLS  :  LOCATE  7,  7 

PRINT  :  PRINT  "Data  initialization:":  PRINT 

INPUT  "Run  Designation  for  all  data  files  ?  ",  FileNamc$ 

iflag’/.  *  0 

dt$  *  DATE$ 

time*  «  TIMER 

CLS  :  LOCATE  7,  7 

INPUT  "Current  temperature  (deg.  F)7  ",  Deg! 

Temp!  «  (Deg!  -  32!)  *  5!  /  9! 

rholow!  =  1.3947 

drho!  =  .4666 

Tlov!  =  23.15 

dTemp!  =  100! 
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CLS  :  LOCATE  7,  7 


rho!  *  rholou!  -  drho!  •  (Temp!  ♦  Tlow!)  /  dTemp! 

PRINT  "Current  air  density  is  rho!;  SPC(l);  "kg/m"3" 

LOCATE  12,  7 

PRINT  "  Press  any  key  to  continue..." 

DO:  LOOP  WHILE  INKEY$  «  "" 

CLS  :  LOCATE  7,  7 

PRINT  :  PRINT  "Data  Initialization  Complete." 

*****tt*  +  **Mt*#**t  +  44#***t**t>Mt*t#  +  *#*  +  + ft*********##*# +  +  +  #*## 

’Assign  a  unique  identifier  to  the  HP38S2  for  use  with  the  GPIB 
CALL  IBFIND("HP3852" ,  dvm*/.) 

IF  dvmy.  <  0  THEN  CALL  GPIBERR 
CALL  IBWRT(dvm'/.,  "RST") 

’Routine  Branching  for  the  MAIN  HEMU 
DO 

CLS  :  PRINT  :  PRINT  "  MAIN  MENU":  PRINT 

PRINT  STRINGS (80,  PRINT 

PRINT  "01)  Temperature  Run" 

PRINT  "02)  Pressure  Run" 

PRINT  "03)  Turbulence  Frequency  Run" 

PRINT  "04)  Temperature  and  Frequency  Run" 

PRINT  "05)  Temp.,  Freq.,  and  Press.  Run" 

PRINT  "06)  Frequency  and  Pressure  Run" 


PRINT  "07)  Fast  Fourier  Transform  of  Data" 


PRINT  "08)  Currant  Temperature  of  tha  blada" 

PRINT  "09)  Secondary  Prassura  in  Manifold" 

PRINT  "10)  EXIT" 

PRINT  :  PRINT  "Type  your  selection  (01-10):" 
chs*  -  INPUT* (2) 

'  Use  SELECT  CASE  to  process  response. 

SELECT  CASE  chs* 

CASE  "01":  CALL  TRUNCdvmy,,  Tamp!,  FileName*,  rho!) 

CASE  "02":  CALL  PRUN(dvmy.,  rho!,  dt*.  FileName*,  iflag'/., 
Temp!) 

CASE  "03":  CALL  FRUN (dvm‘/, ,  FileName*) 

CASE  "04":  CALL  TRUNCdvmy.,  Temp!,  FileName*,  rho!) 

CALL  FRUN  (dvm'/, ,  FileName*) 

CASE  "05":  CALL  TRUNCdvmy.,  Temp!,  FileName*,  rho!) 

CALL  PRUN (dvmy, ,  rho!,  dt*,  FileName*,  iflag'/., 
Temp!) 

CALL  FRUN  (dvm'/. ,  FileName*) 

CASE  "06":  CALL  PRUN(dvmy.,  rho!,  dt*,  FileName*,  iflag'/., 
Temp!) 

CALL  FRUNCdvmy.,  FileName*) 

CASE  "07":  CALL  FFT 
CASE  "08":  CALL  TM0N(dvmy.) 

CASE  "09":  CALL  SECPRES(dvm*/,,  sfa!) 
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CASE  "10":  EXIT  DO 


CASE  ELSE 


BEEP 

END  SELECT 

LOOP 

END 

SUB  CONVERT  <VALS!(),  NNX) 

DEFINT  I-M 
DIM  VALS! (2049) 

N  ■  2  •  NN 
J  «  1 

FOR  i  «  1  TO  N  STEP  2 
IF  J  <*  i  THEN  GOTO  2200 
TEMPR  «  VALS(J) 

TEMPI  *  VALSCJ  +  1) 

VALS(J)  «  VALS(i) 

VALSCJ  +  1)  «  VALSCi  +  1) 

VALS(i)  «  TEMPR 
VALSCi  +  1)  «  TEMPI 
2200  M  *  N  /  2 

2210  IF  M  >*  2  AND  J  >  M  THEN  GOTO  2240 
J  *  J  +  M 
GOTO  2310 
2240  J  *  J  -  M 
M  *  M  /  2 
GOTO  2210 
2310  NEXT  i 
MMAX  =  2 

2330  IF  N  <=  MMAX  THEN  GOTO  2700 
ISTEP  *  2  *  MMAX 

THETA*  «  6.28318530717959#  /  ISIGN  /  MMAX 
WPR#  *  -2  *  CSIN (THETA#  /  2)  "  2) 

WPI#  «  SINCTHETA#) 

WR«  *  1# 

WI#  «  0.000000000000000 
FOR  M  «  1  TO  MMAX  STEP  2 
FOR  i  -  M  TO  N  STEP  ISTEP 
J  «  i  +  MMAX 

TEMPR  «  CSNGCWR#)  *  VALS(J)  -  CSNG(WI#)  *  VALSCJ  +  1) 
TEMPI  =  CSNGCWR#)  *  VALSCJ  +  1)  +  CSNG(WI#)  *  VALS(J) 
VALSCJ)  =  VALSCi)  -  TEMPR 
VALSCJ  +  1)  =  VALSCi  +  1)  -  TEMPI 
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VALS(i)  -  VALS(i)  +  TEMPR 
VALS(i  +  1)  «  VALS(i  +  1)  4  TEMPI 
NEXT  i 

WTEMPi  »  WRi 

WRi  «  WR>  *  WPRi  -  WI»  ♦  WPI#  4  WRi 
WI#  -  VIC  •  WPRi  4  WTEMPi  •  WPIi  4  VIC 
NEXT  M 

MMAX  «  ISTEP 
GOTO  2330 

2700  'ALL  DATA  IS  NOW  CONVERTED  AND  IS  READY  TO  BE  STORED 
IF  ISIGN  «  -1  THEN  RETURN 
FOR  i  ■  1  TO  NN  *  2 
VALS(i)  «  VALS(i)  /  NN 
NEXT  i 

END  SUB 

DEFSNG  I-M 

SUB  DVMERROR  STATIC 

PRINT  "DVM  GPIB  ERROR" 

END  SUB 

SUB  FFT 

•THIS  ROUTINE  TRANSFORMS  DATA  FROM  THE  TIME  TO  FREQUENCY  DOMAIN 
'OR  BACK  USING  THE  DANIELSON-LANCZOS  FFT  METHOD 
•THE  DATA  FORMAT  REQUIRED  FOR  READING  IS  AS  FOLLOWS: 

'  TOP  LINE: 

'  NN  C  R  ISIGN 

•WHERE  NN  *  •  OF  SAMPLED  DATA  POINTS  <  1025 

'  C  «  •  OF  CHANNELS  READ 

•  R  *  (1)  IF  REAL  DATA,  (2)  IF  COMPLEX  DATA,  (3)  IF  MAGNITUDE  DATA 

'  ISIGN*  (41)  TF  DATA  IS  IN  TIME  DOMAIN,  (-1)  IF  DATA  IS  IN  FREQ  DOMAIN 

I 

•  REMAINING  DATA: 

»  TIME  OR  FREQ  CHI  CH2  . * 


>*********************  ppx  ROUTINE  ********************* 
<***44*44444****4**44  MAIN  PROGRAM  ************************ 
DIM  VALS! (2049) 

DIM  VALU! (2049,  6) 

DIM  VALM!(10?r>,  6) 

'GET  DATA 

PRINT  "AVAIL  DATA  FILES:" 

FILES  "\LELLO\*.FFT" 

INPUT  "INPUT  DATA  FILE  NAME  W/0  EXTENSION  (.FFT):  ";  ND$ 
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OPEN  "\LELLO\"  +  ND$  +  ".FFT"  FOR  INPUT  AS  «1 

INPUT  *1,  NN,  C,  R,  ISIGN 

Sl$  «  "CONVERTING  FROM  TIME  TO  FREQ  DOMAIN" 

S2$  «  "CONVERTING  FROM  FREQUENCY  TO  TIME  DOMAIN" 

IF  R  «  3  THEN 
PRINT  "NO  CAN  DO" 

GOTO  3000 
END  IF 

IF  ISIGN  *  1  THEN  PRINT  Sl$  ELSE  PRINT  S2$ 

FOR  i  «  1  TO  NN  *  R 
FOR  J  -  1  TO  C  +  1 
INPUT  fll,  VALU(i,  J) 

NEXT  J 

NEXT  i 
CLOSE  *1 

IF  R  «  1  THEN  DELTA  *  1  /  NN  /  (VALU(2,  1)  -  VALU(1,  1)) 

IF  R  *  2  THEN  DELTA  *  1  /  NN  /  (VALU(3,  1)  -  VALU(1,  1)) 

•'STRIP  DATA 
FOR  JJ  *  2  TO  C  ♦  1 
CC  «  1 

FOR  i  «  1  TU  NN  *  2  STEP  2 

VALS(i)  *  VALU(CC,  JJ) 

IF  R  *  1  THEN  VALS(i  +  1)  ■  0  ELSE  VALS(i  +  1)  *  VALU(i  +  1,  JJ) 

CC  *  CC  +  R 
NEXT  i 
CC  =  1 

'GO  TO  FFT  SUBROUTINE 
CALL  CONVERT(VALS!  ()  ,  NN’/.) 

'REPLACE  DATA  (UNSTRIP) 

FOR  i  *  1  TO  NN  *  2 
VALU(i,  JJ)  *  VALS(i) 

NEXT  i 

'GET  MAGNITUDE  DATA 
FOR  i  »  1  TO  NN  *  2  STEP  2 

VALMCCC,  JJ)  «  SQR(VALS(i)  *  VALS(i)  +  VALS(i  +  1)  *  VALS(i  +  1)) 

CC  *  CC  +  1 
NEXT  i 

NEXT  JJ 

I 

'DATA  STORAGE  ROUTINE 

'DATA  STORED  IN  TWO  SEPARATE  FILES  BASED  ON  THE  ORIGINAL  FILE  NAME 
'THE  TWO  FILES  WILL  BE  AS  FOLLOWS: 

’  (1)  NAME.TSD  (TIME-SPLIT  DATA)  OR  NAME.FSD  (FREQUENCY  SPLIT  DATA) 
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»A  (2)  NAME.TMD  (TIME  MAGNITUDE  DATA)  OF  NAME.FMD  (FREQ  MAG  DATA) 

i 

‘STEP  1:  CONVERT  FIRST  COLUMN  (TIME  <— >  FREQUENCY) 

VALM(1,  1)  «  0 
VALU(1,  1)  ■  0 
VALU(2,  1)  «  0 
QTY  «  DELTA 
CC  «  2 

FOR  i  **  3  TO  2  ♦  NN  STEP  2 
VALU(i,  1)  ■  QTY 
VALU(i  +  1,  1)  «  QTY 
VALM(CC,  1)  «  QTY 
QTY  «  QTY  *  DELTA 
CC  «  CC  ♦  1 
NEXT  i 

i 

•STEP  2:  STORE  DATA 

IF  ISIGN  »  1  THEN  EXT1$  *  "-FM.DAT"  ELSE  EXT1S  *  "-TM.DAT" 

IF  ISIGN  *  1  THEN  EXT2S  *  "-FS.DAT"  ELSE  EXT2$  *=  "-TS.DAT" 

OPEN  "\LELLO\KFT\"  +  ND$  +  EXT2$  FOR  OUTPUT  AS  *2 
OPEN  "\LELLO\FFT\"  +  ND$  +  EXT1$  FOR  OUTPUT  AS  *3 
PRINT  «2.  NN;  C;  2;  ISIGN  *  (-1) 

PRINT  13,  NN;  C;  3;  ISIGN  *  (-1) 

FOR  i  «  1  TO  NN  *  2 
FOR  J  *  1  TO  C  +  1 
PRINT  #2,  VALU(i,  J); 

NEXT  J 
PRINT  *2, 

NEXT  i 
CLOSE  #2 
FOR  i  =  1  TO  NN 
FOR  J  *  1  TO  C  +  1 
PRINT  13,  VALM(i,  J); 

NEXT  J 
PRINT  13, 

NEXT  i 
CLOSE  #3 
3000  END  SUB 


SUB  FINDERR  STATIC 
PRINT  "IBFIND  ERROR" 

END  SUB 

SUB  FRUN  (dvm'/, ,  FileName*) 
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*  This  subroutine  performs  the  voltage  scan  on  the  hot  wire  for 

'  conversion  to  velocities.  Tvo  channels  are  scanned  consecutively 

*  with  the  HP447023  and  the  FET  Multiplexer.  The  FET  is  connected 

*  internally  to  the  HP44702B  via  a  ribbon  cable. 

CLS 

DIM  chanl(2047)  AS  SINGLE,  chan2(2047)  AS  SINGLE 
DEFINT  I-M 


LOCATE  7,  7 

INPUT  "Use  the  same  root  filename  as  in  MAIN  ?  ",  ans$ 

IF  ans$  <>  "y"  THEN 

CLS  :  LOCATE  7,  7:  INPUT  "  Input  the  new  file  name  ",  FileNamoJ 
END  IF 

OPEN  "A:\"  +  "TU"  +  FileName*  +  ".dat"  FOR  OUTPUT  AS  #1 
COLOR  0,  3 

CALL  SECPRES(dvm'/,,  sfa!) 

CLS  :  LOCATE  10,  7 

PRINT  "  RIBBON  CABLE  CONNECTED  INTERNALLY  TO  FET  MULTIPLEXER?" 

PRINT  : 

PRINT  "  WARNING!  TURN  OFF  HP3852A  PRIOR  TO  CONNECTING  RIBBON  CABLE" 
PRINT  "" 

PRINT  "  Press  any  key  to  initiate  scan" 

DO:  LOOP  WHILE  INKEY$  «  "" 

CLS  :  LOCATE  7,  7:  PRINT  "Turbuence  scan  in  process" 


CALL  IBWRTCdvmy..  "RST  400") 

CALL  IBWRTCdvmX,  "USE  400") 

CALL  IBWRTCdvmy, ,  "REAL  WAVE(4095) .0UT1C2047) .0UT2C2047)") 
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CALL  IBWRTCdva?.,  "DISP  OFF") 


i 


' Configure  the  HP44702B  ifor  high-apeed  acanning. 

CALL  IBWRTCdvaX,  "SCAMMODE  0N;C0NF  DCV;ARMODE  BEFORE") 

CALL  IBWRTCdvaX,  "TERM  RIBBON;range  9;RDGSM0DE  COMPLETE") 

CALL  IBWRT (dvm% ,  "SCDELAY  0") 

CALL  IBWRT(dvn7,,  "SPER  10E-6") 

CALL  IBWRTCdvmy,,  "PRESCAN  2048; POSTSCAN  0") 

CALL  IBWRTCdvmy.,  "CLWRITE  SENSE, 321-322;ASCAN  ON;SCTRIG  SGL") 

CALL  IBWRTCdvmy.,  "XRDGS  400,  4096  INTO  WAVE")  'Trigger  tho  scan 


CALL  IBWRTCdvm*/,,  "SUB  SEPARAT")  'Separate  the  readings  into  2  arrays 
CALL  IBWRTCdvm'/,,,  "INTEGER  I,J") 

CALL  IBWRTCdvm'/,,  "J  =  0") 

CALL  IBWRT(dvm‘/.,  "FOR  I  *  0  TO  4095  STEP  2") 

CALL  IBWRTCdvm'/,,  "0UT1(J)  =  WAVE(I)") 

CALL  IBWRTCdvmy,,  "0UT2(J)  =  WAVE(I+1)") 

CALL  IBWRTCdvmy,,  "J  *  J+l") 

CALL  IBWRTCdvmy.,  "NEXT  I") 

CALL  IBWRTCdvmy,,  "SUBEND") 

> ****************************************************************** 


CLS  :  LOCATE  7,  7:  PRINT  "Processing  Data" 


I 


CALL  IBURT(dvmX,  "CALL  SEPARAT") 
CALL  IBWRTCdvnX,  "SCTR1G  HOLD") 


CLS  :  LOCATE  7,  7:  PRINT  "Writing  data  to  temporary  disk  storage" 


♦output  to  the  controller 


’Transfer  all  readings  to  the  computer  controller. 


CALL  IBWRTCdvro*/., 
CALL  IBRDF(dvm‘/,, 
CALL  IBWRT(dvm'/, , 
CALL  IBRDFCdvm*/., 


"VREAD  0UT1") 
"C:\LELL0\FREQ1.DAT") 
"VREAD  0UT2") 
"C:\LELL0\FREQ2.DAT") 


'Read  from  temporary  files  into  memory. 

OPEN  "C:\LELL0\FREQ1.DAT"  FOR  INPUT  AS  13 
OPEN  "C:\LELL0\FREQ2.DAT"  FOR  INPUT  AS  #4 
FOR  i  ■  0  TO  2047 
INPUT  «3,  chanl(i) 

INPUT  #4,  chan2(i) 

NEXT  i 
CLOSE  #3 
CLOSE  #4 

’ ♦*****♦♦♦**♦*****♦****♦#***♦♦#*♦♦*♦*♦*♦******♦♦*•***♦♦**♦*«**♦«•**«* 


CLS  :  LOCATE  7,  7 

PRINT  "Writing  data  file  to  disk  A:" 


time!  =  0! 

FOR  i  *  0  TO  2047 

PRINT  #1,  time!,  chanl(i),  chan2(i) 
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time!  ■  tint!  +  .00002 
KEXT  i 
CLOSE  II 


CALL  IBWRT(dvmX,  "DISP  ON;RST  400”) 


CLS  :  LOCATE  7.  7 

PRINT  "Turbulence  data  acquisition  completed" 

LOCATE  9,  8 

PRINT  "DISCONNECT  THE  RIBBON  CABLE  IF  FINISHED" 

PRINT  :  SLEEP  2:  CLS  :  LOCATE  7,  7 

PRINT  "press  any  key  to  return  to  main  menu" 

DO:  LOOP  WHILE  INKEY$  «  "" 

CLS 

END  SUB 

DEFSNG  I-M 

SUB  GPIBERR  STATIC 

PRINT  "GPIB  ERROR" 

END  SUB 

SUB  PLOTT  (xd&t()  AS  SINGLE,  ydat()  AS  SINGLE,  xmin!,  xraax!,  ymin!, 
ymax!,  xlab$,  ylab$,  pltlabt,  ntotX) 

>44**** ****************** ********************************************* 

’Subroutine  PLOTT 

’This  subroutine  creates  a  plot  of  ydat  vs  xdat 
CLS 

SCREEN  9  ’Hi-res  graphics  mode 
VIEW  (120,  10)-(570,  290),  ,  1 
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WINDOW  (xain!,  ynin!)-(xmax! ,  yaax!) 

style*  *  ItHFFOO 

LOCATE  23,  42:  PRINT  pltlab$ 

LOCATE  1,  9:  PRINT  pax!:  LOCATE  10,  7:  PRINT  ylabl:  LOCATE  21,  9 
PRINT  yain! 

LOCATE  22.  14:  PRINT  xmin!:  LOCATE  22,  43:  PRINT  xlab$:  LOCATE  22,  71 
PRINT  xaax! 

VIEW  PRINT  24  TO  25: 

PRINT  SPC(14):  ,  "hit  any  key  to  return  to  Main  Menu" 

CIS 

LINE  (1,  0)-(0,  0),  ,  ,  stylo*/. 

FOR  i  «  1  TO  ntot4/, 
x  *  xdat(i) 
y  *  ydat(i) 

LINE  -<x,  y) 

NEXT  i 

DO:  LOOP  WHILE  INKEY$  =  «" 

CLS 

SCREEN  0:  COLOR  3,  0 
END  SUB 

SUB  PRUN  (dvm*,  rho!,  dt$,  FileName*,  if  lag*/,,  Temp!) 

********************************************************************* 

*  A  scanivalve  driver  program  for  the  HP  3852A  data  acquisition  system. 

*  Programmed  by  Lt  Galassi  6/21/89 

*  This  subroutine  performs  a  pressure  scan  oi;  the  turbine  blade 

’  using  the  HP3852A  data  acquisition  system  to  drive  a  scanivalve 
'  step  motor  over  30  pressure  ports. 
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The  ports  »  <  assigned  as  follows: 


0 

atmospheric 

1-10 

blade  static 

11 

atmospherric 

12-21 

blade  static 

22 

atmospheric 

23,24,25 

blade  static 

26 

atmospheric 

27 

pitot  head 

28 

pitot  static 

29-36 

atmospheric 

SHARED  chord()  AS  SINGLE 
DEFINT  I-H 

DIM  xdat(23)  AS  SINGLE,  ydat(23)  AS  SINGLE 
DIM  pscan(30)  AS  SINGLE 
CLS  2 

IF  ifJ^gX  «  1  THEN 

INPUT  "File  name  to  store  Cp  data  ?  ",  FileNamo$ 

END  IF 

iflag'/,  =  1 

OPEN  "A:\CP\"  +  "CP"  +  FileNameS  +  ".cp"  FOR  OUTPUT  AS  #1 
OPEN  "A:\CPV  +  "GR"  +  FileName*  +  ".DAT"  FOR  OUTPUT  AS  #3 
CLS 

'CALL  TM0N(dvm'/,) 

SLEEP  4 
COLOR  0,  3 
LOCATE  7,  7 
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INPUT  "Secondary  flow  on  (y)  or  (n)?",  ans$ 

IF  ans$  ■  "y"  THEN 
CLS  :  LOCATE  7,  7 
CALL  SECPRES (dvmy, ,  sfa!) 

PRINT  "Secondary  flow  pressure  (kPa)  ■  ",  sfa! 

END  IF 

LOCATE  8,7:  PRINT  "BLADE  PRESSURE  SURVEY  (approx  115  sec.  duration)" 

> ***•»»«*«******************************»*******•+*»******+**+******** 
'set  up  the  HP44701  for  voltage  measurements 

CALL  IBWRT(dvm’/.,  "RST  600") 

WRT$  *  "REAL  A ( 1 ) ,B(l) ,  PSCAN(30) ,0UTPUT(30) ,C(49) ,L,H,M,S" 

CALL  IBWRTCdvm'/.,  HRT$) 

'  Arrays  are  designated  as  follows: 

I 

’  A  -  voltages  from  the  pressure  transducer  calibration 

’  B  -  corresponding  pressures  to  the  voltages  listed  in 

’  array  A . 

’  PSCAN  -  actual  pressures  in  kPa  units 

’  OUTPUT  -  a  temporary  holding  array  for  HP3852  buff or. 

CALL  IBWRT(dvm'/.,  "VWRITE  A ,-10.000,10. 000") 

CALL  IBWRT(dvm'/.,  "VHRITE  B, -6. 8242, 9. 8271") 

’  Perform  the  actual  pressure  scan 

CALL  IBWRTCdvm'/,,  "USE  600;AZER0  ONCE;  NPLC  16") 

CALL  IBWRT(dvm'/,,  "CLOSE  101")  'this  call  homes  the  scanivalve 

CALL  IBWRTCdvm'/,,  "OPEN  101;CONF  DCV") 

CALL  IBWRTCdvm'/,,  "NRDGS  50") 

WRT$  =  "HEAS  DCV,  221, INTO  C"  ’read  atmospheric  press. 

CALL  IBWRTCdvm'/,,  WRT$) 
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CALL  IBWRTCdvm’/.,  "ST AT  L,H,M,S,C") 

CALL  IBWRTCdvn’/,,  "PSCAN(O)  «  K") 

CALL  IBWRTCdvmy,,  "DISP  PSCAN(O)") 

CALL  IBWRTCdvmy,,  "SUB  LOOP")  'HP  3852A  subroutine 

CALL  IBWRT(dvny.,  "INTEGER  I") 

CALL  IBWRTCdvay,,  "FOR  I  «  1  TO  30") 

CALL  IBWRTCdvm'/,,  "CLOSE  100") 

CALL  IBWRTCdvmy,,  "OPEN  100;C0NF  DCV;NRDGS  50") 

CALL  IBWRTCdvay.,  "MEAS  DCV,  221,  USE  600,  INTO  C") 

CALL  IBWRTCdvmy.,  "STAT  L.H.H.S.C") 

CALL  IBWRTCdvmy.,  "PSCANCD  =  M") 

CALL  IBWRTCdvm*/.,  "DISP  PSCANCD") 

CALL  IBWRTCdvm'/..  "NEXT  I") 

CALL  IBWRTCdvm*/.,  "SUBEND") 

CALL  IBWRTCdvm*/.,  "CALL  LOOP") 

SLEEP  119  'allow  120  seconds  for  scan 

WRT$  =  "CONV  A.B.PSCAN,  INTO  OUTPUT" 

CALL  IBWRTCdvm*/.,  WRT$) 

CALL  IBWRTCdvm*/,,  "VREAD  OUTPUT") 

>******************************************************************** 
'Read  output  into  a  temporary  file  on  the  hard  disk 
CALL  IBRDFCdvm'/,,  "C:\LELLO\PRESS.DAT") 

OPEN  "C:\LELLO\PRESS.DAT"  FOR  INPUT  AS  #2 
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FOR  i  «  0  TO  30 
INPUT  12,  pscan(i) 
NEXT  i 

CLOSE  12 


*  Output  the  data  into  FileNamet 

’  Array  locations  ar«  assignad  as  follows: 

1  0,11,22,29-36  -  Atmospharic  prassura 
1  1-10,12-21,23-25  -  Blada  static  prassuras 

*  27  -Pitot  head  prass. 

’  28-  Pitot  static  press. 

CLS 

COLOR  0 

PRINT  "Prassura  Measurements  on  Turbine  Blade  #3" 

PRINT  II,  "Presure  Measurements  on  Turbine  Blade  I  3" 

PRINT  :  PRINT  "Date  of  run:",  dt$ 

PRINT  :  PRINT  #1,  "Date  of  run:",  dt$ 

PRINT  :  PRINT  "Time  of  run:",  TIMER 

PRINT  :  PRINT  II,  "Time  of  run:",  TIMER 

PRINT  :  PRINT  II,  "Secondary  flow  pressure:",  sfa! 

PRINT  :  PRINT  II,  "Head  and  static  pitot  press.  (Kpa):",  pscan(27),  pscan(28) 

*  Note:  q  is  in  kPa 

q!  =  pscan(27)  -  pscan(28) 
deltah!  *  q!  *  1000!  /  249! 

PRINT  :  PRINT  "The  difference  in  height  of  the  manometer  columns" 

PRINT  "is  computed  by  scanivalve  measurements  to  be:  " ;  deltah! 
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INPUT  "  Is  this  correct?  (y  or  n)  ",  ans$ 

IF  ans$  <>  "y"  THEN 
BEEP 

PRINT  :  PRINT  "SCANIVALVE  readings  are  faulty" 
PRINT  "Returning  to  Main  Henu":  SLEEP  10 
CLOSE  «1:  CLOSE  »3:  COLOR  4,  11s  RETURN 
END  IF 


Vel!  *  SQR(2000!  •  q!  /  rho!) 

PRINT  :  PRINT  "Freestream  velocity  (m/soc):  ";  Vel! 
nu!  *  1 . 144E-05  +  (Temp!  +  23. IS)  *  9.48E-06  /  100! 

Re!  *  Vel!  *  .114  /  nu!  '0.144  *  chord  length 
PRINT  "Reynold's  no.  «  ",  Re! 

PRINT  :  PRINT  «1,  "Freestream  velocity  (m/soc):",  Vel! 
PRINT  #1,  "Reynold's  number  :",  Re! 

SLEEP  S 

’  Perform  shift  on  PSCANO 

FOR  i  =  11  TO  20 

pscan(i)  =  pscan(i  +  1) 

NEXT  i 

FOR  i  =  21  TO  23 

pscan(i)  *  pscan(i  +  2) 

NEXT  i 

PRINT  STRINGS (70,  "_") 

PRINT  #1,  STRING$(70,  "_") 

PRINT  #1,  "  Port  No.  CHORD  LOCATION  Cp" 

LOCATE  8,  1 
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PRINT  "  CHORD  LOCATION",  SPACES (12),  "Cp" 

LOCATE  9,  1 
PRINT  STRINGS (80, 

LOCATE  22,  1 
PRINT  STRINGS (80, 

VIEW  PRINT  8  TO  20 
CLS  2 

FOR  i  *  1  TO  10 

PRINT  i,  chord(i),  (pscan(i)  -  pscan(28))  /  q! 
PRINT  II,  i,  chord(i),  (pscan(i)  -  pscan(28))  /  q! 
PRINT  13,  i,  chord(i),  (pscan(i)  -  pscan(28))  /  q! 
NEXT  i 

PRINT  "  Hit  any  key  to  continue" 

DO:  LOOP  WHILE  INKEYS  =  "" 

FOR  i  =  11  TO  20 

PRINT  i,  chord(i),  (pscan(i)  -  pscan(2B))  /  q! 
PRINT  il,  i,  chord(i),  (pscan(i)  -  pscan(28))  /  q! 

PRINT  13,  i,  chord(i),  (pscan(i)  -  pscan(28))  /  q! 

NEXT  i 

PRINT  "  Hit  any  key  to  continue" 

DO:  LOOP  WHILE  INKEYS  =  "" 

FOR  i  =  21  TO  23 

PRINT  i,  chord(i),  (pscan(i)  -  pscan(28))  /  q! 
PRINT  II,  i,  chord(i),  (pscan(i)  -  pscan(28))  /  q! 

PRINT  #3,  i,  chord(i),  (pscan(i)  -  pscan(28))  /  q! 

NEXT  i 

PRINT  "Hit  any  key  to  continue" 

DO:  LOOP  WHILE  INKEYS  =  "" 

’This  section  drives  the  plotter. 
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ntotX  *  23 


'This  no.  corresponds  to  the  total  no.  of  suction  sido  ports 


ylab$  *  "Cp" 
xlab$  «  "x/c" 

pltlabt  «  "Cp  on  the  turbine  blade" 

xmin!  *  0! 

xnax!  *  1! 

yrain!  *  -4.5 

yraax!  *  4.5 

FOR  i  *  1  TO  ntoty, 

ydat!(i)  »  (pscan(i)  -  racai/28))  /  q 
xdat!(i)  *  chord(i) 

NEXT  i 

CALL  PLOTT(xdat! () ,  ydat!(),  xmin!,  xmax!,  ymin!,  ymax!,  xlab$, 
ylab$,  pltlabJ,  ntot*/,) 

CLOSE  II:  CLOSE  #3 

CALL  IBWRT(dvm‘/.,  "RST") 

END  SUB 

DEFSNG  I-H 

SUB  SECPRES  (dvra'/,,  sfa!) 

’Measure  the  secondary  flow  manifold  pressure 
CALL  IBWRT(dvm'/, ,  "REAL  Pl(l)  ,V1(1)  ,V2(1)  ,0UTPUT(i)") 

CALL  IBWRT(dvm'/, ,  "vwrite  VI,  .000,  .319") 

CALL  IBWRT(dvm*/,,  "vwrite  PI,  0.000,633.6505") 

CALL  IBWRT(dvra’/.,  "RST  600;USE  600;  AZERO  0NCE;NPLC  16") 

CALL  IBWRT(dva'/, ,  "CONFMEAS  DCV, 222, INTO  V2") 
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CALL  IBWRTCdvn*/.,  "DISP  V2(0)") 


CALL  IBWRT<dVRi7.,  "CONV  VI ,P1 .V2.INT0  OUTPUT") 

CALL  IBWRTCdvn*.  "VREAD  OUTPUT(O)") 
rd$  -  SPACES (16) 

CALL  ibrd(dv»X,  rdS) 

*f&!  «  VAL(rdS) 

CLS  :  LOCATE  7,  7 

PRINT  "Secondary  pressure  (kPa)  sfa! 

PRINT  :  LOCATE  9,  7 

PRINT  "Proas  any  key  to  continuo" 

DO:  LOOP  WHILE  INKEYS  -  "" 

CLS 

END  SUB 

SUB  TMON  (dvm*/.) 

••******************•*******************•****«*****************'*«*+* 

’Configure  tho  HP3852  to  sense  blade  temperature 
CALL  IBWRT(dvmy.,  "USE  600;RST  600;REAL  TjAZERO  ONCE") 

CALL  IBWRT(dvmy,,  "CLOSE  104") 

CALL  IBWRT(dvmy,,  "CONFHEAS  TEMPJ.206,  INTO  T") 

CALL  IBWRT(dvoy, ,  "OPEN  104") 

*********************** ♦***♦****♦♦*♦ ******************** +************ 

’Read  datum  from  HP3852 
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CALL  IBWRT(dvm'/, ,  "DISP  T") 


CALL  IBWRT(dvn7,,  "VREAD  T") 
rd$  «  SPACES (16) 

CALL  ibrd(dvm*/i,  rd$) 


T?  ■  VAL(rdS)  'convert  the  string  value  into  a 

’ number 
CLS 

LOCATE  7,  7 

PRINT  "Blade  Temperature  (F)  is  currently:  ",  (T!  •  9!  /  5!)  +  32! 
PRINT  :  PRINT 
IF  T!  >  65.56  THEN 

CALL  IBWRTCdvm’/.,  "CLOSE  104")  'this  call  turns  power 

BEEP:  SLEEP  1:  BEEP:  SLEEP  1:  BEEP 

PRINT  "Blade  temperature  too  high!":  PRINT 

PRINT  "Check  Power  Supply":  PRINT 

PRINT  "  Press  any  key  to  continuo" 

DO:  LOOP  WHILE  INKEYS  *  "" 

CALL  IBWRTCdvm*/. ,  "open  104") 

CLS  2 
END  IF 

LOCATE  8,  7 

PRINT  "Press  any  key  to  continue" 

DO:  LOOP  WHILE  INKEYS  *  "" 

CLS 

END  SUB 

SUB  TRUN  (dvmy, ,  Temp!,  FileNameS,  rho!) 
»**+**********************+*****+*♦*+*****♦**♦*♦****♦*♦****♦♦♦♦****♦* 

’This  subroutine  performs  the  temperature  scan  on  the  turbine  blade 
'and  stores  the  data  in  A:/heat/FileName$.dat.  The  HP3852  is  used  to 
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’acquire  data  along  with  the  HP  44713  FET  Multiplexer. 

DIM  nu(28)  AS  SINGLE 

DIM  rhol(28)  AS  SINGLE 

DIM  t2«can(28)  AS  SINGLE 

DIM  ar«a(23)  AS  SINGLE 

DIM  Vel(23)  AS  SINGLE 

DIM  S(23)  AS  SINGLE 

DIM  Nalt(23)  AS  SINGLE 

DIM  St(23)  AS  SINGLE 

DIM  xdat! (23) ,  ydat!(23) 

SHARED  chordv)  AS  SINGLE 
DEFINT  I-M 
CALL  TMONCdvn*/.) 


’Read  the  blade  temperaturos 


CALL  IBWRT(dvra'/,, 
CALL  IBWRKdvrnV., 
CALL  IBWRT(dvn'/,, 
CALL  IBWRT(dvmX, 
CALL  IBWRT(dvoy,, 
CALL  IBWRT(dvnX, 
CALL  IBWRTCdvoy,, 
CALL  IBWRT(dvm*/,, 


"RST  600;  AZERO  0NCE;NPLC  1") 

"SUB  TBLADE") 

"REAL  A (21)") 

"CLOSE  104") 

"CONFMEAS  TEMP J, 200-220, 223, INTO  A") 
"OPEN  104") 

"VREAD  A") 

"SUBEND") 


* 
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I 


*  Interactive  input  for  tho  power  paraaeters  on  tho  turbino  blade. 
DO 

CLS  :  LOCATE  7,  7 

INPUT  "Input  the  delta  H  of  the  tunnel:  ",  deltah! 

Vinf!  -  SqR(2!  /  rho!  •  deltah!  •  249!) 

LOCATE  9,  7 

PRINT  "  The  freestrea  velocity  is  Vinf!;  "  n/s" 

SLEEP  3 

CLS  :  LOCATE  7.  7 

INPUT  "Power  supply  current  (dof.  ■  38.5  aaps)  ?  ",  Curr! 

IF  Curr!  «  0!  THEN  Curr!  *  38.5 
LOCATE  8,  7 

Resist!  *  .13  'Ohms;  Blade  area  *  0.26192  m  x  O.llOOO  a 
7  INPUT  "New  filename  y  or  n  ?  ",  ans$ 

IF  ans$  *  "y"  THEN 

INPUT  "Enter  New  file  name  FileName! 

END  IF 

CLS  :  LOCATE  7,  7 

INPUT  "Is  the  secondary  injection  on  (y  or  n)  ?  ",  answer! 

IF  answer!  «  "y"  THEN 
CALL  SECPRES(dvay,,  sfa!) 

ELSE 

sfa!  =  0! 
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END  IF 


'  Compute  the  total  energy  flux  input  into  the  blade 
Qheat!  ■  (Curr!  "  2)  *  Resittt! 


CLS  :  LOCATE  3,  7:  PRINT  "Temperature  Scan  Menu" 

PRINT  STRINGS (80,  PR 'NT 

PRINT  "i)  Temperature  Scan" 

PRINT 

PRINT  "2)  EXIT" 

PRINT  :  PRINT  "type  your  selection  (1-2)" 
chsS  «  INPUTS(l) 

SELECT  CASE  chsS 
CASE  "1" 

CLS  :  LOCATE  7,  7 

'Turbulence  ON  temperature  scan  on  the  blado 
PRINT  "Press  any  key  to  initiate  scan" 

DO:  LOOP  WHILE  INKEYS  -  "" 

CALL  IBWRTCdvoy.,  "DISP  OFF") 

CALL  IBWRTCdvoy,,  "CALL  TBLADE") 

CALL  IBRDFCdvoy,,  "C:\LELL0\T2SCAN.DAT") 

CALL  IBWRTCdvoy,,  "DISP  ON") 

» ***************************************** **************** *********** * 
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•Read  data  from  storage 


OPEN  "C:\LELI.0\T2SCAN. DAT"  FOR  INPUT  AS  #1 

INPUT  *1.  t2ican(28) 

INPUT  II,  tinf! 

FOR  i  «  23  TO  13  STEP  -1 
INPUT  «1,  t2scan(i) 

NEXT  i 

FOR  i  *  12  TO  9  STEP  -1  'These  nodal  ara  daad 
t2scan(i)  ■  tinf! 

NEXT  i 

t2scan(0)  *  tinf! 
t2scan(l)  ■  tinf! 
t2scan(2j  *  tinf! 
t2scan(24)  «  tinf! 


FOR  i  «  8  TO  3  STEP  -1 
INPUT  II,  t2scan(i) 

NEXT  i 

INPUT  *1,  t2scan(2S) 

INPUT  II,  t2scan(26) 

INPUT  II,  t2scan(27) 

CLOSE  II 

•*»********************+ ******************************** **** ********** 

'Read  in  the  Prassura  Coefficients  on  the  blade  to  be  used  in  computation  of 
'Nusselt  and  Stanton  numbers. 

OPEN  "a:\cp\gr"  +  FileName$  +  ".dat"  FOR  INPUT  AS  II 

Vel(O)  *  0! 

FOR  i  *  1  TO  23 

INPUT  #1,  dummy,  dummy,  Vel(i) 
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NEXT  i 


CLOSE  «i 

'read  in  thermocouple  distances  along  the  blade  S() 

OPEN  "C:\lello\data\S.dat"  FOR  INPUT  AS  II 

FOR  i  -  i  TO  23 
INPUT  II,  S(i) 

NEXT  i 

CLOSE  *1 

’perform  energy  balance  calculations 
’conduction  constants 

uk!  *  .026  ’urethane  foam  thermal  conductivity 
ak!  *  ,0263  ’air  thermal  conductivity 
’convection  constants 
Cp!  «  1000.7 
Pr!  *  .707 

’compute  nu(i)  and  rhol(i)  on  the  blade 
FOR  i  =  0  TO  28 

nu(i)  *  1. 144E-05  +  ((t2scan(i)  +  tinf!)  /  2!  +  23.15) 

*  9.48E-06  /  100! 

rhol(i)  =  1.3947  -  .3797  *  ((t2scan(i)  +  tinf!)  /  2!  + 
23.15)  /  100! 

NEXT  i 

’radiation  constants 

eps!  *  .17  'emmisivity  for  polished  stainless  steel  at  300  K 

sigma!  =  5.67E-08 
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’compute  local  velocities  on  the  blade 


FOR  i  *  1  TO  23 

Vel(i)  «  SQR((Vinf!  ‘  2)  *  ABS(1!  -  Vel(i))) 

NEXT  i 

’compute  local  Reynold’s  numbers  on  the  blade 

FOR  i  «  1  TO  23 

Vel(i)  «  Vel(i)  *  S(i)  /  nu(i) 

NEXT  i 

’compute  the  Nurielt  and  Stanton  numbers  for  the  blade  only  in  rogion 
’of  active  thermocouples. 

ATot!  «  .0288112 

FOR  i  *  0  TO  23  ’clean  the  arrays 

Nslt(i)  «  0! 

St! (i)  *  0! 

NEXT  i 

dl!  *  .009 

FOR  i  *=  3  TO  4 

FI!  *  uk!  *  (t2scan(i)  -  t2scan(25))  /  dl! 

F2!  *  eps!  *  sigma!  *  (t2scan(i)  "  4  -  tinf!  4) 

F3!  *  (t2scan(i)  -  tinf!) 

h!  =  (Qhoat !  /  ATot!  +  FI!  -  F2!)  /  F3! 

Nslt(i)  *  ABS<h! )  *  S(i)  /  ak! 

St(i)  =  Nslt(i)  /  (Pr!  »  Vel(i)) 

NEXT  i 


i  *  5 

FI!  *  uk!  *  (t2scan(i)  -  t2scan(26))  /  dl! 

F2!  *  eps!  *  sigma!  *  (t2scan(i)  4  -  tinf!  “  4) 

F3!  *  (t2scan(i)  -  tinf!) 

h!  *  (Qheat !  /  ATot!  +  FI!  -  F2!)  /  F3! 

Nslt(i)  «  ABS Ch ! )  *  S(i)  /  ak! 

St(i)  "  Nslt(i)  /  (Pr!  *  Vel(i)) 

FOR  i  =  20  TO  23 

FI!  =  uk!  *  (t2scan(i)  -  t2scan(25))  /  dl! 
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F2!  *  eps!  ♦  sigma!  •  (t2scan(i)  “  4  -  tinf!  *  4) 

F3!  *  (t2scan(i)  -  tinf!) 

h!  «  (Qhaat!  /  ATot!  +  FI!  -  F2!)  /  F3! 

Nslt(i)  «  ABS(h! )  *  S(i)  /  ak! 

St(i)  «  Nslt(i)  /  (Pr!  *  Vel(i)) 

NEXT  i 

dl!  *  .012 

FOR  i  «  17  TO  19 

FI!  *  uk!  *  (t2scan(i)  -  t2scan(26))  /  dl! 

F2!  ■  ops!  *  sigma!  *  (t2scan(i)  “  4  -  tinf!  “  4) 

F3!  *  (t2scan(i)  -  tinf!) 

h!  *  (qhaat!  /  ATot!  +  FI!  -  F2!)  /  F3! 

Nslt(i)  *  ABS(h!)  ♦  S(i)  /  ak! 

St(i)  =  Nslt(i)  /  (Pr!  ♦  Vel(i)) 

NEXT  i 

dl!  =  .005 

i  =  6 

FI!  =  uk!  *  (t2scan(i)  -  t2scan(27))  /  dl! 

F2!  =  eps!  ♦  sigma!  *  (t2scan(i)  ‘  4  -  tinf!  “  4) 

F3!  =  (t2scan(i)  -  tinf!) 

h!  =  (Qhaat!  /  ATot!  +  FI!  -  F2!)  /  F3! 

Nslt(i)  =  ABS(h!)  *  S(i)  /  ak! 

St(i)  =  Nslt(i)  /  (Pr!  *  Vel(i)) 

dl!  =  .01 

FOR  i  =  15  TO  16 

FI!  =  uk!  *  (t2scan(i)  -  t2scan(27))  /  dl! 

F2!  =  eps!  *  sigma!  *  (t2scan(i)  “  4  -  tinf!  *  4) 

F3!  =  (t2scan(i)  -  tinf!) 

h!  =  (Qhaat!  /  ATot!  +  FI!  -  F2!)  /  F3! 

Nslt(i)  =  ABS (h ! )  »  S(i)  /  ak! 

St(i)  *  Nslt(i)  /  (Pr!  *  Vel(i)) 

NEXT  i 

dl!  =  .0045 

FOR  i  =  7  TO  8 

FI!  =  uk!  *  (t2scan(i)  -  t2scan(28))  /  dl! 

F2!  =  eps!  *  sigma!  *  (t2scan(i)  "  4  -  tinf!  "  4) 
F3!  =  (t2scan(i)  -  tinf!) 
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h!  «  (Qheat!  /  ATot!  +  FI!  -  F2!)  /  F3! 

Hslt(i)  *  ABS(h! )  *  S(i)  /  ak! 

St(i)  =  JJslt(i)  /  (Pr!  ♦  Vel(i)) 

NEXT  i 

dl!  «  .005 

FOR  i  «  13  TO  14 

FI!  *  uk!  *  (t2scan(i)  -  t2scan(28))  /  dl! 

F2!  *  eps!  *  sigma!  *  (t2scan(i)  "  4  -  tinf!  “  4) 

F3!  «  (t2scan(i)  -  tinf!) 

h!  *  (qheat!  /  ATot!  +  FI!  -  F2!)  /  F3! 

Nslt(i)  «  ABS(h! )  •  S(i)  /  ak! 

St(i)  *  Nslt(i)  /  (Pr!  *  Vel(i)) 

NEXT  i 

'save  data  on  disk 

CLS  :  LOCATE  1,  10:  PRINT  "Temperature  output" 

PRINT 

OPEN  "A:\HEAT\"  +  "HT"  +  FileNamol  ♦  ".DAT"  FOR  OUTPUT  AS  #2 
PRINT  #2,  "File:";  SPC(l);  FileNameS 

PRINT  "Secondary  flow  pressure  =  SPC(l);  sfa! ;  SPC(2);  "kPa" 

PRINT  «2,  "Secondary  flow  pressure  =  SPC(i);  sfa!;  SPC(2);  "kPa" 
PRINT  #2,  "Electrical  energy  input  =  ";  SPC(l);  Qheat!;  SPC(l);  "Watts" 
PRINT  ,  "x/c",  "Re  no.",  "Nusselt",  "Stanton" 

PRINT  12,  "x/c",  "Re  no.",  "Nusselt",  "Stanton" 

FOR  i  =  1  TO  23 

PRINT  #2,  chord(i) ,  Vel(i),  Nslt(i),  St(i) 

PRINT  chord(i) ,  Vel(i),  Nslt(i),  St(i) 

NEXT  i 

CLOSE  #2 

PRINT  "Press  any  key  to  continue" 
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DO:  LOOP  WHILE  INKEY$  *  "" 

'♦*********♦******•***•******♦****♦*•**•****♦***♦***♦•*♦**♦*♦♦♦♦♦♦♦♦♦* 

'Plot  out  the  heat  transfer  solution  for  the  blade  using  the 
'PLOTT  subroutine. 

'This  section  drives  the  plotter  for  turbulence  ON  conditions. 

ntotlX  «  8  'This  no.  corresponds  to  the  total  no.  of  pressure  side 
'thermocouples 

ntot2‘/,  «  12  'This  no.  corresponds  to  the  total  no.  of  suction  side 
'thermocouples 

ylab$  «  "Temp.  (C)" 

xlab$  *  "x/c" 

pltlaM  =  "Temp,  on  turbir.o  blade  " 

xmin!  *  0! 

xmax!  =  1! 

ymin!  =  20 

ymax!  =  65 

FOR  i  =  1  TO  ntotiy. 
ydat(i)  *  t2scan(i) 
xdat(i)  =  chord(i) 

NEXT  i 

FOR  i  *  ntotlX  ♦  1  TO  ntot2’/. 
xdat(i)  «  chord(i) 
ydat(i)  »  0 
NEXT  i 

FOR  i  *  ntot2X  +  1  TO  ntotl'/,  +  ntot2'/,  +  3 
ydat(i)  *  t2scan(i  -  4) 
xdat(i)  *  chord(i) 

NEXT  i 

ntot'/,  =  ntotl*/,  +  ntot2'/,  +  3 
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CALL  PLOTT(xdat! () ,  ydat!(),  xraini,  xnax!,  ymin!,  ymax!,  xlab$, 
ylabS,  pltlab$,  ntot*/,) 

CASE  "2":  EXIT  DO 

CASE  ELSE:  BEEP 

END  SELECT 

LOOP 

CALL  IBWRTCdvmX,  "RST") 

END  SUB 
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